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eee 


Zur Behandlung der Strémung durch einen Voith-Schneider-Propeller 
mit kleinem Fortschrittserad 


Von W.-H. Isay 


1. Einleitung. Wir betrachten die Strémung durch einen rotierenden Voith-Schneider-Propeller 
mit sechs Fliigeln. Die Strémung sei eine inkompressible Potentialstromung, und wir bauen das 
Strémungsfeld auf, indem wir die Fliigelprofile nach der in der Schaufelgittertheorie bewahrten 
Weise mit Zirkulationsverteilungen belegen. Der Propeller liege in einer konstanten Anstrémung 
mit der Absolutgeschwindigkeit u, in Richtung der positiven x-Achse. Da es sich um ein instatio- 
nares Problem handelt, das sich auch nicht quasistationar machen laBt, andern sich die auf den ein- 
zelnen Propellerfliigeln angeordneten Zirkulationsverteilungen dauernd, und durch die zeitliche 
Anderung der festen Wirbelelemente werden freie Wirbel induziert, die mit der Strémung in Rich- 
tung der positiven x-Achse fortgetragen werden!. Das (absolute) Strémungsfeld in einem bestimmten 
Punkt (x = r cosy, y =r sing) der Ebene besteht also neben der konstanten Strémung u, aus den 
beiden von den gebundenen und den freien Wirbeln induzierten Geschwindigkeitsanteilen. Wir 
werden beide Anteile getrennt berechnen. 


2. Das Geschwindigkeitsfeld der gebundenen Wirbel. Zur Berechnung des von den gebundenen 
Wirbeln auf den Propellerfliigeln induzierten Geschwindigkeitsfeldes kénnen wir voraussetzen, dal 
die Konturkurve des Propellerfliigels angenahert ee AP 
mit der Kreisbahn des Propellerradius R zusam- a = 
_ menfalle (Abb. 1), jedenfalls insoweit, als wir die X 
Zirkulationsverteilungen der Propellerfliigel auf ie x 
dem Kreisumfang anbringen kénnen. y) 

Die sechs Konturkurven der Fligelprofile sind / 


dann in der Form | 
: te ; 4 
1 (v+ ra m+ ow é 


Y 


“2 
Cz= Re \ 
G1... 0; aS y=) \e 
gegeben. Das gilt fiir jeden beliebigen Zeitpunkt. \ 


Dann ist das von den auf den sechs Propeller- 
fliigeln angebrachten Zirkulationsverteilungen p(y, ¢) 
im Punkte z =r e'” induzierte Geschwindigkeits- 


feld gegeben durch 


‘ i 5; ( gal 7 R dy ‘ (1) 
Wye Uy 5 fr(v.s } 3 _ w (v4 phat ot) 
galery pee e | : 


Dabei haben wir sofort beriicksichtigt, daB ja, wenn y,(y,t) die Zirkulationsverteilung auf dem 


q-ten Fliigel ist, offenbar gilt 
q—l1a2\ q—1lia 
YolPot) = 11 (pst rec he y (vst oe | 


undenen Wirbeln ist hier die aus der instationaren Trag- 


i schei ischen frei d geb 
ee BW. Bimnbour Math. Mech. 4 (1924) S. 277; u. W. Schmeidler, 


fliigeltheorie bekannte. Vel. z. B. W. Birnbaum, Z. angew. 


Z. angew. Math. Mech. 19 (1939). 
26 
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Die Zerlegung von (1) in Real- und Imaginarteil liefert 


rsing —Rsin(y +2> "2 + ot) 


6 a 
i! we q—la 
ult, ~) = a5 Jo (ver = Rdy, 
: ra | oe RE 41? — 27 Roos (p—y—! x —ot) 
. (2) 
1 6 « ecti ge Reos (y +15" a +o) r cos Pp 
tops Sf r(ner te) AES, 
3 m2 | ree Rt +r? —2r R cos (g—y—15 *2—o) 


3. Das Geschwindigkeitsfeld der freien Wirbel. Wir betrachten ein festes an die Stelle y eines 
Propellerfliigels gebundenes Wirbelelement y(y,t). Durch die zeitliche Anderung von y lést sich in 
jedem Moment ein Wirbelelement —y(y,t)dt ab, das mit der konstanten Geschwindigkeit ug in 
Richtung der positiven x-Achse fortgetragen wird. Das feste Wirbelelement selbst beschreibt mit 
der Winkelgeschwindigkeit ~ eine Kreisbahn. Die Umlaufzeit betragt 2 z/w. Zur Zeit t befinde sich 
unser Wirbelelement y(y,t) an der Stelle @t-+ y der Kreisbahn. Dann befand sich y zur Zeit 
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Abb. 2. AbflieBende freie Wirbel eines Voith-Schneider-Propellers mit dem Fortschrittsgrad 1/16. (Idealisiert, d. h. ohne Beriicksichtigung 
der gegenseitigen Beeinflussung der freien Wirbel). 


t—V/w an der Stelle w@t+yw—@# der Kreisbahn. Dort hat sich ein freies Wirbelelement 
— y(p,t— B/w) dt von ihm abgelést, das sich mittlerweile an der Stelle 


Y= Rsin(w@t+y—V), X = Reos (t+ p—8) + uy— 


befindet. Ein solches freies Wirbelelement induziert an einem beliebigen festen Punkt z = r e'? die 
Geschwindigkeit 


: : oO 
) ee one 
i ? (v, @ dd ‘ dd 
mite — 
oO o o 
o 


2% ; Ax pee R Oh ta) ear 


Nun befinden sich als Folge unseres festen Wirbels y (an der Stelle w t + w) auf der ganzen Bahn bis 
nach X-> oo, d.h. auch #—> 00, lauter freie Wirbel, die alle zusammen an der Stelle z — re’? die 
Geschwindigkeit 


oo A Oo 
4 i pel? — Rie Olds Gane) eae Oo 
a) 
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induzieren. Alle festen Wirbel der sechs Propellerfliigel liefern dann ein abgehendes Wirbelband 
von freien Wirbeln, das die Geschwindigkeit 


7% 


& Ly pean q—l Tee 
ees i R y | Hite ota 
faa fame ge at ee = 


spine ed lat JRE (ote re pas) —- yy) 
induziert. Das ist die strenge Formel, wenn man von der schwer zu erfassenden Zusatzbewegung 
der freien Wirbel infolge ihrer gegenseitigen Beeinflussung absieht. Es hat sich nun aber gezeigt, 
da® einer theoretischen Verwendung dieser Formel groBe Schwierigkeiten entgegenstehen. Diese 
Schwierigkeiten sind dadurch begriindet, da®B die Kurven, auf denen sich die abflieBenden freien 
Wirbel bewegen, zykloidenférmig sind. Dieser Umstand macht eine Auswertung der vorkommenden 
Integrale sehr schwer. 

Man kann diese Schwierigkeit umgehen, wenn man sich auf die Behandlung von Propellern 
mit einem Fortschrittsgrad beschrankt, der kleiner als etwa 1/,, ist’, (Forischrittsgrad ist das Ver- 
haltnis u)/Rw.) Dann werden namlich die zykloidenférmigen Kurven, auf denen sich die abflieBen- 
den freien Wirbel bewegen, fast Kreise. Man kann dann so tun, als ob sich alle wahrend eines Um- 
laufes des festen Wirbelelementes y(y,t) erzeugten freien Wirbelelemente auf einmal ablésten und 
zwar zur Zeit t. Dadurch entsteht zur Zeit t eine freie Kreiswirbelverteilung }(p,3/w) d/o 
OS 7 =< 27): 

Von demselben festen Wirbelelement haben sich natiirlich auch zu den Zeiten t— 2a n/w 
(n = 1, 2, ...) lauter freie Kreiswirbelverteilungen abgelést. Diese Wirbelringe kénnen in der hier 
betrachteten Naherung als nahezu unabhangig von y (d. h. von der Kennzeichnung des festen Wir- 
belelementes im Verband des Propellerfliigels) angesehen werden. Denn die Uberlagerung aller 
Wirbel des Propellerfliigels bewirkt ja lediglich eine geringfiigige Verdickung der freien Wirbelringe, 
die im Rahmen der hier gemachten Vernachlassigung nicht beriicksichtigt zu werden braucht. 


22 w 


. % 
Wir vollziehen also den Ubergang von 7(p,t) zu ['(t) = i) v(p,t) dy. (Zu diesen Betrachtungen 
vgl. Abb. 2.) sa 

Von einem Propellerfliigel stammen somit unendlich viele, freie Wirbelringe (Kreiswirbelver- 
teilungen), die sich zu den Zeiten t— 2 n/w von dem Propellerfliigel abgelést haben, und die zur 
Zeit t mit ihren Mittelpunkten an den Stellen uy 27 n/w auf der x-Achse (n = 0, 1, 2, ...) ange- 
langt sind. (Dabei wird wieder abgesehen von der Zusatzbewegung der freien Wirbel infolge ihrer 
gegenseitigen Beeinflussung.) Nun haben wir es ja nicht mit einem sondern mit sechs Propeller- 
fligeln zu tun, so daB wir eine dichte Folge von Kreiswirbelverteilungen haben, die zur Zeit t mit 
ihren Mittelpunkten an den Stellen uy 27 n/6 (n= 0, 1, 2,...) angekommen sind. Da wir einen 
Propeller mit kleinem Fortschrittsgrad behandeln, so ist Ax* = 22 u)/6 Kein. Somit kénnen 
wir im Rahmen der hier angesetzten Theorie den letzten konsequenten Schritt tun: Wir voll- 
zichen den Grenziibergang Ax*—> 0 und betrachten nur noch eine kontinuierliche Folge von ab- 
flieBenden, freien Kreiswirbelverteilungen der Starke — = lo (5) dO( 0 O22 7). 

Somit ergibt sich also endgiiltig das in einem beliebigen Punkt 2 = r e'? von den freien Wirbeln 
induzierte Geschwindigkeitsfeld zu 


Us— Lvs = — 


i / [ L'(O/a) dd dx* — (3) 


rev R eo =x 
nt) 
Mit dieser Formel werden wir in Zukunft arbeiten.. Die Zerlegung von (3) in Real- und Imaginarteil 


liefert 
o 2n 


1 ; I (9/w) (R sin 8 —r sin P) dd CE a ee 
uf(r, 7) Set are R2 +r? + x*?— 2r R cos (p — 9) —2rx* cosy + 2R x* cos 3’ 


ey (4) 
oO 20 : 
i! j I’ (9/w) (r cos p — R cos & — x*) dd dx* 
v,(7, ?) any R? +r? + x*2 — 2r R cos (vy — 9) —2 rx* cosy+2 Rx* cosd” 
2am +r? + (p ) 
0 0 


1 Méglicherweise ist die hier entwickelte Methode auch fiir groBere Fortschrittsgrade noch eine brauchbare 
Naherung. Ich hoffe, darauf in einer spateren Arbeit zurickzukommen. _ 

Zusatz bei der Korrektur: Diese Vermutung bestatigt sich, wie in einer welteren Arbeit gezeigt 
werden wird. Dort wird auch Formel (3) noch einmal in einer fiir beliebige Fortschrittsgrade gultigen 
Weise begriindet werden. 


On 
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4. Randbedingung und Aufstellung der Integralgleichung. Die Randbedingung fur die Strémung 
an den Propellerfliigeln lautet in bezug auf das ruhende Absolutsystem!: An jeder Stelle des Fligels 
soll die Normalkomponente der Absolutgeschwindigkeit tibereinstimmen mit der Geschwindigkeit 
des Fliigelelementes normal zu diesem selbst. Es muf somit 


w r cos (7— ®) = (v, + vs) cosy — (Uy + U, + Uy) Sin H 


sein lings der Fliigelprofilkontur r(®). Dabei ist mit 0 =y + wt 
ee r(®) sin B + r(P) cos © 
81] ~ +) cos ® — r(P) sin © 


der Anstieg der Profilkontur (Abb. 3). 

Ware die Profilkontur der Propellerfliigel genau ein Kreisbogen, so wiirde tg 7 = — ctg ® und 
infolgedessen cos (7 — ©) = Osein. Wir diirfen natiirlich mit der Vereinfachung nicht so weit gehen; 
es soll lediglich r(p,t) = R gesetzt werden, aber r’(g,t) ist von Null 
verschieden. Uber r’(y,t) ist noch folgendes zu sagen: Zunachst 
kann r’ als unabhangig von py angesehen werden (—a [<p Xa; 
p ist wie y jetzt die Winkelkoordinate zur Bezeichnung der ein- 
zelnen Punkte des Propellerfliigels). Langs eines Propellerfligels 
bleibt ja r’ praktisch konstant. Dagegen variiert r’ wahrend der 
Umdrehung des Propellerfliigels mit der Zeit t, da dieser ja 
wahrend einer Umdrehung von positiven Anstellwinkeln auf 
negative und wieder zuriick auf positive umgesteuert wird. Auf 
diese Weise kann je nach der Lage des Anstellwinkels Null ein 
} Schub in einer beliebigen Richtung erzeugt werden, und das ist 

ja einer der Vorztige des Voith-Schneider-Propellers. Somit ist 


, , —l1 
r= —rhc0s(mt+ 2a—a +4), 
poe Tro ==_konst.\ (piers 16). 


Dabei ist fiir den ersten Propellerfliigel p = 1, usw. zu setzen und A wird durch die Lage des An- 
stellwinkels Null bestimmt; soll dieser z. B. wie in Abb. 1 bei w t = 7/2 und bei w t = 3 77/2 liegen, 
so istA = 0 zusetzen. Damit lautet die Randbedingung fiir die Strémung am Propellerfliigel 1 
—w Rrocos(wt-+A) + uy R cos (py + w t)— uy 1 cos (wt +) sin (p+ ot) = 
= —[v,(R,p+ ot) +, (R, py +o 2B) [ro cos (wt +/) cos (p +t) + Rsn(p+o?)] (5) 
+ [uy (R,p + ot) + u, (R,p + 2)] [ro cos (wt + A) sin (p + t)— Reos(p+ort)]. 


Gleichung (5) stellt eine Integralgleichung fiir die Bestimmung der Zirkulationsverteilung y(y,t) dar. 


Sie lautet, wenn wir uy, vs, U,,, Vv, gemaB (2), (4) einsetzen, zur Abkiirzung — = y schreiben und 
die elementare Zwischenrechnung hier iibergehen, 


oo ln 
@ 79 cos ( t +A) — uy cos (p + @ t) + wy cos (0 t +A) sin (p + ot) + 7 — pei ‘ 
@ 
0 0 


/ 


x sin (p-+w t)+sin (p+ t—?)-4 “a cos (w t-++A) cos (p+ t)— F008 (ot tA) [l1—cos(y+wt—#)] 


xX 


1 di dy 
1 +> x° — cos (p + wt — #) — x cos (p + wt) + x cosd 


1 i S ‘ q—l1 It 
er as ps (Oba A) DS a = 


123) 
q=1 


eke ; qe Sines l a5 
=a2> | 7(vt+ 3 7) ez (7 2 Bs n) dy 


q=1 


—% 


* Vel. W.-H. Isay, Ing.-Archiy 22 (1954) S. 203. Dabei sind u und v die x- und y- Komponente der Absolut- 


geschwindigkeit c. 
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Gleichung (6) ist die Integralgleichung unseres Problemes; sie gilt fur jeden beliebigen Zeitpunkt. 
Es ist 


e 
TQ) =f y(y.2) dy (ass. 
athe y 
Gleichung (5) war die Randbedingung fiir den ersten Propellerfliigel. Die Randbedingung fiir den 
p-ten (p = 2, ... 6) entsteht aus (5), wenn man dort 
p—la 
t durch t + <a as 


ersetzt. Die aus ihr hervorgehende Integralgleichung wird durch die Funktion y wt eae = 


gelést, wie es ja auch sein muB. Deshalb brauchen wir uns um diese weiteren Randbedingungen 
nicht mehr zu kiimmern. 


5. Lésung der Integralgleichung. Fiir die Auflésung der Integralgleichung wollen wir eine Nahe- 
rungsmethode anwenden, die sich unter den verschiedenen Lésungsmethoden, die der Verfasser ge- 
priift hat, als die einzige fiir praktische Berechnungen gangbare erwiesen hat. Wir machen den 
Ansatz 


an(yp) eneas > C4) = a, r) 


: (7) 
A, ermot (AST p= Ave 


Wenn man mit dem Ansatz (7) in die Integralgleichung (6) hineingeht, benétigt man, wie leicht er- 
sichtlich, die Kenntnis einiger Integrale, die wir nachfolgend zusammenstellen. Dabei fihren wir 
noch die Abkiirzungen 

ri 


R 008 (Wt +A) = Alt), wot+p=@ (8) 


ein. Dann ist 


o 27 


1 inl be LEE OG BON ot 0 Eo Saeed) 0h Pe 


ae L +f — cos (@—9) —z cos @ + 7 005 9 
0 0 
1). i(n—1)® Be Me ee <3 
— (i + h(t) et ) fiir eS 5? (9) 


CY (i hle)) eto) = (i b(n) tm +00 4 NG + HQ oto ne 
| fir —[<@<F, (n2)) 


und entsprechend ergibt das Integral mit e—'"® den konjugiert komplexen Wert. Aus (9) folgt in 
elementarer Weise ferner 


o 20 


fiir n gerade: lal cos n } — db dy = 
0 6 


eee cos(n—1)@ _ fiir 5S oe 
= + fsin (n— 1) ® + Bin (n +1) B] —“O feos (n— 1) @—2 co (n + 1) O] (10) 
fiir —_2<O<-<, 
wo 2x 
fur n ungerade: eal | cos ni} — di dy =— ~sin (n—1) 6 + eos (n—1)@, 


0 06 
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co Qn 
fiir n gerade: = | | sin n = di dy = 2 cos (n— 1) B+ sin (n— 1) @, 
66 
co ln 
fiir n ungerade: = | | sin n o— di dy = 
60 (11) 
h(t). _ Tt 32 
= cos (n— 1) ® + ou sin(n—1)@_ fir z= =< a 


= —~ [eos (n— 1) ® + 2 cos (n + 1) 6) —" [sin (xn — 1) 6— 2 sin (n + 1) P| 
fir — >< C=. 


Um den Gedankengang der Auflésung unserer Integralgleichung (6) hier nicht zu unterbrechen, soll 
der Beweis der Relation (9) in einem besonderen Abschnitt 6 gebracht werden. 

Interessant ist, daB diese Integrale, durch die ja im wesentlichen der EinfluB der abflieBenden 
freien Wirbel charakterisiert wird, verschieden sind, je nachdem ob sich der betrachtete Propeller- 
fliigel im IT. und III. oder I. und [V. Quadranten befindet (Winkelkoordinate ©). Das liegt natir- 
lich daran, daB im letzteren Fall der Propellerfliigel durch die abflieBenden freien Wirbel hindurch- 
tritt, im ersteren aber nicht. 

Die Integralgleichung (6) gilt nun fiir alle Zeiten, d. h. fiir jedes ®. Zu ihrer Auflésung mussen 
wir deshalb fiir das Doppelintegral auf der linken Seite von (6) eine fiir alle ® giiltige Darstellung 
suchen. Man kénnte zunachst daran denken, fiir jedes einzelne der Integrale (9) bzw. (10), (11) eine 
fir 0< O< 2a giiltige Fourierentwicklung herzustellen. Es hat sich jedoch gezeigt, daB dieses 
praktisch nicht durchfiihrbar ist. Denn eine solche Fourierentwicklung diirfte, damit der Ansatz (7) 
zur Lésung der Integralgleichung ausreicht, in keinem Falle héhere Potenzen als e°'®' haben. An- 
dererseits tritt z. B. fiir das Integral in (10) grade auf Grund des Ansatzes (7) fiir n = 6 auf 

— sin ® fiir Lue 
(12) 
~sinSD+—sin7® fir —L a Oat. 


Dieser Ausdruck kann aber unméglich durch eine Fourierentwicklung approximiert werden, die 
nur bis sin 6 @ reicht. 

Es bleibt also nur die Méglichkeit, das Doppelintegral auf der linken Seite von (6) als ganzes 
[bei Verwendung des Ansatzes (7)] durch ein Fourierpolynom zu approximieren. Denn das Doppel- 
integral als ganzes hat natiirlich nicht die hohe Periodizitaét wie z. B. der Ausdruck (12), so daB man 
mit einem Fourierpolynom bis zum sechsten Grade auskommen diirfte. Es hat sich gezeigt, daB 
dieser Weg zum Erfolg fiihrt. Aber fiir diese Approximation ist jeweils die Kenntnis der Konstanten 
A,, aus (7) erforderlich, so daf$ unsere Methode in der Praxis auf ein Iterationsverfahren hinaus- 
lauft. 

Der Weg zur Auflésung der Integralgleichung (6) ist also folgender: In erster Naherung wird das 


Doppelintegral auf der linken Seite von (6) weggelassen und auf diese Weise die erste Naherung 
(1) (1) 
y(y,t) baw. L(t) bestimmt, die dann in das Doppelintegral eingesetzt wird; mit dieser linken Seite 


(2 2 
von (6) wird die zweite Naherung a t) bzw. T ‘1 bestimmt, diese wieder in das Doppelintegral ein- 
gesetzt und so weiter. Da die aufzulésenden Integralgleichungen fiir jede Iteration den gleichen 
Kern haben, kann das Iterationsverfahren weitgehend schematisiert werden; das ist fiir die nume- 
rische Rechnung angenehm. : 

Wir wollen nun die Durchfiihrung eines solchen Iterationsschrittes effektiv verfolgen; es sei der 
m-te Schritt. Durch Einsetzen der (m— 1)-ten Naherung (7) in das Doppelintegral gewinnen wir 
unter Verwendung der Formeln (9) bzw. (10), (11) und nach Durchfiihrung einer harmonischen 
Analyse mit der Schrittweite 7/6 folgende Darstellung: 

27 (m—1) 


iar ayn Niners orate pep Ih, c_,=€ 
fall | it a a dd Wim De Cy el? 1 h(t) > k,, eH? | H 7 ks (13) 


00 f=—6 p=—6 ky aa k,,| 
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Damit nimmt die Integralgleichung (6), wenn man den Ansatz (7) einsetzt, die Form an 


/ 


a a mee + 2 uy-3 ea + @ t) 


+ Se Op  ginloton) 4 — R cos (mt +A) > oy eite(p + ot) sth e > cos (@ t + A) x 


H=—6 u=—6 


sR Ronda k S fiir , 


n=—7 


a 
~ — y —— _o8 p—Yy—-— 
< p—y v = n e bis ; t “ae n 3 j 
(ug? OG rae a ae DS era eas Ay Tenor, 
Es ist nicht schwer zu zeigen, daB bei unserem Ansatz Ye 


s (m) q— la (m) 
P(t +132) Eee ety iat 4 CAnee 
q=1 
wird. Durch Koeffizientenvergleich der e'"* zerfallt (14) in 15 Integralgleichungen, von denen nur 
8 wesentlich sind; denn mit seul ist ja zugleich auch a_n(y) =4,(y) bekannt. Diese Integral- 
gleichungen sind fiir n = 0, 1,2,..., 7 von der Form 


a a ; . aA ote ae 2% 
(m) (m) re) i—n a - i—n ae =e 
yp 3 3 ; 
Sip) = el an(y) (oe? ee ee ee 
a (15) 
DEG IL 
2% ee es im ae ae 
+ ean eg ee a aes oe id El + ae > Djs ck ) dy 
2 9 ? 
mit In(P) = f(y): 
(m) ia 1 z (m) rf (m) a 6 (m) ro ) 
f(~) = “oR 3; OLE ig ied aa eae Rh OO erates aR De 
(m) ‘ - (m) 7 (m) / ie) aye Ae) 
f(y) = @ 16 e'* — uy e'? + cc, e'? + = e'4 ky + st SEY CR A ea ES + Clade 
m , m / (m , (m) 
Deore on hot pote 
i 
(m) my en 3 ae / we) 
JAG) a ¢,e%'? +5 — ek, e729 aS om Cah ere 
> (16) 
(m) (m) : (m) ; 
fo) = i, e ete + pale tlt e eet” +55 et k; em ; 
(m) (m) Doe (m) ; t ‘ (m) : 3 1 (m) } 
As(~) = ¢5 eP? + oe e4 k, ett? + ff Ces ben’ ae 2 Ag Ca. 
(m) (m) ri en) 
6(P) = Ce eff? IR ek; e'? 
(m) Ee ré (m) ™ 
f(y) = Ro’ - ebip ae — Ge el, | 


In der Praxis ist ” meist schon so klein, da es nicht mehr beriicksichtigt zu werden braucht. 
Wir haben nun die Auflésung der Integralgleichungen (15) durchzufiihren. Dazu nehmen wir 
die Transformation 


O Ss 
tgp =— tg 5 C08 s , ig 5 =— tg-5 coso (0s $s) (17) 
vor. Dabei wollen wir festsetzen, dai 
tgs S01 db. a< 12° (18) 


sei, und zwar grundsatzlich fiir die ganze hier dargelegte Theorie. Fir die Praxis bedeutet das kaum 
eine Einschrankung, denn ein sechsfliigeliger Voith-Schneider-Propeller, dessen Fliigel sich je tiber 
mehr als 24° erstrecken, diirfte selten vorkommen, 
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Mit der Transformation (17) werden nun die einzelnen Teile des Kernes von (15) umgerechnet. 
Dabei ist wegen (18) 


m 
1 + tg? a It 
und auch sonst ergeben sich Vereinfachungen durch (18). 
aes 
Es ist z. B. dy =2tg>sinado, 
OW eal 1 1 
Ce roe at cone Secon es 
tg o8 
bie je Fell tee (cos 6 — cos s) 
de ee as 
1 (cos 6 — cos s) 
ie 
2 
a vee 2 
ar & > 
Tt a 2 a 2 
a 1+ tg>te54 = (cos og — cos s) + tg? 5 }- = (cos g — cos s)? 
Aas tea tg 


Wir wollen diese Transformation auf die neuen Variablen s und g hier nicht weiter verfolgen, da die 
Rechnungen durchaus elementar sind. Wegen der Kleinheit von tg 4 kénnen Nenner in geome- | 
trische Reihen entwickelt und mit dem zweiten oder dritten Glied abgebrochen werden. Das Er- | 
gebnis ist: Wenn wir noch | 
(m) 1 &&) (m) 1 ™ 
InP) = Say Sal) » An(p) = Sag Inl2) 
setzen, gehen die Integralgleichungen (15) in solche von der Form 


24 


N N (n) 


(m) 1 f@™ 
8n(8) =— [ y«(o) | _— SS, b,, Sin fu s cos » 7 dg (19) 


cos 0 — COS $ 
0 


tiber, die aus der Theorie der Schaufelgitter wohlbekannt sind!. In unserem Falle ist | 


c— ii 


sin ze n sin. eae sl eu oe ] 
= — a —— ig) sin —n 
; 3 3 3 a 3 3 
b@) —— 2it te 3 | 
10 i e9 Rae Cars eet lar nl eonteR eae 2 
&% 83 eae hay Ss 
PARE Tt 
cos —n cos ——=-n cos —n 
n 1 n 3 3 m 3 
te | (1) 4 bose |b 2 6g ee 
: ae 
sin2 YA sin? oa tg? — sin? és 
cos—n co eu is 
: 1 ay mow fe i eae 
be) =—2 tg? 2\ 9 ( 1s ie | ae 3 tg* 5) pe ° 
sin2 7 sin? = tg? ae sin? % 
neo ie co & oy 
bh) = L to 3 | 3 (n) ae rg 
ia aes a) = I = a > by) — — te 3 - , 
US Fie tg — sin? — tee sift 
6 6 aes 3 SVG ten 
sin az n paid 
ou sin—_ n cos —n 
b{") — 27 tg? 2 AC aioe en 
4 nM > i Ay d 13 g 9 ’ 
ge sin’ = tg 3 sin’ tg? 6 an = 
ay a P \ 
sin —n sin —n Cosme ny 
Hives tee a ee yes 3 
i 2 te = sin? 2 3 tg * sin? 2 Me = ‘8 2 
CaO: sue aih ae 


It I 
tee sin — 
3 KG 6 


‘Vel. W.-H. Isay, Z. angew. Math. Mech. 33 (1953) S. 405. 
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Alle weiteren bj”) sind Null oder vernachlassigbar klein (< 10-8). Selbst die Beriicksichtigung aller 
in (20) cdeetalirten bi", wird in vielen Fallen nicht nétig sein. Die Auflésung von (19) erfolgt nach 
der Theorie von Schineidlen auf dem tes tiber die Kernmatrix 1. 

Wir wahlen die beiden im Bereich 0... vollstindigen Orthonormalsysteme 


3) =//2 Sin 24. se, ¥,(0) = cosfo, Yo= oe (6,x =1, 2,...). 


Wird der Kern von (19) mit K(s, o) bezeichnet, so ergibt sich fiir die Kernmatrix 
Ko) = | K(s,0) ®,(s) Y9(0) do ds = 8,9 +5 0} (B,x=1,2,.--) 
6 0 
und _ speziell 


. n 1 
es Sie (Cad Re erro 


Die Auflésung von (19) ist damit auf die Lésung des Gleichungssystems 
M (m) — (m) 


ake Ypn = Sxn (Sel prceye) (21) 
: " () (m) 
zuriickgefiihrt. Die Konstanten g,,, und yz, stammen aus den Entwicklungen! : 
aE (m) (m) 1 (™) yes (ms) 
gals) =\/2 > 2, & sin x s, Vo(C) = x Yon \- 2% cosfo. (22) 


(m™) (m) 
Die g,,, sind also bekannt, die yz, sind zu bestimmen. 
In unserem Fall der ia eee (19) hat die Kernmatrix KY} die Form 


1 1 
geass ioe amen sale te 
1 1 
Ve bi"), 5: bm), 1+ > 8, 0, Osea. 
oe) 4 orien ere 
y2 bf), oy bt”), 0, 1, 0, 
0, 0, 0, 0, 1, 


Die Auflésung des Gleichungssystemes (21) ist somit sehr leicht méglich. Dabei bleibt eine Kon- 


(m) 
stante, Yon, frei; sie wird wie bei der Schaufelgitter- oder Tragfliigeltheorie durch die AbfluBbe- 


dingung an der Hinterkante des Propellerfliigels festgelegt. 
Das Ergebnis der Auflésung von oe ree tiie 7 == Os Nes oo Hh 


(m) 1 /(™ (m) 1 (m) il (m) 
: (Ein — 5 HP Ben ——y HP Bs ns HL Yon) 


Y Aiaiomeet A ig 82n 7 
(m) (m) 1 (m) 1 @& 1 eae 
Van Tein 2 A. Dips v5 onl Hh PeoEan ii 2) (23) 
PIE EELS oa ie Laem por) 
Yan = 83n— 97, . Se me 5 Yonl? 1) — aq Ot 31]? 
(m) (m) 
Vpn = Spn fiir B > 4 
Dabei ist 
22 
cos_n cos n 
A, =1—tg | 5 (— 1)" + ——— (24) 
’ ee Ne sin? = sin? 


1 Die Vertrautheit mit der von W. Schmeidler in seinem Buch: Integralgleichungen Bd. I, Leipzig 1950, ent- 
wickelten Theorie muB hier vorausgesetzt werden, speziell S. 50 des Buches. 
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In (23) und (24) sind bereits alle Glieder weggelassen, die wegen (18) so klein sind, daB sie keine Be- 
deutung haben. 
Die Lésung von (19) lautet also 


(m) a) ie M_ (m) ; eee (25) 
A = — n COSPO N= VU,4,+6- . 
v (a) Vx Yo n jas Xp 


(m) 
In den meisten Fallen wird M nicht groBer als 3 oder 4 sein. Zur Berechnung der yz, eae 


(m) (m) 
man die Kenntnis der g,,, aus der Entwicklung (22). Dafiir miissen die Funktionen te = a aC ) 


aus (16) auf die Variable s umgerechnet werden. Den Weg dazu deuten wir kurz an: Aus 


tg 4 =— tg > 008 s 


folgt mit der Abkiirzung tg = 6 2D. 


; 1—tecoss A 
Ch SS re COStS — es COSaS 
l-+diecoss 


Auf die Wiedergabe der an sich elementaren Umrechnung wird hier verzichtet. Aus der Kleinheit 
von ¢ ergeben sich wieder erhebliche Vereinfachungen. Im Endergebnis erhalt man folgende fiir die 
oben erwahnte Umrechnung wichtige Formeln: 


2 ~2 
ae BPN ots e242, 
eivsins =(1— 7) sins —ie sin 2s—->sin3s, 


e”'? sin s = (1 — 2 €”) sin s — 1 (2 e— 3 é*) sin2s—2e?sin3s+— ietsin4s, 

e?'” sin s = (1 — 4,5 €?) sin s— i (3e— 9,5 €°) sin2 s— 4,5e? sin 3 s + 4,757 3 sin 4s, 

et? sin s = (1— 8 2) sins — i (4¢ — 22 ¢°) sin 2 s — (8 e? — 36 €*) sin3 s (26) 
+lliesin4ds-+12¢sin5s, 

e'? sin s = (l— 12 ¢?) sin s—i (5 ¢— 41 ¢°) sin 2 s — (12,5 e? — 84,4 e*) sin 3 s 
+ i (21,25 63 — 125,25 ¢°) sin 4s + 28,1 e4 sin 5 s— 31,312 sin6s, 

e° sin s = (l— 18 ec? + 114 ¢*) sin s—i (6 ¢ — 708°) sin 2 s —(18 e? — 171 4) sin 3 s 
-+ i (36,5 6? — 295,5 ¢*) sin 4s + 57 et sin 5 s—73,9ie> sin 6s . 


(m) 
Die Konstanten yo, in (25) werden, wie schon gesagt, aus der AbfluBbedingung bestimmt. Diese 
lautet 


v(i—a,t)=0, dh. also a,(—o)=90, bzw. y,(0)=0. 
Unter Beriicksichtigung von (23) ergibt sich dann 


tt (m) 1 a) M (m) 
l (m) i (eu) S55; bP) fo + EBn 
= Mn ia = 5 
2 1 1 
aap ( xis bsp) bp — byp — bgp 


Nun ist 


a (™) 


(mF) ~ fom) = 
An =| a(y) dy =2 tg | ¥,(0) do =2 [ox te > Fon, 
a Oo 


und somit erhalten wir die wichtige Formel 


1 1 (m) 1 (m) — M(m) 
(m) _ giz ate a (1-27, en + Sem 
A, =—2y2n ; e=*—te5 (n=0,1,...7). (27) 
t= ve (1 Dees 5 oP bo — bY — hw 
Aus (25) erhalt man dann 
(m) 1 
an(y) = = — yale 
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durch Riicktransformation auf die Variable y. Wegen (18) kann man dabei ohne Bedenken den 
Tangens durch den Winkel ersetzen, so daB gilt 


2 3 
cosg~—, cos2a~ 241, cos3q~ 32 — 4h usw. (28) 


. Ci . . 3) (m) 
Damit sind die Integralgleichungen (15) gelést, und y(y,t) sowie /'(t) sind berechnet. 
Bei der effektiven Durchfithrung des Iterationsverfahrens zur Auflésung der urspriinglichen 
Integralgleichung (6) ergibt sich noch folgende Erleichterung: Bis zum vorletzten Iterationsschritt 


(m 
einschlieBlich braucht man nur jeweils aus Formel (27) 4, berechnen, waihrend auf die Bestimmung 


(m) 
der a,(y) verzichtet werden kann; denn fiir die Durchfiihrung des Iterationsverfahrens selbst braucht 
(m) (m) 
man ja nur /‘(t), nicht aber y(y,t). Erst beim letzten Schritt muB alles berechnet werden. 
Im konkreten Zahlenbeispiel berechnet man zunichst die bei jedem Iterationsschritt gleich- 
bleibenden GréBen bi"). 
Das Iterationsverfahren kann dann in folgender Weise schematisiert werden: Es sei bekannt die 


linke Seite der Integralgleichung (6) in der Form (14) durch Einsetzen der (m— 1)-ten Naherung 
(m—1) (m) 

I’ (t) in das Doppelintegral. Dann sind auch die Funktionen f,(p) nach (16) bekannt. Diese werden 
darauf mit Hilfe der Formeln (26) auf die Variable s umgerechnet, d. h. es werden die Funktionen 


(m) (m) 
gn(s) gebildet; danach kann man unmittelbar die Konstanten g,,, ablesen. Diese werden in (27) ein- 


(m) (2) 

gesetzt; damit ist /(t) bestimmt. Die Ableitung /'(t) wird dann in das Doppelintegral auf der linken 
Seite von (6) eingesetzt und dieses mit Hilfe von (9) bzw. (10), (11) ausgewertet. SchlieBlich wird 
durch eine harmonische Analyse mit der Schrittweite 7/6 eine fiir 0 < ®< 2 z giiltige Darstellung 
(13) des Doppelintegrals hergestellt. Und darauf kann der nachste Iterationsschritt beginnen. 

Die Konvergenz des Iterationsverfahrens hangt von der GréBe tg ¢ =e ab. Dieses kann man 

m) (m—1) 

leicht erkennen, wenn man etwa die Differenz ae A, nach der Formel (27) und unter Benutzung 
von (16) und (26) bildet. Man erhalt dann etwa einen Ausdruck der Form 
(m)  (m—1) ((m) (m—1) , (m) = (m—1) ff (m) — (m—1) 
‘vas Ay, SS ate G,,(€) ee — Cy =o oR n—1(€) ome an | oe E me oan i(é) Ve ere Kens i . 


Dabei sind G,,, Gz, Gi* fiir jeden Iterationsschritt gleich. Da rj/2R klein (etwa 1/8 bis 1/4) ist, 
spielen der zweite und dritte Summand fiir die Konvergenz nur eine geringe Rolle, wie sich auch 
bei numerischen Berechnungen gezeigt hat. 

Unschwer kann man nun die pe bzw. iE nach (13) und mit Hilfe der Relation (9) wieder mit 

m—1 
aoe i des vorigen Iterationsschrittes in Verbindung bringen, und auf diese Weise eine Ab- 
schatzungsformel fiir den Konvergenzbereich ¢ < ¢) gewinnen. Auf die Wiedergabe dieser Unter- 
suchung wird verzichtet, denn bei der Abschatzung der recht umfangreichen und umstandlichen 
Ausdriicke ist es unvermeidlich, daB der Konvergenzbereich ¢, kleiner herauskommt als er in Wahr- 
heit ist. An unserem ausfiihrlich behandelten Beispiel werden wir sehen, daB fiir ¢ = 0,0875 das 
Iterationsverfahren ohne weiteres konvergiert; bei Proberechnungen zeigte sich, daf} auch fiir 
¢ = 0,1 noch eine wenn auch langsamere Konvergenz vorliegt. Somit kann im Rahmen der fiir die 
ganze Theorie geltenden Einschrankung (18) die Konvergenz des Iterationsverfahrens als gesichert 
angesehen werden. 


6. Der Beweis fiir die Integralformel (9). Um den Gedankengang der Auflosung von (6) nicht zu 
unterbrechen, hatten wir die Integralformel (9) ohne Beweis mitgeteilt. Dieser soll nun nachgeholt 
werden. Wir fiihren zunichst die Integration iiber # aus und setzen zur Abkiirzung 

20 
een ein 9d) (nee 1). 
27. ae 
f 
Durch die Transformation e'” = z, d) = dz/iz geht J, uber in 


1 Ze? __%i yzsin OW —e'® —ih(t)[z2e—'? +24 2cosP—2z+e'*| 
I,.= — fa! = d 


Bs 
: : 1 
27 GfK) (yx en, '*) (z et? f- x) (: oS 


Jn 
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Der Integrationsweg ist der Einheitskreis der z-Ebene. In dieser Form kann J, leicht nach der 
Residuenmethode ausgewertet werden. Dabei sind zwei Falle zu unterscheiden : het 

1) Es liegt z, = 1/(e—'® — y) im Einheitskreis; dann liegt z, = e'® — x nicht im Einheitskreis, 
und es ergibt sich 

ae 1 
Urea UNL MA pears great 

2) Es liegt z, = e'® — y im Einheitskreis; dann liegt z,; = 1/(e~'® — x) nicht im Einheits- 

kreis, und es ergibt sich 
Jn =— (i— h(t)) e® (et? — yn) 
Der Grenzfall |z,| = |z.| = 1 bedeutet 
1+ 77—2ycoss®=1, 

dh. y = 0 oder y = 2 cos @®. Dann ergibt sich J, als arithmetisches Mittel der beiden unter 1) 
und 2) gefundenen Ergebnisse, was man auch durch direkte Integration fiir y = 0 und y = 2 cos ® 
bestatigt. 

Fiir die nun noch folgende Integration iiber y sind zwei Falle zu unterscheiden: 

1) Es ist cos@ negativ, d.h. 7/2<®<3 7/2, dann ist fiir 0<y< oo immer |z|>>1, und es wird 

oe Ber ey ee dy a die ey (ean) 
wv ~=—(i+ (t)) e (ay = (et ( ))e = aya 
0 ) 
2. Es ist cos ® positiv, d.h. —2/2< ®<a/2, dann ist fir 0<y<2cos@ lza|<< 1, und 


es wird somit 


2 2cos@ a6) 
| Jn dy —— 1 — hier | (ey a alts) ome | ( ~at yt 
0 0 2cosP a 

= oars (i— h(t)) e—i(n—l)® — (i— h(t)) e+ Ho 4 os (i ae h(t) ) e—i(n+1)o (n = 1) ‘ 


Fir ® = 7/2 und ® = 37/2 stimmen beide Werte iiberein. Damit ist Formel (9) bewiesen. 

An dieser Stelle mégen noch zwei weitere Integralformeln angegeben werden, die mit (9) in 
einem gewissen Zusammenhang stehen, und deren Beweis ganz analog verlauft, so daB wir ihn hier 
nicht wiedergeben brauchen. Es ist fiir n > 1 

co 27 


o | / e'”® (sin & — sin ©) dd dy 
2 l 
see Pay X= 008(P — 8) =— 7 cos O-> 7-cos@ 


0 0 
ae eth, fire ea ad 4 
n ep 
Pelgite eae . hie me & ee: 
re ae) e€ fiir Sp es 
und 
o 22 
Io i"? (cos D — cos # — x) ddd 17 
al - (o98 COSA, =—— er? fir 02 O227. CO) 
J a 1 + x? — cos (BP — 3) — y/cos B + x cos B o 


Aus (29) und (30) folgt ohne weiteres 


o 2n 


1 a cos n # (sin & — sin OD) dd dy 
2a Lae = 
J 1+ x* — cos (P — 6) — x cos ® + x cos 


a ges 2 
Linn BF (— 1)" Sinn ® fiir —FSeO<F 
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oO 27 


ae sin n @ (sin & — sin ®) dd dy 
2 7 1 
Bei. Dae fos D 8) y con® Fy cos 0 


0 0 
1 _ Tt 3% 
——cosn® fir =< G6< — 
n eect ee ee 
aan 2 Tt I (32) 
pcos n BP —— (— 1)" cos n® fiir Fat Oe 
o 22 
1 cos n & (cos D — cos & — yx) ddd ih 
at . Ee a cosn® fir 0X G@<2z, (33) 
1+— 7? — cos (® — 3) — x cos ® + x cos - 
0 0 2 ; 
oan 
inf sin n @ (cos D — cos — y) ddd il 
eal i 1) a9 diy -: sinn® fir 0< @<2n. (34) 
5 => 
a Me t= 00s (P= 0) y cos Dy cos 0 


Diese Formeln werden benétigt, wenn man die Geschwindigkeiten berechnen will, die von den ab- 


flieBenden freien Wirbeln am Ort der Propellerfliigel induziert werden [Vegl. (4)]. 


7. Zahlenbeispiel zur Berechnung der Zirkulationsverteilung. Wir wollen nun das Iterationsver- 
fahren zur Auflésung der Integralgleichung (6) an einem Zahlenbeispiel durchfiihren. 

Dabei wollen wir annehmen, da® der Anstellwinkel Null der Propellerfliigel wie in Abb. 1 bei 
m/2 und 37/2 liege. Dann ist jad = 0. Die Fliigel des Propellers mégen sich tiber 20° erstrecken; 
also ist a = 10° und tga/2 = 0,0875. Der Fortschrittsgrad unseres Propellers sei u,)/Rw@—= 1/16. 
SchlieB'ich sei r5/R = 0,25. Dann ist der maximale Anstellwinkel der Propellerfliigel gegen die 
Kreisbahn gleich 14°. Fiir die Rechnung ist es dann noch zweckmaBig, sich iiberall y durch y/uy 
ersetzt zu denken und dafiir wieder y zu schreiben; d. h. mit anderen Worten y in Einheiten von 
Uy zu berechnen. Man hat so reine Zahlen und muB sich den Faktor uy, jedesmal dazu denken. 

Mit den angenommenen Werten berechnen wir nach (20) und (24) 


ee aoa, b= — b= 0,180, bb bi) 0, 
oS!) — 6) — 0,007, o2 = of) =— 0,017, {9 = bf) = 0,046, bi’ = — 0,026, (35) 
ee  -0.0027,. b= 0,002 1, BY) = bb 0, 
by = — 4 by), (5) = 2 by), 
AA e004 An — A, 1,017, Ag = 1, — 0,990,450 240 (36) 
Aus (26) erhalten wir 
e'? sin s = 0,996 sin s — 0,0875 i sin 2 s— 0,004 sin 3s , 
e2'” sin s = 0,985 sin s — 0,173 7 sin 2 s— 0,015 sin 3 s 4+ 0,001 isin4s , 
ei” sin s = 0,966 sin s — 0,256 i sin 2 s— 0,034 sin 3s + 0,003 7 sin 4s, (37) 


et? sin s = 0,939 sin s — 0,335 sin 2 s— 0,059 sin 3 s + 0,007isin4s + 0,001 sin5s, 
&'” sin s = 0,908 sin s— 0,410 i sin 2 s— 0,091 sin 3 s + 0,014i sin4s + 0,002 sn 5s, 
e'? sin s = 0,869 sin s — 0,478 7 sin 2 s— 0,128 sin 3 s + 0,023 7 sin 4s + 0,003 sin5s. 


Die Durchfiihrung des Iterationsverfahrens erfolgte genau in der angegebenen Weise (vgl. 


S. 389). Das Ergebnis der ersten 6 Iterationsschritte fiir die Zirkulationsverteilung /’(t) (m = 1, ... 6) 
sowie fiir die Funktionen 
6 (m+1)_ 6 (m+) QQ). (1) 
Signe See C= 0) 
p=—O p=—O 

aus (13) sind in Abb. 4, 5 und 6 graphisch dargestellt. Die neben den Kurven befindliche 
Zahl gibt die Nummer der Naherung an. Man sieht deutlich, wie das Verfahren konvergiert. 
Um das Verfahren etwas abzukiirzen! haben wir als Ausgangspunkt fiir die 7. Iteration 


Ue 


1 In der Praxis kann man das Iterationsverfahren wesentlich abkurzen, indem man von vorne herein solche 
Mittelwerte bildet und das Verfahren gleich mit einer geschickten Schatzung beginnt. 
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Abb. 4. Die ersten sechs Naherungen der Funktion Ze,eil! ® heim Iterationsverfahren. 
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Abb. 5. Die ersten sechs Naiherungen der Funktion Zk eu P beim Iterationsverfahren. 
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eh 


Abb. 7. Die letzten drei Naherungen der Funktion Be elie beim Iterationsverfahren. 
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(6) ; 6 (7) ae 
nicht (wie eigentlich vorgesehen) die aus der 6. Naherung /"(t) entstandenen eee a Cu CA 


ke"? gewahlt, sondern einen Mittelwert von 


ia 6 (6) 
6 (7) 6 (6) 6 ae 
5, epee und De, eH ta Dawe > k,e#® und ale et HO, 


p=—6 p=—6 


We 9 liegt 2wischen 8 und 6,7 


Abb. 8. Die letzten drei Naherungen der Funktion Zk ,eu? beim Iterationsverfahren. 


Abb. 9. Die letzten drei Naherungen der Zirkulationsverteilung I beim Iterationsverfahren. 


(Vgl. Kurve 6,7 in Abb. 7 u. 8). Ebenso wurde als Ausgangspunkt fiir den 9. Iterationsschritt 


ein Mittelwert von 


6 (9) 6 (8) 6 (9) 6 (8) 
iD) 6 eee and DP lee TANG Dy ke! eennd Di hnere 
w=—O6 u=—6 u=—O6 u=—6 


genommen. Mit der 9. Iteration kann das Verfahren abgebrochen werden, wie sich auch aus Abb. 7, 
8 und 9 ergibt. Die sich bei dem 9. Iterationsschritt ergebenden Funktionen stellen wir nachfolgend 
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ausfithrlich zusammen. Dabei lassen wir die obere Kennzeichnung (9) weg, da die 9. Naherung als 
die endgiiltige Lésung anzusehen ist. 
Es wird nach (13) 


co 27 
1 eet vaee 6 | 6 ; 
2 Ww | | 1 2) aa oe di dy af DS Pt e He -|- h(t) pa ue eb uP = 
0 0 i w=— 


= — 3,122 cos O— 1,412 cos 2 d— 0,232 cos 3 © + 0,134 cos 4D +. 0,050 cos 5 OD 
— 0,026 cos 6D + 0,072 sin 2 O + 0,016 sin 3 6 — 0,042 sin 4 ® — 0,046 sin 5 @ 


1 
JE 7 coswt (— 0,383 — 0,042 cos 2 ® + 0,016 cos 3 © + 0,046 cos 4@ + 0,020 cos 5 @ 
— 0,014 cos 6 ® + 0,484 sin @ + 1,412 sin 2 ®+ 0,870 sin 3 d — 0,134 sin 4 © — 0,448 sin 5 @). 


Ferner ist nach (16) (alles in Einheiten von Up) 
Joly) = — 0,155 i ef” + 0,155 i e—i#, 
3 ; 
Sil) = 3,952 + = Ay— 2,561 ef — (0,003 + 0,088 i) ei”, 


f(p) = — 0,155 i e'? — (0,706 + 0,036 i) e2'” +. (0,001 — 0,054 i) e*” , 
f3(~) =— (0,116 + 0,008 7) e'” — (0,003 +. 0,088 i) e2'” + (0,003 -- 0,008 i) et”, 
Sa(p) = (0,067 + 0,021 i) et’ +. (0,001 — 0,054 i) e3' + (0,001 + 0,028 i) eF#2, 


fs() = (0,025 + 0,023 i) e54 + (0,003 + 0,008 i) et” — 0,001 e” + > 4,, 
fe() =— 0,013 e%'? + (0,001 + 0,028 i) edé#, 
fp) ~ 0. 


Die Rechenreihenfolge ist so, da f,(~) vor f;(~) behandelt wird, da ja A, bekannt sein muB, bevor 


fs(p) weiter behandelt werden kann. Die Umrechnung auf die Variable s mit Hilfe der Formeln (37) 
ergibt 


&(s) =— 0,027 sin 2 s, 

&1(s) = (1,400 — 0,087 7) sin s + (— 0,015 + 0,224 i) sin 2 s + (0,010 + 0,001 i) sin 3s, 
&2(s) = — (0,694 + 0,242 7) sin s + (— 0,034 + 0,122 i) sin 2 s + (0,011 + 0,003 1%) sin3s, 
&3(S) = — (0,112 + 0,087 z) sin s + (— 0,015 + 0,029 7) sin 2s + (0,004 + 0,001 2) sin3s, 
&1(s) = (0,065 — 0,007 i) sin s +- (0,005 — 0,023 i) sin 2 s — (0,004 + 0,002 7) sin-3 s, 
g5(s) = (0,025 + 0,027 7) sin s + (0,012 — 0,011 7) sin 2 s — (0,002 + 0,002 i) sin 3 s, 


&6(5) = (— 0,010 + 0,025 7) sins + (0,011 + 0,006 i) sin 2 s + (0,002 — 0,003 i) sin 3s. 
Unter Beriicksichtigung der AbfluBbedingung ergibt sich entsprechend Formel (25) 
¥o(a) = 0,029 (1— cos 20) , 
¥,(a) = — 1,290 + 0,336 i + (1,288 — 0,551 i) cos o + (— 0,007 + 0,216 1) cos 20 
+ (0,009 — 0,001 7) cos 30, 
¥(o) = 0,691 — 0,013 i — (0,680 + 0,116 t) cos o + (— 0,022 4 0,125 7) cos 26 
+ (0,011 + 0,004 i) cos 30, 
y3(0) = 0,118 + 0,054 i — (0,110 + 0,085 7) cos o + (— 0,012 + 0,030 1) cos 26 


+ (0,004 + 0,001 i) cos 3c, 

y4(c) =— 0,067 + 0,020 i + (0,067 + 0,005 i) cos o + (0,004 — 0,023 i) cos 26 
— (0,004 + 0,002 i) cos 3c. 

y;(0) =— 0,027 — 0,024 i + (0,017 + 0,037 i) cos o + (0,012 — 0,011 i) cos 2.6 
— (0,002 + 0,002 i) cos 3c, 

y¥e(0) =— 0,002 — 0,030 i -+ (— 0,011 + 0,026 i) cos o + (0,011 + 0,007 i) cos 20 


+ (0,002 — 0,003 i) cos3c. 
Mit Hilfe der Formeln (27), (28) erhalten wir entsprechend unserem Ansatz (7) die endgiiltige 
Zirkulationsverteilung ; es ist 
I(t) = 0,016 — 1,418 cos m t— 0,370 sin wt + 0,760 cos 2 wt + 0,014 sin 2 wt 
0,130 cos 3. wt— 0,060 sin 3 wt— 0,074 cos 4.@ t— 0,022 sin 4 t— 0,030 cos Sat ¢ (38) 
+ 0,026 sin 5 m t— 0,002 cos 6wt + 0,034 sin6mt, (in Einheiten von uy) . 
al 
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A Sees ; 
y(y, t) = 0,058 1 v4 Vi a ( 2,566 -+ 0,044 — 0,072 = cos t 


2 Oo Oo 


2 2 
— (0,240 + 0,856” + 0,008 5 oes (1,426 — 0,088 3) foster 


' 2 
3 (0276 een seoe tae 3) ne ree (0,260 — 0,016” — 0,032 i) cost 


~ 


4 
ide cos5at 


) 2 2 
(0,048 + 0,128” — 0,008 ») anise — (0,142 + 0,016”. — 0,032 3 cos 4at| (39) 


0,086 — 0,108 © + 0,016 5] ld oe 0,078 + 0,032" 4 0,016 


a 


2 


2 

) 
; 2 ‘ 2 
5 (0,026 + 0,060” — 0,016 s Sey one (0,026 — 0,060". + 0,016 a) coset 
i (0,074 0,004” — 0,024 7) SG: | (—aSypSa; a 


= 40°)e 


Die Formeln (38) und (39) gelten fiir jeden beliebigen Zeitpunkt t. Ist z. B. fiir t = t) der Mittelpunkt 

des Propellerfliigels 1 an der Stelle wt, der Propellerkreisbahn angekommen, so ist eine momentane 
—l1 ; ; 3 

Zirkulation y(y, t)). Zugleich ist dann y (v, ty + . =) (p = 2,...6) die momentane Zirkulation 


@ 


der anderen fiinf Propellerfligel. 


8. Bestimmung der resultierenden Kraft des Voith-Schneider-Propellers. Um die resultierende 
Kraft des ganzen Propellers zu bestimmen, miissen wir zunichst die resultierenden Krafte der ein- 
zelnen Propellerfliigel berechnen. Dazu wenden wir den Kutta-Joukowskischen Satz sinngemaS 
an: Die resultierende Kraft K steht senkrecht auf der resultierenden Geschwindigkeit w (relativ 
zu dem als ruhend gedachten Profil) und hat die GréBe K = 9 w I’. Dabei ist J’ die Gesamtzirku- 
lation der an das Profil gebundenen ,,tragenden“‘ Wirbel. w ist die Gesamtgeschwindigkeit, die am 
Ort des Profils relativ zu diesem herrscht, abziiglich der von den gebundenen Wirbeln des Profils 
selbst induzierten Geschwindigkeit. Die von den abflieBenden, freien Wirbeln induzierte Ge- 
schwindigkeit ist natirlich in w enthalten. (Kine Aufteilung der resultierenden Kraft in Auftrieb 
und Widerstand wie in der instationaren Tragfliigeltheorie, wo der Widerstand durch die abgehen- 
den, freien Wirbel hervorgerufen wird, hat beim Propeller keinen Wert.) 


Die genannte Geschwindigkeit w = jv? + V? setzt sich in unserem Fall in einem Punkt @t der 
Kreisbahn aus folgenden Anteilen zusammen. 

1. Die konstante Anstrémung Up. 

2. Die Umfangsgeschwindigkeit mit den Komponenten 

(— u,) =o Rsinwt, (—v,) =—w Reoswt 

relativ zu dem als ruhend gedachten Propellerfliigel. 

3. Die von den freien Wirbeln am Ort der Propellerfliigel induzierten (absoluten) Geschwindig- 
keitskomponenten uy, vy. 

4, Die von den gebundenen Wirbeln der anderen fiinf Propellerfligel am Ort des betrachteten 
Fliigels induzierte (absolute) Geschwindigkeit mit den Komponenten 3/'ur, >)’ vp. 

Somit ist 


U=u+o Rsinwt-+ us(wt) +d” u(t); V=—o Reoswt-+ vs(wt) + Dd” v/-(w 2). 
Da sich w tiber die Lange des Propellerfliigels der GréBe und vor allem der Richtung nach andert, 
soll als die fiir die Berechnung der Kraft nach dem Kutta-Joukowskischen Satz maBgebende Ge- 


schwindigkeit w diejenige genommen werden, die jeweils im Druckmittelpunkt des Fligelprofiles 
herrscht. Dabei ist der Druckmittelpunkt y, definiert durch 


Vol!) = iG | p vy,t) dp . (40) 


——y 
Die Auswertung des Integrales in (40) wird ganz einfach, wenn man auf die Variable ¢ zuriickgeht. 
Die Bestimmung des unter 3) genannten Geschwindigkeitsanteiles erfolgt nach Gleichung (4) 
unter Verwendung der Integralformeln (31), (32), (33), (34). Sie ist ohne weiteres moglich. 
Fir die Bestimmung des unter 4) genannten Geschwindigkeitsanteiles kénnen wir uns jeweils 
die fiinf Propellerfliigel durch Punktwirbel mit den entsprechenden Gesamtzirkulationen ersetzt 
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Jenken. Die Geschwindigkeit, die ein im Punkt Re!” der Propellerkreisbahn angebrachter Wirbel 


R.J‘(t*) in einem anderen Punkt R ce‘! der Kreisbahn induziert, ergibt sich nach der elementaren 
Theorie zu 


T(t*) cos + @ + t*) T(*) sin > (t + t*) 
ur(wt) = = , vp(wt) =- = : 
sin > (¢ —t*) sin iy (t —t*) 


Unter Benutzung dieser Formeln kénnen die unter 4) genannten Geschwindigkeiten leicht be- 
rechnet werden. 

Hat man nun auf diese Weise die resultierenden Krafte der einzelnen Propellerfligel bestimmt, 
so kann man aus diesen ohne weiteres die resultierende Gesamtkraft des Voith-Schneider-Propellers 
zusammensetzen. 

Wir wollen nun die eben allgemein dargelegte Bestimmung der Krafte an unserem Beispiel von 
Abschnitt 7 durchfithren. Die Zahlenwerte der Geschwindigkeiten sind auf Einheiten von u, be- 
zogen. In Tabelle 1 sind die einzelnen oben erwahnten Anteile der Gesamtgeschwindigkeit w 
(relativ zum Propellerfliigel) zusammengestellt. 


Tabelle 1. 

ot a’up Lon sf Pf —w, =v, Uy 
0 | -+0,120 —0,056 | +2,319 +0,189 0,000 —16,000 --+1,000 
7/6 | +0,035 O08 | SEILRRt +0,179 -+-+ 8,000 —13,856 +1,000 
r/3 —0,011 —0,270 | +0,891 —0,206 +13,856 — 8,000 +1,000 
me /2 +0,096 | —0,443 | —0,558 —0,754 +16,000 0,000 +1,000 
22/3 -+0,272 —0,482 | +0,068 —1,037 +13,856 + 8,000 +1,000 
5 7/6 +0,321 —0,305 +0,698 —0,757 + 8,000 +13,856 +1,000 

It +0,132 | —0,154 +1,001 —0,215 0,000 --16,000 +1,000 
72/6 —0,136 | —0,242 +0,923 +0,431 — 8,000 +13,856 +1,000 
Az/3 —0,161 | —0,408 +0,441 +0,887 —13,856 + 8,000 +1,000 
3 7/2 —0,051 | —0,441 —0,274 +0,824 —16,000 0,000 +1,000 
5 2/3 +0,119 | —0,359 | +1,326 +0,277 —13,856 — 8,000 +1,000 
ll 2/6 +0149 | —0,145 +2,175 +0,179 — 8,000 —13,856 +1,000 


Aus Tabelle 1 entnehmen wir leicht die Geschwindigkeiten U, V und w =) U? + V?. Die Richtung 
von w gibt f = arc tg V/U an (gegen die positive x-Achse gemessen). Das Ergebnis ist in Tabelle 2 
zusammengestellt. 


Tabelle 2. 

wt U V w | B | aot | U V w B 

0 + 3,439 | —15,867 16,23 —17,717° IU + 2,133 | +15,631 15,78 + 82,23° 
/6 +10,885 | —13,759 | 17,54 —91,65° 72/6 | — 6,213 | +14,045 13,37 +113,86° 
r/3 +15,736 | — 8,476 17,87 I BWAe Ag/3 | —12,576 | + 8,479 LSS +146,01° 
m/2 —- 16,538 || — 1,197 16,58 eB 3Y* 32/2 | —15,325 | + 0,383 15,33 +178,57° 
2/3 | +15,196 | + 6,481 16,52 +23,10° 52/3 |—11,411 | — 8,082 13,98 —144,68° 
5 2/6 | +10,019 | +12,694 | 16,17 | +51,72° |1l 2/6 |— 4,676 | —13,822 14,59 —108,69° 


Die Berechnung des Druckmittelpunktes nach Formel (40) liefert das in Tabelle 3 wiedergegebene 
Ergebnis 


Tabelle 3. 


a  ———  ———_—_—__OO 
wt | Po| | ot | Wol & | ot | Po|& | ot Po[ O ot Wo| | ot | Po| O 
Om 0,51 60° | 0,62 | 120° | 0,24 | 180° | 0,50 | 240° | OS2e I s002s 20,28 
iis" 0,60 RN OSOLe | eLaceueOs45e 9521 0s500 |e2oogm 0, cleriolomsN0;25 

30° 0,64 90° | 0,59 | 150° | 0,47 | 210° | 0,51 | 270° | 0,34 | 330° | 0,31 

45° | 0,63 | 105° | 0,60 | 165° | 0,49 225° | 0,52 | 285° | 0,32 | 345° | 0,41 


Nun haben wir alle GréBen beisammen, die fiir die Berechnung der resultierenden Krafte an den 
Propellerfliigeln benétigt werden. In Tabelle 4 sind sie iibersichtlich zusammengestellt. Dabei ist 
die resultierende Geschwindigkeit w im Druckmittelpunkt nach GréBe und Richtung mit Hilfe von 
Tabelle 2 interpoliert. Die angegebenen Richtungswinkel sind wie ublich gegen die positive x-Achse 


2708 
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gerechnet. Ferner sind, wie schon gesagt, /’ und w in Einheiten von u, und entsprechend K in Ein- 
heiten von 0 u; R angegeben. 


Tabelle 4. 
i nn UE EE EET EEEEEEEEEEEEEEEEEE 

Lage de | rage des Druck. | sult, Ceschorin~| Richtune von w | profiliskulation | cout ‘rafe | result, Kraft | des Propeller 
mittelpunktes mittelpunktes Druckmittelpkt. mittelpunkt Tr K K fliigelprofiles 
( SRS 16,45 — 73,4° —0,618 10,17 —163,4° — 1716,0° 
G5} 21,0° 17,21 | — 59,5° —0,740 IPA) —149,5° — 61,4° 
30° 36,4° 17,68 | — 46,7° —1,006 Wit) —136,7° — 47,8° 
45° les 17,92 — 35,1° —1,324 23,73 —125,1° — 35,0° 
60° 66,2° 17,60 — 23,3° —1,495 26,31 —113,3° — 22,9° 
on BIL 16,75 — 11,4° —1,416 23,12 —101,4° — 11.3° 
90° 95,9° 16,57 -+—-- 1,2° —1,100 18,23 — 88,8° 0,0° 
105° EEO 16,53 + 14,9° —0,559 9,24 — 75,1° + 11,3° 
120° 122,4° 16,49 + 25,4° +0,151 2,49 +115,4° + 22,9° 
Son 13979) 16.33 + 41,7° --0,862 14,08 +131,7° + 35,0° 
150° 154,7° 16,11 ="D0,0— +1,412 22519 +146,5° + 47,8° 
165° 169,9° 15,90 se UO +1,787 28,41 +161,9° + 61,4° 
180° 185,0° Sail se doy" +2,018 31,70 -+-177,5° + 76,0° 
195° 200,0° 15,49 SENOS 3 +2,059 31,89 —166,7° + 91,4° 
210° INS I 15,33 +-119,3° +1,862 28,54 —150,7° +107,8° 
229" 230,27 15,18 4135,4° +1,464 22,22 —134,6° +125,0° 
240° 245,2° 15,20 = 151,67 -+-0,899 13,66 —118,4° +142,9° 
ASS) 262,15 15,32 sl Os Ole +0,177 Zell —100,0° +161,3° 
RAO? 273,4° 15,18 —177,3° —0,932 8,08 + 92,7° 0,0° 
285° 288,2° 14,65 —159,1° —0,918 13,45 +110,9° —161,3° 
300° 302,8° 14,04 —141,4° —0,871 12,23 +128,6° —142,9° 
Se ST 14,23 —1123,7- —0,642 9,13 +146,3° —125,0° 
330° Bosh” 14,76 05,0 —0,528 TW) +164,5° —107,8° 
345° 349,1° 15,60 — 89,0° —0,558 8,70 —179,0° — 91,4° 


Um nun die resultierende Gesamtkraft des Voith-Schneider-Propellers fiir verschiedene Propeller- 
stellungen aus den Kraften der einzelnen Fliigelprofile zusammensetzen zu kénnen, berechnet man 
am besten aus den in Tabelle 4 gegebenen Werten die Kraft K zerlegt in ihre x- und y-Komponente 
K, und Ky. 

Das Ergebnis gibt Tabelle 5 an. 


Tabelle 5. 


Mittelpunkt des i un e i n e 
Flit ecipratiles ete 2 : Ky VE ioee — : K;, Ky 
0° — 9,75 — 2,90 120° — 1,07 = 2,25 240° —6,50 —12,01 
GY? ——10,97 — 6,46 iLasss° — 9,36 +10,51 Ds)” —0,47 — 2,67 
30° —12,95 —12,20 150° —18,97 +12,56 DOs —0,38 + 8,07 
45° —13,64 —19,41 IOS” —27,00 + 8,83 285° —A4A,80 +12,56 
60° —10,41 —24.16 180° —31,67 + 1,38 300° —7,63 + 9,56 
oie — 4,69 —23,25 195° —31,03 — 7,33 Biles? —7,60 + 5,07 
90° + 0,38 —18,22 210° —24,89 —13,97 5308 —7,50 + 2,08 
105% + 2,38 — 8,93 2256 —15,61 —15,82 345° — 8,70 — 0,15 


Somit ergibt sich fiir die resultierende Gesamtkraft unseres Voith-Schneider-Propellers durch ein- 
fache Addition der Werte aus Tabelle 5 folzendes Endergebnis: 


in der Stellung 0 des Propellers: K(?) —— 67,03 pup R; K{?) =—25,88 ouzR; K(?) —71,84 ou2R; 


Richtung: Bx = — 158,9°; 

in der Stellung 7/12 des Propellers: K(?) =— 64,12 ou R; K) =— 24,13 9pujR; KP) = 68,51 9 uZR; 
Richtung: px = — 159,4°; 

in der Stellung 77/6 des Propellers: K(?)=— 64,31 ougR; K{) =—21,68 pu2R; K(?) 67,86 ou,R 
Richtung: bx = — 161, 4°; 

in der Stellung 77/4 des Propellers: K(?) — —67,370ugR; KY) = —22,920u2R; K®)=71,170u2R; 
Richtung: By == 161,2°, 

Die Stellung 7/3 ist wieder identisch mit der Stellung 0 usw. 

Als Mittelwert ergibt sich K”) — — 65,71 9 uy R; K?) = — 23,65 9 u2R; K() = 69,84 ou.R; 


Richtung: Bx = — 160, 2°, 
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Man erkennt, da®B die resultierende Kraft nicht genau entgegengesetzt zur Anstrémrichtung ge- 
richtet ist, sondern mit dieser einen Winkel von etwa 19,8° bildet. Bei der hier betrachteten Lage des 
Anstellwinkels Null der Propellerfliigel bei z/2 und 3 7/2 liefert der Propeller also keine genaue Vor- 
wartsfahrt, sondern eine Fahrtrichtung, die um 19,8° im Sinn der Rotation des Propellers gedreht ist. 


9. Lage des Anstellwinkels Null der Propeilerfliigel bei genauer Vorwirtsfahrt. Wir wollen nun 
untersuchen, an welcher Stelle der Anstellwinkel Null der Propellerfliigel gegeniiber der Kreisbahn 
legen muB, um genaue Vorwartsfahrt (d. h. resultierende Propellerkraft entgegengesetzt zur An- 
strémung) zu erreichen. Um diese Frage streng zu beantworten, miiBte man die ganze Iterations- 
rechnung zur Bestimmung der Zirkulationsverteilung mit einem offen gelassenen 4 wiederholen 
und 4 erst zum SchluB aus der obigen Bedingung bestimmen. (Wir hatten in unserem Beispiel 
A= 0 gesetzt; vgl. 5.391). Dann werden J’ und auch w und f =: arc tg V/U Funktionen von J. Die 
Rechenarbeit (man denke an arc tg!) ware sehr umfangreich, und deshalb soll dieser Weg hier nicht 
gegangen werden. 

Kine angenadherte Lésung fiir die gestellte Aufgabe kann man sich jedoch verschaffen, indem 
man der Grundanstrémung uy in x-Richtung noch eine kleine Komponente v, in y-Richtung iiber- 
lagert. Unter der Voraussetzung, dafs dann in der Integralgleichung (6) lediglich auf der linken 
Seite noch die Glieder 
ro 


— 2 v9 sin (p + © t)—2 % > 


cos w t cos (p + © t) 
hinzutreten, kann man den durch v, hervorgerufenen Zusatzanteil /‘* der Zirkulation getrennt be- 
rechnen und fir die Hauptzirkulation /' das bisherige Ergebnis (38) beibehalten. 

Diese Voraussetzung trifft aber eben nicht streng zu, denn die freien Wirbel flieBen jetzt nicht 
mehr in Richtung der x-Achse ab, sondern um den Winkel arcig vj/u) gegen die x-Achse geneigt. 
Solange v)/u, klein ist, kann indessen die Beriicksichtigung dieses Umstandes bei einer Naherungs- 
rechnung unterbleiben, um so mehr als nur die in der Nahe des Propellers befindlichen abflieBenden 
freien Wirbel eine wesentliche Rolle spielen. 

Ferner wollen wir den EinfluB, den die zusatzliche Zirkulation /** auf GréBe und Richtung der 
resultierenden Geschwindigkeit w hat, vernachlassigen. Seine Beriicksichtigung wirde dieselben 
umstandlichen Rechnungen wie oben geschildert erforderlich machen. AuBerdem ist dieser EinfluB 
wegen des Uberwiegens der Umfangsgeschwindigkeit klein; und auch wegen des in der Integral- 
gleichung begangenen prinzipiellen Fehlers mit der AbfluBrichtung der freien Wirbel hat eine zu 
genaue Rechnung jetzt keinen Zweck. 

Zur Bestimmung der zusatzlichen Zirkulationsverteilung wird das bekannte Iterationsverfahren 
angewendet. Das Ergebnis ist 


I*(t) = vp/u,y (0,14 — 0,255 cos w t 4+ 1,045 sin w t + 0,255 cos 2m t + 0,03 sin 2 @ t 
+ 0,03 cos 3 wm t— 0,03 sin 3 w t— 0,01 cos 4m t— 0,01 sin 4.@ t) . 


Es sei erwahnt, daB zur Berechnung von /™* lediglich drei Iterationsschritte erforderlich waren, da 
gleich mit einer geschickten Schatzung begonnen wurde. 

J* ist in Einheiten von u, berechnet. In Tabelle 6 ist die zusatzliche Zirkulation J‘* zusammen 
mit den von ihr an den Propellerfligeln hervorgerufenen Zusatzkraften K* nach GroBe und Rich- 
tung angegeben. Bei der Berechnung von K* wurde (wie oben erlautert) der Wert w aus Tabelle 4 
benutzt. K* ist in Einheiten von 9 uj R angegeben. 

Die Zerlegung von K* in seine x- und y-Komponente liefert Tabelle 7. 

Fiir die resultierende Gesamtkraft des Voith-Schneider- Propellers (einschl. der eben betrachte- 
ten Zusatzkraft) ergibt sich somit folgendes Ergebnis: 


in der Stellung 0 des Propellers: Kk‘) — — (67.03 SE 354 a ou, R, 


ouy R, 


K®) = —(25,88—50,9 = ous R; 
) : 
in der Stellung 7/12 des Propellers: K ] 


xX 
= — (6412 + 14,62 
x Uy) 
Ke a 


94:13 51.4 a o u2 R; 
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Tabelle 6. 


Se ee re 


Mittel kt A 4 2 OsNa - Mittelpunkt TAAL GraBe de Richt ms 
Ges Huger | ction re | Zosnateate K* | Zonmaurat ke | TEMES | zitlation P* | Zosatahratt K* | Zostera Ke 
ee a 2 ee eee eee 
0° +0,16 “9 2,6 — + 16,6° 180° +0,61 —% 9,6 ~2. 417752 
Ug Up Ug Uy 
ise +0,39 “0 642" aE 3 055° 195° +0,34 2 oe —166,7° 
> u 
Up Ug Ug 0 
30° +.0,56 2 9,9 0 4+ 43,3° 210° —0,04 22. 0,6 - 4 29,3° 
Uy Ug Ug 0 
v ° 
45° EO TO eo alma 2 a5 + 54,9° 225° =0,34=2 522 + 45,4 
Uy Uy Ug i) 
60° +0,80% | 141 +. 66,7° 240° 0,10 2 10,6—2 + 61,6° | 
Ug 0 Ug Ug | 
15° 1 0;88—% te al4,6—” 786° 255° Sa hie 1s5 9 + 80,0° wi 
Up Uy Uo Ug 
90° +0,95-2 | 15,7—2 + 91,22 270° = 120 18,24 + 92,7° all 
Ug Up Ug Up | 
105° Se OMee elo =" +104,9° 285° 123 oe +110,9° 
0 Ug Uy Uy 
120° 41,05 2 > 17,328 +115,4° 300° = 0s? 15,228 +128,6° 
0 Uo Uy 0 
5 Uo 6 ° vo Vo ° 
135 +1,04 “2 17,0 18 17 315 —0,80 2 Wee +146,3° I 
Ug Ug Up Up 
150° eo o7 ee 15,6 ~° 4+146,5° 330° Soy 6,6 —2 4-164,5° ay 
Up Ug Ug Uy | 
165° +0,82 28 13,0 —° +161,9° 345° —0,12 1,9 = —179,0° 
0 Ug 0 0 
Tabelle 7. 
Mittel- Mittel- Mittel- 
nkt d. * * unkt d. * * nkt d. * * 
PP liigel- & ey Zs Hiigel Ry Ky Piagel: Ke ES 
profiles profiles profiles 
| 
0° 2,49 = 0,742 | 120° | — 7,422 +15,62 —2 240° 5,04. 1 9,32 . 
0 0 Up 0 Uo 0 
15° 5 7 Vo ° Vo v fo} v 
71 3,402 | 135° | —11,31-2 | +12,69"2 | 255 2,69 +1527 
| 0 0 Uo Uo 0 0 
30° | 47,212 6.79 20) 150% de 13,00 2s eee ole = 270° | —0,86 +18,18 2 
Uy Ug Uo Ug Uy 
45° | 7,19 =! 10,23 8 165° 12,36 so" 285° 6,42 116,81 2 
0 0 0 0 Uy Uo 
60° 5,57 —2 12,942 | 180° | — 9,59 + 0,427 | 300° | —9,48 +11,88 2 
0 0 Up Uo 0 Uy 
75° | +2,92 “0 14,502 | 195° 5 Toe 1,22 | 315° 948-2 | 46,3322 
ug Ug Ug Ug Uo Ug 
90°") =0,33 — 15,70 =2 | 210° |; 4. 0,52 28. 10.99.28) e30°) |} 6 36 ents 
Uy Ug Ug Uy Uy Uy 
105°) 4,298 | 16,1328 | 295° | 43,6525 | 4 3.99 2 | g4geuml eeeteog <meta 
Ug 0 Uy Uy u 
| 0 0 


in der Stellung 7/6 des Propellers: 


Ke) =—(64, 31 + 12,8) 9 u R, 


KY) = ee p u2 R; 
in der Stellung 27/4 des Propellers: K() = (6%, yeaa 2 
Uo Ona 


KY?) = — (2, 92 — 50,9 2 


mM 
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Fir die gesuchte genaue Vorwartsfahrt (resultierende Gesamtkraft entgegengesetzt zur Anstré- 
mung) mu8 K(?)/K(?) = vp/uy sein. Diese Bedingung ergibt 


in der Stellung 0 des Propellers: = =(,21, arcte eS ROR 
u u 


0 0 
in der Stellung w/12 des Propellers: a 0,20, arctg ae THe 
ut 0 


0 u 


in der Stellung 7/6 des Propellers: —2 = (el arctg —° ==) ORR, 
0 Ug 


in der Stellung z/4 des Propellers: 2 = 0,19, aretg i OH 
Ug 


0 


Als Mittelwert erhalten wir a = 0,195, aretg Se TaOe, 
0 Uy 


Legen wir nun eine Anstrémung mit v,)/uy = 0,195 zugrunde, so wird 


in der Stellung 0 des Propellers: K(?) = — 69,6 0 ujR, KP) = —15,9 9 uR, KY) = 7140 uk, 
Richtung: Bx =— 167°; 

in der Stellung 7/12 des Propellers: K‘?) = — 67,0 9 ujR, K?)=—14,1 0 uj R, K®)= 68,4 0 up R, 
Richtung: Bx = — 168°; 

in der Stellung 7/6 des Propellers: K?) = — 66,8 0 ujR, KY)=— 11,7 9 upR, K®) = 67,8 9 uy R, 
Richtung: Bx = — 170°; 

in der Stellung 7/4 des Propellers: K?)—— 70,1 9 uzR, KY) = — 13,0 pug R, K”) = 71,3 9 up R, 
Richtung: Bx = — 170°. 


Zuriickkommend auf die Fragestellung am Anfang von Abschnitt 9 kénnen wir jetzt das Er- 
gebnis formulieren: Bei einer Anstrémung in x-Richtung mu8 der Anstellwinkel Null der Propeller- 
fliigel etwa an den Stellen 79° und 259° liegen, wenn genaue Vorwartsfahrt erreicht werden soll. 
Dabei schwankt die Vortriebsrichtung um etwa + 1,5° um den Mittelwert, je nach der momentanen 
Stellung des Propellers. Die GréBe der resultierenden Kraft schwankt um etwa 4.3% um den 
Mittelwert. 


(Eingegangen am 5. Februar 1955.) 
Anschrift des Verfassers: Dr.-Ing. W. Isay, Berlin-Dahlem, Schweinfurthstr. 90. 
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The Equations of Motion for Curved Elastic Bars deduced by the use of 
the ,.Method of Internal Constraints“ @ 


By E. Volterra 


1. Introduction. The equations of motion for curved bars are derived taking into account the 
influence of shear and of rotatory inertia by assuming that during motion the sections originally 
normal to the axis of the bar remain plane. 

The elastic isotropic homogeneous bar is supposed to be an elongated solid with a curvilinear 
baricentric axis. It will be referred to a fixed system of orthogonal cartesian coordinates (QE577,6) 
and to a system of curvilinear right-handed coordinate system (0, x, y, 2) in which the origin 0 
coincides with a generic point on the central axis of the bar, x with the central axis (and has as positive 
direction the positive direction of the central axis), y coincides with the principal normal (7) and is 
directed positively towards the center of curvature, 


ff 
| and z coincides with the binormal (6) at 0 (see fig. 1). 

The condition that the sections originally nor- 
mal to the axis of the bar remain plane during 
motion implies that the elastic displacement vector 


V at a generic point P of the bar depends linearly 
only on the coordinates y and z or that the following 
vector equation be verified: 


V(x, ¥, 2 t) = 9(x,t) + yA(x,t) + 2 u(x,t) (1) 
where v, A and i are three vectors, a priori unknown, 
functions of the variables x and t. 

This method of analysis, which was used more 
than twenty years ago to study static deformations 
of curved elastic bars!, will be referred to as the 
,,Method of Internal Constraints‘, since the funda- 
mental hypothesis on which the method is based 
rests on the assumption that the elastic displacement 
of the bar occurs under the special condition of internal constraint expressed by equation (1). 

The geometric significance of the internal constraint will be first discussed, then the expression 
for the square of the element of length (ds?) will be derived in the more general case in which the 
central axis of the bar has both first and second curvatures (which will be denoted respectively by 
the letters c and rt). Thus, from the fundamental (Euclidean) metric tensor it is possible to calculate 
the components of the strain and then the expressions for the Potential and Kinetic Energies of the 
bar in terms of the nine (unknown) components 9x, 1’), 25 Ax Ays/z5 [xs [ys [lz of the three vectors 7, 
A and jl. 

The application of Hamilton’s Principle furnishes the nine differential equations and the boun- 
dary conditions which univocally determine the nine unknown functions 7,, Vy, 3 Axs AysAgy Wen Ute 
in the case of free vibrations. 

The case in which external forces are acting on the bar can be dealt with by introducing into the 
variational principle the work done by these external forces. 

Although the method of analysis presented in this paper applies to the general case in which 
the curvature of the central axis and the inertia characteristics of the cross-section of the bar are 
known functions of the variable x, for simplicity sake, the equations of motion will be derived only 
for the particular case in which the curvature and cross-section area of the bar are constant. 


ivr 


Fig. 1. 


2. Geometric significance of the internal constraints. Consider a generic element P of the bar in 

the natural (undeformed) state and call x, y, zits coordinates, while call x’, y’, 2’ its coordinates in 
iste Wot ; 

a new position P’ when deformed, and u, v, w, the components of the elastic displacement along the 


* See Enrico Volterra, Verhandlungen des Internationalen Mathematiker-Koner iri 
: r » Verhandh - greB Zurich 1932 II. Band, 
p. 247; Rendiconti Accademia Lincei 16 (1932) p. 220; 20 (1934) p. 424, 463; 21 (1935) p- 14; 23 (1936) p. 329, 
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x, y and z axes. From equations (1) it follows that 


K= 4 -SU =*i , [x % Lan Ys 
y=ytu=y+r, by B+ Ay ¥ 5 (2) 


, 


2 =z +w=2z+y»p, fs 2 +A, ¥ 


The passage from the position P(x,y,z) to the new position P’(x’,y’,z’) can be imagined decom- 
posed in two successive transformations: 
1) A first transformation T, which brings point P to an auxiliary point P. ACGME REA 
2) A second onct ecineon, T, which brings point P, to its final position P’. 
These two transformations are defined by the isilowmne equations: 


i tee 3 (a), 

Di Oa a ye Bye Ay y (b), (3) 
By Ve a AY (c), 
BB Vg dh A Vs 

Tey = Yin (4) 
Bete ibs 


In fact, by substititing into equations (4) instead of x,, y,, z, their values given by equations (3) 
and observing that the quantities W2,, Ay Vy, Vy(%1)3 Ms2% Az¥y, V(%,) differ from the quantities 
My 2, Ay Vo Vy(X)3 Uz23AzY, V2(x) of very small quantities, which in the present order of approximation 
may be neglected, equations (2) are obtained. 

Let us analyze separately the two transformations T, and T,. Since in the transformation T,: 
x = x,, it follows that the transformation T, operates in a generic plane normal to the axis of 
the bar. Accordingly in this transformation not only does every normal section of the bar remain 
plane, but, moreover each normal section transforms into itself, according to formulae (b) and (c ) of 
equations (3). Transformation T, is a linear transformation whose Soe eae significance is an 
affinity since the axes at infinity arena to each other. 

As it is well known, its geometric interpretation is a pure deformation accompanied by a rigid 
motion. In fact, by using the eS notation: 


[= [pte a + (A, + fly) x ¥ + Ayy*] , 


ay =Vyy5 “i, 


equations (b) and (c) of equations (3) can be written under the form 


df 
M1 Y= dy + O21 Gy, 

df (6) 
ee Se I che OE 


where a, and a, are the components in the y z plane of a vector 4, while w is the component in the 
x direction of a vector © normal to the y z plane. 

If the deformation is infinitesimal, then the three terms of the right-hand side of equations (6) 
may he considered as the result of three successive transformations: 


1) A pure deformation which displaces point P to point P’ defined by 


(P’ — 0) = (P—0) + grad. (2) 
2) A rotation which displaces point P’ to point P” defined by 
(P’ — 0) =(P’— 0) 4+ @x (P’— 0). (8) 
3) A translation which displaces point P” to point P, defined by 
(Bye Opies (BS — 0) pa’, (9) 


In fact, by substituting equations (7) and (8) into equation (9) the former equations (6) are 
obtained. 
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As a consequence of the transformation T, defined by equation (4), each longitudinal fiber of 
the bar (i. e. parallel to the x axis) is subjected to a rigid displacement along itself, the amplitude 
of the displacement «’ — x, varying from fiber to fiber. 

It thus follows from equations (4) that 


—— x, — An N41 + [x 21 -+ Vx (10) 


and that 
x’ — x,=0 
when 


Ax Vy xt +x = 0. (11) 


With reference to the 7, plane normal to the x axis of the bar (z, being the original plane x 
normal to the axis of the bar and which has been transformed in accordance with the T, trans- 
formation) the fibers parallel to the x axis (and therefore normal to the z, plane) can be identified 
by their intersection Q with the z, plane. 

Since equations (10) can be written as 


= Peace i i Ax y 
hi KV ee ae psa : ‘ 12 
eee Wek roe 
it follows that the above displacement is proportional to the distance of point Q from the axis r of 
the zr, plane, the axis r being defined by equation (11). 
Due to the fact that the displacement 


x’ — xX = Vx ot Ux] ae Au Vy (10) 


in the plane 7, is supposed to be infinitesimal, it can be considered as a rigid rotation around the 
r axis in the same plane. 

Since the transformation T, operates in the normal sections of the bar by transforming into 
themselves, while the T, transformation produces a rotation of the transformed sections around 
the r axes, it follows that as a consequence of the deformation produced by the constraints and 
expressed by equations (1), or the analogous transformations T, and T, represented by equations (3) 
and (4), the sections normal to the axis of the bar remain plane. 


3. Expressions of ds? Components of strain in orthogonal curvilinear coordinates. In order to 
derive the expression of the square of the line element ds? in the orthogonal curvilinear system 
X,Y, 2, consider a generic point P of the bar of coordinates x, y, and z. Denoting by 0 the inter- 
section ofthe normal from P to the central axis of the bar (axis x) with this x axis (see fig. 1) call 
t, n and b the unit vectors of the tangent normal and binormal to the central line at point 0. 

A generic displacement dP of point P referred to the fixed coordinates system (2,&,7,C) can 
always be considered as the resultant of the two displacements dP, and dP,, where the first dis- 
placement dP, is relative to the mobile coordinate system (0, x, y, z), while the second displace- 
ment dP, is produced by the displacement of point P considered rigidly connected with the coordi- 
nate system (0, x, y, 2). 

Since 


P—0=yn+2z6 (13) 
and due to the fact that 0 7 and 6 in the first displacement dP, have to be considered as constants, 
it follows that 

dP,=dyn+dzb. (14) 


The quantity dxt represents the displacement of the point 0 [origin of the mobile system of 
coordinates (0, x, y, z)] as a consequence of the displacement of point P, 
The instantaneous rigid motion of the system (0, x, y, z) (by assuming 0 as instantaneous center) 


can be imagined characterized by an elementary translation dxt and by an elementary rotation 
dx W, where, in view of Frenets formulae 


®O=—ttten. (15) 
From the fundamental formulae of kinematics of rigid motion, it follows that 


dP, = dx [t +x (P—0)]. (16) 
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By substituting for @ and (P— 0) the expressions given by equations (13) and (15) one obtains 
dP=dP,+dP,=t(l— yc)dx +n (dy + 21 dx) +b (dz— yr dx). (17) 
Since the expression for ds? is given by 
aR ES i te (18) 
From equations (17) and (18) the expression for ds? 
ds? = dP. dP = dx* [(l— yc)? + (yt)? + 2(r)?] + dy? + d2&—2yrdxdz+2zrdxdy (19) 


follows. If it is supposed that, in a first approximation, the second powers of the quantities (pure 
numbers) y c, zc, y t and zt can be neglected, equation (19) becomes 


ds? = dx? (l1— 2 yc) + dy? + d2—2yrdxdz+2zrdxdy. (20) 


The fundamental tensor will in this case be expressed by the metrics! 


a,,=1—2ye, GS ak T > ae al 
Ay = 22T, a= 1, dy, = 0, (21) 
a3) =— 2yT, a3,= 0, da, = 1. | 
From which it follows with the same degree of approximation that 
a=1—2yc, ja=1—ye, is Saar (22) 


Ja 
In the particular case in which the central axis of the bar is supposed to be a plane curve t = 0 and 
equations (20), (21) and (22) become 


ds? = (1— 2yc) dx? + dy? + d#, (23) 
a,=—1—2yc, a,—0, Gj, — 0, | 
a,—0, A,=—1, a, = 0, (24) 
a3,=0, dag = 0, a,—1, 
al {I 
a=1—2yc, ja=1—ye, a= Sey es 
Vay 
The components of strain are calculated from the formulae? 
es 0 U; ag 1 < 0a;; Ut, (25) 
(=F) 2 aj; ol, OX}, Vay, ¥ 
1 @ if hp 0 / uj )| 
fe aes |r, = aj es 24 
2 Va; 4;; | 55 (Tex) # (7 (28) 
(uj =u, pyc tr, LE Ce SE C26 TSS Meg X3 = 2). 
It follows that 
OA, Gib, LOD. 
c= (1+90(15, aa 52) — ery + py + Vy) , 
éy Ay 5 
Ez = Uz» 
OA, OM, , Ov 
@=t&y + by + & = (1 + ¥c) (va: z Ox + 52) —6 yy Hy% t Vy) 4 hy + Ms ° (27) 


0 
Yay = (1+ 9) (9 GEE GEL Be) te (bay + bat +) +e 


02, OL, , Wz 
Va, = (1+ ye) (y Ox Sina +5) + Hes 


Vys = My +A: + 


1 The notation used in the text is the same as the one used by T’.. Levi-Civita in his book: The Absolute 
Differential Calculus, London 1926. | 
2 See A. E. H. Love, The Mathematical Theory of Elasticity pp. 51 — 58, 4th ed., Cambridge 1934. 
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4. The equations of motion. The Total Potential Energy of the bar is expressed by 


Revi, 2 G 2 2 2 
=| ya [rte +6 (ete +e) +2 vty +e t vid] de dy de , (28) 


where V is the volume of the bar, Ais Lames constant and G the modulus of rigidity of the material 
of the bar related to Youngs modulus E and to Poissons ratio y by the relationships 
vE E 


ASC py Oty alee) 


By substituting equations (27) into equation (28) and performing the integrations with respect to 
the variables y and z, calling 


Basti = 2 
pees ios hi belee ie cA 
A A 
the principal moments of inertia of the cross-section of the bar and A the area of the cross-section 
of the bar. 
Since 


[ 9 dd= fo2dA=0 
A A 


one obtains the following expression for U: 
vibe git 9 Bf eda 
0 
5+ G AY +51 (2) +552) +24 +74 
(ass we} ae ep + A (ua)? +S Ay)? 


() 
a Pa ee Ov, 


Os 
a 
& 
2 


AQ) yay peer es =4G6Apyd, (29) 


OA, 0 
+ (A+ 26) I, 52s Up ie mame 
AA Chis ¥y— (A+ 26) € A vy a NOT cei eee ieee EEC her 


Ov, OA, Ovy 


0 
Cd eC Le rei 28 Cte 


The Total Kinetic Energy of the bar is expressed by 
0 — |/du\2 dv\2 Ow \2 
T= {Va (7) + (] ! (Fe) | ae ay a (30) 


2 
V 
where g is the density of the material of the bar. By substituting into equation (30) the expressions 
for the components of the elastic displacement u, v, and w, by performing the integrations over 


the cross-section of the bar, one obtains for the Total Kinetic Energy of the bar the following 
expression; 


l 
Oey Na (eAe\ Oly Ov,\2 Ov, Ohy OA 0 
ee ee € Jy +5(3 rich Ale) 2 Ize = +E fy + 45(Fe) +4(R) 
0 


Avy OA an,\2 a 
2 We Fae + LF tte F) + AG ae Me 0m | a. 


(31) 


Ot ot 


On the assumption that no external forces (body and surface forces) are applied to the bar, and 
supposing also, for simplicity sake, that the curvature c of the axis of the bar and the inertia cha- 
racteristics of the cross-section are constant, by applying Hamiltons Principle under the form 


df (U+T) dt =0 
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to the whole bar, supposing the displacements to be zero at t, and ty, 


nine equations of motion in the nine unkno wn functions 9 Aes > hes 5 29 Vy 7) Vos Vz of x and ts 
hes Ay bi by Us x 4 


(A + 26) Ste CH + (A +26) eS (A436) er 1 AOU Ye. 
2 Oh, av, : 


2 


= OTe OTC an 


OA 2 
Ce sae eo (A 8G) yee ei Got, Gri cS = 0 one, 


Gr oo Qe 26y,. ae Ore sey} ee 4 Dery +6 ro? 


9 OAy O*vy 
= rt eg teh, (32) 


ru, 
Cry et — Ady —(A + 26), — ASE + Levy = 18 Me, 


e OPA, 
an +A + 26) es — Ge 2% (A 4 26) 02% Cea, 


uy 
Cry 53 


Pee? Ohne 


Guy—GA.+ (A+ 36) ¢ a 
iS (33) 
lee ets 


Pca... 2 


(+36) or9 ME 9 1 Mh 


A,  o%, » a, 
(aa eo, 18 os 


—c1,\0/A,=90, 


OV x 
x 
On 
cae + er, 4 Ja] Ory = 0, 
=) 


9 


| 
| 
Fe tog t+ ets] d2y = 9, (3%) 
| 
| 


Geezer (ce Cyt —(A+ 26) evn, + Ady + Ap] dm = 0, 


Bs. : ah A= 0, 
ea eUs| 0 My = 9» 
[os 


mie Oe 


re 
| 


[ee tae + ores] =— ‘5%. =0. 


In these equations appear together with the nine unknown functions hie Ayo Aas pecs [bys far Yxe Vy and 
y, the radii of inertia r, and r, of the cross-section of the bar, and the curvature, c of the axis 
of the bar which have been supposed, to be constant. However, the more general case in which 
the curvature of the axis of the bar and the inertia characteristics of the cross-section are functions 
of the variable x can be dealt with by the use of the same variational principle. As equations (32), 
(33) (34) and (35) show, during free vibrations of a curved bar the functions /,,/,, Us, Vx and vy; 
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“A, lx, [ly and y, are coupled. In the case of straight rods c = 0 and the above equations become 
es Ar Z 0A. 0’y OA, O*vy 
(A+ 26) r: ae OT He” Car te a, =o ae ey) 
OA Ov, 0rA 
Griag — (A+ 26)Ay i= A Ms Olea > 
0? Ov, aun 
Gr} oe — Ady —(A+ 26) ts —A 5 = 0 Tae» (37) 
A a) Ov, 
(£26) Se + Sy 4 als 9 Se, 
OA, OA, a? au 
Grea. =— Gi, Guy = Trae Gry aq — Gy —CA,= Ory ae 2 (38) | 
iy Cb Ov, OU Ov, OMe 0”, 
(A+ 26) 15 "Y Ox? Cus G One: ry ow” G Ox? | Ox Ot? (39) 
with the following boundary conditions at x = 0, x =I: | 
gx] 34s = 9. GE + Ae] y= 0, (40) | 
aa OM. Ov. 
Fe |r = 0 [Se] om = 0. G | 2G) = tAAy+Apz|d%x = 0, (41) 
Ou OA, 
ge] Ou = 9 [as |o% =o. (42) 
OM. av, 
pe] Ou = 0 ge + ae] do = 0. (43) 


In the case of plane strain, the components u and v of the elastic displacement in the x and y coor- 
dinates will be given by | 


ae A v= Vy + Vy > (44) | 


where /,, A), vx and vy are four functions of the variables x and ¢ (respectively components of the 


vectors 2 and y on the « and y coordinates). x coincides in this case with the central (baricentric) 
curvilinear axis if the bar, while y coincides with the normal to the central axis of the bar and is 
directed positively towards the center of curvature. Onthe assumption of constant curvature c and 
inertia characteristics of the bar, one obtains the following four equations of motion in the four un- 
known functions /,, Ay, 7, and v, of x and t: 


7A, o2y E 01, E E ov IBGE 0” 
peed 2 LER — 2 x ¥ sh 2 x 
Or e8 wT 0 atone (= 98) Gxt Lo) ae EE Oe ony ghea eee ts 
Ec(3—») dy come 
D (Da), de Ue) ae 


Gye Oy AEN yy Bele ES og 
Goon 0 8 oe — D9) Gah 1 2 (Depp) Ga (Lg) oe ecu ere 
i =?) Eno, 
2(1—*) ° Ox 
(ey, ote OL hve ty BRA eu, .) 
0 OR oot  —*) Ox? " 2) Ox 0)" ° Ox 
GB—r”)E cae Ec 2 
2 (1 8) © HT 
ar a E aA E ine Ge TO 
9 Ay , sree 2 Phy v Ue Ne eee ie 
Rt Oi Wane bu Ort ay (peg) a Oa 91) (191) dx 209) oy 
Ey E 0’ 
Se | re x 
F(a) oo eee es J 
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with the following boundary conditions at « = 0, x =I: 
OA, Ov, - 
ie doe + ¢hy] 02s =the 
Ovy A 2 Ody es 
lar +4 +r rain cn] 3 =0,. 
Ov, OA, I) 
Ee vay 4 roa — cry] bre =O, 
Ody Ovy © = 
age Che] Oy = ue 
In the case of a straight bar c = 0 and the above equations become 
polls. HS 0, Ess E avy 
Or aa i (T Sayt) dx? 21+)" 2(1+7) | i 
Ory E dy Ea, | mp) 
Or 2) os LID) on > 
Ce E 0%», Ev dhy 
Cae ~ A —*) dx? | A — 2%) Ox? i 
shy E Oly E_, Ev ed 
Cr at 20+) ont” 1 — 08) A — 0) Ox 
with the following boundary conditions at x = 0, x = I: 
A 
Fe yea 
: (49) 
yy ah, 
lez +4.| O70 
OV» a 
= +9 iy] r =0, 
a (50) 
ia DA 20 
x 


(Eingegangen am 28. Februar 1955.) 


Anschrift des Verfassers: Prof. Enrico Volterra, Troy, New York (USA) Rensselaer Polytechnic Institute. ! 


1 The material presented in this paper is based on an investigation which is being conducted under 
the sponsorship of the Office of Naval Research, Department of the U.S. Navy, Washington, D. C.; 
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A One-Dimensional Theory of Wave-Propagation in Elastic Rods Based 
on the ,,Method of Internal Constraints“ 


By E. Volterra 


1. Introduction. The paper presents a one-dimensional theory of wave propagation in elastic rods 
which takes into account the effects of shear and of rotatory inertia and which unifies a number of 
separate engineering treatments of the problem. 

The equations of motion which were derived in a previous paper! will be used. These equations 
were obtained by assuming that during motion the sections originally normal to the axis of the bar 
(x axis) remain plane, or that the following vector equation be verified: 


V(x, ¥, 23 t) = V(x, t) + y A(x, t) + 2 B(x, 1). (1) 


In equation (1) V represents the elastic displacement vector at a generic point P of the bar,y,/ 
and wu are three vectors, a priori unknown, functions of the coordinate x and ofthe time t. The bar is 
referred to a system of orthogonal cartesian coordinates (0, x, y, z) where x coincides with the axis of 
the bar, while y and z with the principal axes of inertia of the cross-section of the bar. The equations 
of motion which, in the case of a straight bar of infinite length, univocally determine the nine func- 


LIONS Vx, V yy Vz3AxyAysAc3 xs Wy, Uz (components of the vectors y,/ and pu on the axes x, y 2) are: 


Upie noe cL oe noe | ica (2) 
He re ermeee  e e ee o o (3) 
Ore? Otis ity Gr oe Cu CL on 
Crone yee er i ae 
CF ay ae = ES (5) 
weno fy eat nade eM 


In the case of plane strain the equations which determine the four functions /,,A,, 7, and y, are: 


» Px Oye E Ey 
r= = r ana Ae eon > 
ot (l—»?) ax? 2(1 + ») 2(1 +») Ox 
Bing, ER Coy el Oe (6) 
Cae 20+) ox? | 2040) Ox’ | 
Ox a E 02x | E dhy 
Ca ~ I —r) ox? | G—v*) dx’ | 
fGly Ee Oy E Eas (2 


hy 


uf = 2 : a eee ee 
ero 20+) o2  —) (l—?) Ox * 


In the above equations / is Lames constant and G the modulus of rigidity of the material of the 
bar related to Youngs modulus E and to Poissons ratio y by the relationships 
= vE om E 
~ (L+7)(l—2y) ’ ~ 21+)" 


While ry, 7, and r are the radii of gyration of the cross-section of the bar, g is the density of the 
material of the bar. 


A 


1 See: Ing.-Arch. 23 (1955) p. 402. 
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These equations admit a direct treatment of the problem of wave propagation in straight bars of 
infinite lengths. In fact it is possible to derive from them the frequency equations for flexural, tor- 
sional and compressional waves. 

The results for flexural waves will be compared with those obtained by the Elementary!, Ray - 
leighs®, and Timoshenkos® one-dimensional theories and with those obtained from the exact Poch- 
hammer solution’. In the case of compressional waves it will be shown that, in the particular case 


of cylindrical rods, the solution of the problem reduces to the same solution as the one obtained by 
Mindlin and Herrmann’. ; 


2. Flexural Waves. Let the components of the displacement of the bar be 


Ce ¥ A(x, t) ’ v= Vy(X, t) > Oi — 0 (8) 
or 


u's Ufx, 1), = 0, 0 —— (a,b) A (9) 

The motion of the bar will be characterized by equations (2) or (3) respectively or, in the case of 

plane strain by equations (6). Writing y instead of, or y,, r instead of ry or r, equations (2) and (3) 
can be written as 

E (1 —) ae 3—4y\ oly 2(1 et) 2 ary 

e (l+vr)(l—27) Ox4 1 — 2] 0x? dt? E ots 


O*y 
e Ae 0, (10) 


while equations (6) can be written as 


E Te Oty ofs—\ ot 2 (Mi) nee Yn ony 
0 (l—»*) oxt ” ( = dx? Or Eo" OT Oe 


=0. (11) 


Let us compare the results given by equations (10) and(11) with the results given by other one- 
dimensional theories of flexural wave propagation. 

The simplest theory of flexural vibrations (Elementary Theory) assumes that the displacement 
of an element of the bar consists solely of translation perpendicular to the x axis. The corresponding 
equation of motion is 

E 2 dty | dy 
a ORE ot? 


= 0. (12) 


When the elements of the bar are considered to undergo rotation (without distortion) the term 
— r* dty/dx? dt? must be added to equation (12) to take into account the effect of the rotatory 
inertia of the element. The differential equation (12) is modified to the form usually attributed to 
Rayleigh (Rayleighs Theory) 

ee, 02%, ees: o? 

LS SEL ITN I ae (13) 


When the shearing, the rotation and the lateral displacements of the elements of the bar are 
taken into account, the differential equation takes the form given by Timoshenko in 1921 (Timo- 


shenkos Theory)? 


Q ) 
Prescott’ in 1942 has derived equation (14) by a different method from that used by Timoshenko. 
The constant K (Timoshenkos Shear Coefficient) which appears in equation (14) is defined as 
the ratio of the average shear stress on the cross-section to the product of the shear modulus and 
the angle of shear at the neutral axis and depends both on the shape of the cross-section of the bar 
and the frequency of vibration. 
Professor R. D. Mindlin’ derived for rectangular cross-sections the value of K = 7?/12 = 0,822. 
Recently Professors R. D. Mindlin and K. Deresiewicz® have determined the coefficient K for the 


om t EK ae Fn 


Ph 2(1++y)] a 21 +7) Oty dy 
fee | ae. ‘la y or 0. (14) 


1 See H. Kolsky, Stress Waves in Solids, Oxford 1953. ; 

2 See Lord Rayleigh, The Theory of Sound, New York, Vol. I, Chapter VIII. 

3 See S. P. Timoshenko, Philos. Mag. 41 (1921) p. 744. 

4 See A. E. H. Love, The Mathematical Theory of Elasticity paragraphs 199—202 4th ed., Cambridge 1934. 

5 See R. D. Mindlin and G. Herrmann, Proceedings First National Congress Applied Mechanics Chicago 
1951 p. 187. 

6 See J. Prescott, Philos. Mag. 23 (1942) p. 703. 

7 See R. D. Mindlin, J. Appl. Mech. 18 (1951) p. 31. 

8 See R.D. Mindlin and H. Deresiewicz, Office of Naval Research, Project NR—064—388, Technical 
Report No. 10, June 1953. 
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following cross-sections: circle, ellipse, orthogonal parabolas and a variety of ovaloids and have 
shown that the values of K computed for these cross-sections all lie within about 10 percent of that 
for the rectangular cross-section. 

In the particular case of a circular cross-section R. D. Mindlin finds for K the value K= 0,847; 
while R. M. Davies! uses for the same cross-section the value K = 0,90. 

The problem of flexural vibrations in cylindrical rods of infinite length (Exact Theory) was first 
investigated in the year 1876 in terms of the general equations of the Mathematical Theory of 
Elasticity by L. Pochhammer and independently in 1889 by C. Chree. An outline of this exact 
theory is discussed by A. E. H. Love. 

The frequency equation for flexural vibrations which results from this theory was given in deter- 
minant form in 1941 by D. Bancroft? who suggested that it was unlikely to be useful owing to its 
complexity. However, in 1943 G. E. Hudson® was able to compute from this equation the values of 
the phase velocities of flexural waves in a cylindrical bar. Moreover, Hudson showed that the fre- 
quency equation of flexural waves has only one root and therefore that flexural waves in cylindrical 
bars unlike compressional and torsional waves can be propagated in only one mode. Hudsons 
results for different values of Poissons ratio y are presented in tabular form for the ratio of phase 
velocities of flexural waves to the velocity (ce a VGio) of distortional waves in an unbounded 
medium versus the ratio between the wave-length of distortion waves and the circumference of the 
cross-section of the bar. 

In 1948 R. M. Davies by interpolating the values given by Hudson for the case of Poissons ratio 
y = 0,29, gave in graphical form the relationships of the ratio c/c, of the phase velocity to the velo- 
city of longitudinal waves of infinite wave-length and of the ratio c,/cy of the group velocity to the 
velocity of longitudinal waves of infinite wave-length of flexural waves in a cylindrical bar versus 
the ratio a/L of the radius of the cross-section of the bar and the wave-length. 


3. The Phase-Velocity of Flexural Waves. Let 
y(a,t) =A ee e9) (15) 
where A is a constant, y = 27/L, L the wave-length, c the phasevelocity and i = ay By 


substituting equation (15) into equation (12) and (13) one finds the following expressions for the 
phase-velocity: In the Elementary Theory: 


¢. are 
Gaur ebae (16) 
in the Rayleighs Theory 
c 1 
ee oc aa (17) 


where 


By substituting equation (15) into equations (14), (10) and (11) one finds the following frequency 
equations: In the Timoshenkos Theory 
K + 40° 8 [K +2(149)] K 
4 2 25, 
aT 8 (1 + v) x? 6? xy f a1 + v) O: (18) 


in the Theory based on the ,,Method of Internal Constraints“: 
a) in the case of three-dimensional elasticity 
4 (l—27)+4070(3—49) , l—y 
x 


*, 8(1—2rv)Q+rn259 “2 Sa py an = 9 (19) 


b) in the case of two-dimensional elasticity 


hie 
1+ 472 0? ( = 4 ; (20) 
Ls = x | = 0 
8 (lp) ade 3) BT 98 : 


1 See R. M. Davies, Philos. Trans. Royal S 

'. A : yal Soc. London 240 (1948) p. 375. 
2 See D. Bancroft, Physical Review 59 (1941) p. 588. Ps 
’ See G. E. Hudson, Physical Review 63 (1943) p. 46. 
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In equations (18), (19) and (20) 


x= au OS Tr 

The frequency equations (19) and (20) like the frequency equation (18) are quadratic equations in 
x° == (c/c9)? the roots of which are real, the greater giving values of x greater than unity (and there- 
fore corresponding values of ¢ > cy) while the smaller roots giving values of x less than unity (and 
therefore corresponding values of ¢ < cy). Since the larger roots do not appear to have any physical 
significance, they will be ignored. In fact, as already pointed out in the previous paragraph, Hudson 
has shown that the root x of the flexural waves in the Exact Theory is a single valued function of 6. 

From equation (18) one sees that for large values of 6, x, approaches asymptotically the value! 


bee as (pa 
m= eee (21) 


While from equations (19) and (20) one sees that for large values of 6 the roots x, and x, approach 
asymptotically the values 


a ees / i 
Xy = — > — 
cee 2(L-F9)’ se 


ace re 
MCL see ch aa Vaa=s . (23) 
In Tables 1, 2 and 3 the values of the roots x,, x, and x, of equations (18), (19) and (20) are given 


for the following values of the parameters y and 6: 


P— OL0> 0,055 0.105 045:,0,20* 0:25> 0:30: 0.352 .0.40', 
6 = 0.01; 0.05; 0.10; 0.50; 0.1.00; 2.00; 3.00; 6.00; 10.00; co. 


Table 1. Phase velocities for flexural waves according to Timoshenkos Theory (i: = 3) : 


ae 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
0.01 0.063 0.063 0.063 0.063 0.063 0.063 0.063 0.063 0.063 
0.05 (OAPATL 0.2740 0.2729 0.2719 0.2708 0.2698 0.2687 0.2678 0.2667 
0.10 0.4318 0.4278 0.4240 0.4201 0.4165 0.4129 0.4093 0.4059 0.4025 
0.50 0.6425 0.6287 0.6157 0.6034 0.5917 0.5807 0.5702 0.5603 0.5509 
1.00 0.6601 | 0.6447 0.6303 0.6168 0.6041 0.5922 0.5809 0.5703 0.5602 
2.00 0.6650 0.6491 0.6343 0.6204 0.6075 0.5952 0.5837 0.5729 0.5626 
3.00 0.6659 0.6499 0.6350 0.6211 0.6081 0.5958 0.5843 0.5734 0.5631 
6.00 0.6665 0.6504 0.6355 0.6215 0.6084 0.5962 0.5846 OPSTSt 0.5634 

10.00 0.6666 0.6505 0.6356 0.6216 0.6085 0.5962 0.5847 OF iiom 0.5634 
oO 0.6667 0.6506 0.6356 0.6217 0.6086 0.5963 0.5847 0.5738 0.5634 


Table 2. Phase velocities for flexural waves according to Three-Dimensional Theory of Internal Constraints. 


ee 0.00 0.05 0.10 | 0.15 0.20 0.25 0.30 0.35 0.40 
0.01 0.063 0.065 0.065 0.066 0.066 0.068 0.072 0.079 0.091 
0.05 0.2776 0.2772 0.2782 0.2809 0.2858 0.2936 0.3059 0.3251 0.3574 
0.10 0.4409 0.4380 0.4367 0.4373 0.4402 0.4458 0.4553 0.4698 0.4919 
0.50 0.6762 0.6624 0.6497 0.6381 0.6273 0.6254 0.6082 0.5996 0.5917 
1.00 0.6985 0.6825 0.6676 0.6537 0.6407 0.6285 0.6171 0.6063 0.5961 
2.00 0.7049 0.6881 0.6720 0.6579 0.6443 0.6324 0.6194 0.6080 0.5972 
3.00 0.7061 0.6892 0.6734 0.6587 0.6450 0.6324 0.6198 0.6083 0.5974 
6.00 0.7068 0.6898 0.6740 0.6592 0.6454 0.6324 | 0.6201 0.6085 0.5976 

10.00 0.7069 0.6900 0.6741 0.6593 0.6455 0.6324 0.6201 0.6086 0.5976 
oe) 0.7071 0.6901 0.6742 0.6594 0.6455 0.6324 0.6202 0.6086 0.5976 


In computing Table 1 the value K = 8/9 has been assumed for Timoshenkos Shear Coefficient. 
This is the value given by Professor S. P. Timoshenko® using the experimental values of LN. G. 
Filon®. In figure 1 the ratio c/cy between phase velocities of flexural waves in cylindrical bars and 


1 See R. M. Davies, Philos. Trans. Royal Soc. London 240 (1948) p. 375. 
2 See S. P. Timoshenko, Philos. Mag. 43 (1922) p. 125. 
3 See A. FE. H. Love, L.c. paragraph 245, 
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Table 3. Phase velocities for flexural waves according to Two-Dimensional Theory of Internal Constraints. 


ie 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
0.01 0.063 0.067 0.071 0.074 0.077 0.084 0.087 0.095 0.105 
0.05 0.2776 0.2911 0.3060 0.3225 0.3409 0.3615 0.3848 0.4113 0.4428 
0.10 0.4409 0.4600 0.4809 0.5037 0.5289 0.5568 0.5878 0.6268 0.6619 
0.50 0.6762 0.6963 0.7178 0.7410 0.7661 0.7927 0.8236 0.8568 0.8940 
1.00 0.6985 0.7175 0.7379 0.7600 0.7840 0.8104 0.8395 0.8718 0.9080 
2.00 0.7049 0.7234 0.7434 0.7652 0.7889 0.8150 0.8437 0.8757 0.9116 
3.00 0.7061 0.7246 0.7445 0.7662 0.7897 0.8158 0.8445 0.8764 0.9123 
6.00 0.7069 0.7252 0.7451 0.7668 0.7904 0.8163 0.8450 0.8769 0.9127 

10.00 0.7070 0.7254 0.7453 0.7669 0.7905 0.8164 0.8451 0.8770 0.9129 
ee) 0.7071 0.7255 0.7454 0.7670 0.7906 0.8165 0.8452 0.8771 0.9129 


eeien 


Flexural Waves (Elementary Theory) 


Flexural Waves (Rayleigh Theory) 


Extensional Waves (First Mode) 


Flexural Waves 
(Method of Internal Constraints) 


Flexural Waves 
(Exact Theory) 


Eee 


Values trom 
Timoshenko’s Theory 


—x— k=0,9 (Davies) 
—o— k=0,847 (Mindlin) 


S|p 


0 04 O4 06 08 40 46 a4 46 48 2,0 
Fig. 1. Phase velocity of flexural waves in cylindrical bars {vy = 0,29). 


velocity of longitudinal waves of infinite wave-length are given versus the ratio a/L between radius 
of the cross-section of the bar and wave-length in the case in which y = 0,29. 

The curves obtained from the Elementary Theory [equation (16)], from the Rayleighs Theory 
[equation (17)], from the Timoshenkos Theory [equation (18)] and the results obtained from equation 
(19) as well as the curve for the first mode of extensional waves are compared with the results 
obtained by R. M. Davies interpolating the values given by Hudson from the Exact Theory 

In fig. 1 the results from Timoshenkos Theory are given as a series of points, the eS repre- 
senting the values of the roots x, computed by assuming K = 0,9 (according to R. M. Davies) while 
the squares representing the values of the roots x, computed by assuming K = 0,847 (according to 
R. D. Mindlin). As fig. 1 shows, for wave-lengths greater than about ten times ake radius of the bay 
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(a/L = 0,10) all four theories of flexural waves propagation give approximately the same results. 
For increasing values of the ratio a/L the values of c/c) given by the Elementary Theory and by the 


_ Rayleighs Theory are greater than those given by the Exact Theory. When 6— co the phase 


velocity becomes infinite in the Elementary Theory and equal to cy in the Rayleighs Theory whereas 
in the Exact Theory it becomes equal to c, = 0,5764 the velocity of Rayleighs surface waves. 
Instead in the Timoshenkos Theory the limiting values of x, = c/cy (for y= 0,29) are 


ec, = VK/2 (1 +) and therefore for K = 0,9: c/cy = 0,5906, for K = 0,847: ¢/cy = 0,5730. 
In the three-dimensional Theory based on the ,,Method of Internal Constraints“ the limiting 
value of x, = c/cy = 1/2 (1 +), and therefore c/cy = 0,6226. 


4, The Group-Velocity of Flexural Waves. The velocity with which the energy of a pulse is propa- 
gated is called its group velocity c, and is given in terms of the phase velocity c and the wave-length 


L by the equation 


dc 
¢,=c—L aL 
or by 
d d 
eg=e+d=c4 bas. (21) 
where 
r a 
= ah: Oy care [a 
From equation (21) it follows that the group velocity of flexural waves is: 
in the Elementary Theory: 
g=e+— pM = 20, (22) 
in the Rayleighs Theory: 
1 
ee 
ik A 7? r2 L2 
Cy =o + = / F (23) 
yp T fate 
in the Timoshenkos Theory: 
1 
C64 Ve A at (24) 
| eager eee 1 


The Theory of wave propagation based on the method of ,,Internal Constraints“ gives for the group- 
velocity of flexural waves the following expressions?: 


a) in three-dimensional elasticity : 


1 } 
GO, =O ile 2y¢? (25) 
‘ ie Aaa | re) ( 
l L— 2 Co 
1 Equations (25) and (26) can be derived from the frequency equations (19) and (20). In fact writing these 
equations under the form 
xt—y(d)2?+C=0. (27) 
Where C is a constant and differentiating equation (27) with respect to 6 one obtains 
dx dy(6) dx 
Si aeeeN.8 = 
A x i dé 2x y (6) =0 
or 
dx sy dy(0) 


and in view of equation (21): 


| 5 4) | 
Cy dx dé a. 
gar Sears ee Late apeRes oe 5 


from which equations (25) and (26) immediately follow, 
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b) in two-dimensional elasticity : 


1 ) 
ee ile cee, . eC ( (26) 
1+ 4ato|(T—3) 4(1 »(£)| | 
Equation (25) becomes in the case of a cylindrical bar of radius a 
eet fae eal : (28) 
} (eee) TE sees nator (<)| | 
l (29) co} | J 


As equations (25) and (28) show, when 
6=0,—> 0: 


when Jes Coe == (06 


The same applies in the Timoshenkos Theory. 
In Tables 4, 5 and 6 the values of the group velocities c, derived from equations (24), (25) and 
(26) are given for the following values of the parameters y and 6: 


y = 0.0; 0.05; 0.10; 0.15; 0.20; 0.25; 0.30; 0.35; 0.40, 
6 = 0.01; 0.05; 0.10; 0.50; 1.00; 2.00; 3.00; 6.00; 10.00; co. 


In computing Table 4 the value K = 8/9 has been assumed for Timoshenkos Shear Coefficient. 


Table 4. Group velocities for flexural waves according to Timoshenko’s Theory. 


| 

ce 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

0.01 0.125 0.125 0.125 0.125 0.125 0.125 0.125 0.125 0.125 
0.05 0.4888 0.4853 0.4818 0.4785 0.4750 0.4718 0.4684 0.4654 0.4620 
0.10 0.6530 0.6432 0.6338 0.6244 0.6157 0.6071 0.5986 0.5905 0.5826 
0.50 0.6861 0.6686 0.6524 0.6373 0.6231 0.6099 0.5974 0.5857 0.5748 
1.00 0.6728 0.6562 0.6407 0.6263 0.6128 0.6002 0.5883 0.5772 0.5667 
2.00 0.6683 0.6521 0.6370 0.6229 0.6098 0.5973 0.5856 0.5747 0.5642 
3.00 0.6674 0.6512 0.6362 0.6222 0.6091 0.5967 0.5851 0.5742 0.5638 
6.00 0.6669 0.0507 0.6358 0.6218 0.6086 0.5964 0.5848 0.5739 0.5636 
10.00 0.6667 0.6506 0.6357 0.6217 0.6086 0.5963 0.5848 0.5738 0.5635 
co 0.6667 0.6506 0.6356 0.6217 0.6086 0.5963 0.5847 0.5738 0.5634 


Table 5. Group velocities for flexural waves according to Three-Dimensional Theory of Internal Constraints. 


bie 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
0.01 0.125 0.129 0.129 0.131 0.131 O35 0.143 0.156 0.180 
0.05 0.4969 0.4946 0.4943 0.4965 0.5018 0.5107 0.5250 0.5462 0.5780 
0.10 0.6757 0.6673 0.6604 0.6553 0.6518 0.6497 0.6490 0.6481 0.6443 
0.50 0.7301 0.7114 0.6938 0.6771 0.6611 0.6549 0.6312 0.6170 0.6034 
1.00 0.7150 0.6971 0.6805 0.6648 0.6497 0.6363 0.6232 0.6108 0.5991 
2.00 0.7093 0.6920 0.6753 0.6608 0.6467 0.6344 0.6210 0.6091 0.5980 
3.00 0.7081 0.6909 0.6749 0.6600 0.6461 0.63333 0.6205 0.6088 0.5977 
6.00 0.7073 0.6902 0.6743 0.6595 0.6457 0.0326 0.6203 0.6086 0.5977 
10.0 0.7071 0.6902 0.6742 0.6594 0.6456 0.6325 0.6202 0.6086 0.5976 

co 0.7071 0.6901 0.6742 0.659 4, 0.0455 0.6324 0.6202 0.0086 0.5976 


Table 6. Group velocities for flexural waves according to Two-Dimensional Theory of Internal Constraints. 


— Y 
5 a 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

0.01 0.125 0.133 0.141 0.147 0.153 0.167 Oats 0.188 od 
0.05 0.4969 0.5195 0.5442 0.5713 0.6014 0.6347 0.6720 0.7140 0.763 
0.10 0.6757 0.7010 0.7283 0.7577 0.7897 0.8244. 0.8624 0.9109 0.9502 
0.50 0.7301 0.7480 0.7673 0.7883 0.8111 0.8354 0.8641 0.8950 0.9300 
1.00 0.7150 0.7330 0.7524 0.7736 0.7967 0.8223 0.8506 0.8822 0.9177 
2.00 0.7093 0.7275 0.7472 0.7687 0.7922 0.8181 0.8466 0.8784 0.9141 
3.00 0.7081 0.7264 0.7462 0.7678 0.7912 0.8172 0.8458 0.8776 0.9134 
6.00 0.7074 0.7257 0.7455 0.7672 0.7908 0.8166 0.8453 0.8772 0.9130 
10.00 0.7072 0.7256 0.7454 0.7670 0.7906 0.8165 0.8452 0.8771 0.9130 
oa) 0.7071 0.7255 0.7454 0.7670 0.7906 0.8165 0.8452 0.8771 0.9129 
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In fig. 2, which, like fig. 1 was derived from Professor R. M. Davies paper, the ratio c,/c, bet- 
ween group velocity of flexural waves in a cylindrical bar and of the velocity of loneinudivel waves 
of infinite wave-length are given versus the ratio a/L of the radius of the cross-section of the 
bar and the wave-length in the case in which »y = 0,29. 


The curves obtained from the Elementary Theory [equation (22)], from the Rayleighs Theory 
[equation (23)], from the Timoshenkos Theory [equations (24)], and the results abtatuell from the 
»Method of Internal Constraints“ [equation (28)], as well as the curve for the first mode of exten- 
sional waves are compared with those obtained from the Exact Theory by Hudson and Davies. As 
before, the results from Timoshenkos Theory are given as a series of points, the crosses representing 


a) 
Co 
Gil 
40 1 | 
Flexural Waves Flexural Waves 
LV (Elementary Theory) (Rayleigh Theory) 
‘ey | 
| Extensional Waves (First Mode) 
98 | 4 
0,7 = ee Waves 4 ) 
Method of Internal Constraints 
Tez 
06 | Lok pe ex x x { 
yo / Oj QO ae D a ry 0 dt x 3 i. 
x 
O5 ihe 
i 
x 
att : ai sl 
Flexural Waves 
o3t4-(Exaet Theory) Bhs 
02 + t | + iP i 
Values from 
Timoshenko’s Theory 
Of - —x— k=0,9(Davies) + 
—o— k=0,847 (Mind/in) 
t a 
(ee cert a 
0 02 OY 06 08 10 42 U4 46 7,8 2,0 


Fig. 2. Group velocity of flexural waves in cylindrical bars (y = 0,29). 


the values of the ratios c,/cy computed by assuming K = 0,9 (according to R. M. Davies) while the 
squares representing the values of the ratios c,/cy computed by assuming K — 0,847 (according to 
R. D. Mindlin). As fig. 2 shows, for wave-lengths greater than about ten times the radius of the bar 
(a/L < 0,10) all four theories give very closely the same results. For increasing values of the ratio 
a/L the Elementary Theory and Rayleighs Theory give values of c,/cy greater than those given by 
the Exact Theory and when the ratio a/L— o the Elementary Theory gives an infinite value for 
the group velocity while the Rayleighs Theory gives c, = ¢. 

Instead in the Theory of Internal Constraints as well as in the Timoshenkos Theory to large 
values of the ratio a/L correspond values of the group velocities equal to the values of the corre- 
sponding phase velocities. 

The Exact Theory predicts that the group velocity will have a maximum for a = 0,3 L. This 
implies that when a flexural pulse is propagated along a bar, Fouriers components whose wave- 
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lengths are about three times that of the radius of the bar will travel ahead of the components 
having other wave-lengths, and will therefore appear at the head of the pulse ?. 


The same applies for flexural pulses according to the ,,Theory of Internal Constraints. In fact, 
as fig. 2 shows the maximum value of the ratio ¢,/cy given by the approximate theory occurs at 
a/L = 0,3. This good agreement between the results given by the approximate theory presented 
here and the exact theory is perhaps the most interesting feature of fig. 2. 


5. Torsional Waves. Consider a bar with a circular cross-section and let the components of the 
displacement be 


0 = 2 fly =— 2, (29) 
WV AO Wie 
The motion of the bar will be characterized by equations (4), where 
— ply =Ag =O, Ty HM 
In fact, since in cylindrical coordinates u, and uy 
u, =vcos 9+ wsin O = (uy +4) = 0, 
: 2? y? 
Ue =— vsin O + w cos OF Hy at As r 
the condition — wy = /, = w follows. Equations (4) become 
Ce (30) 


From equation (30) it follows for the velocity c, of propagation of torsional waves 


1 

C= V—. 

Q 

6. Compressional Waves. Let the components of the displacement of the bar be 
w= yx 1), 


v= ¥ A(x, 1): (31) 
W == 2% IL,(X, t). | 


The motion of the bar will be characterized by equations. (5) Let 


ram Noite), | 
hy — J, e'v(*—et) F 
lee M eiv(x—et) 


(32) 


where N,L,M are three constants y = 2 7/L and cis the wave-velocity. Substituting equations (32) 
into equations (5), one obtains three homogeneous linear algebraic equations in N, L and M whose 
determinant set equal to zero, yields the following frequency equation: 


/ a+t2, a+t+2 9a+6 5 6 4 4 
Eel barn aes tool +p + page +2045) | 


en ur de ry rey 
oeey  Wboacy Ssigese (33) 
Ha tap tap) tat e=0. 
In equation (33) 
eee ie h 
oar ye o=Ve, ae 

If the bar has a circular cross-section of radius a, then 
Ne (34) 


1 See H. Kolsky, 1. c, 


—— 
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In fact, since in cylindrical coordinates u, and ug 


u,= ves O+wsinO=1 +=, | (35) 
ue =— vsin 0 + wcos O= (—/, + p,) xy 
and since 
ue = 0 
equations (34) follows, and the first of equations (35) can be written 
up =A,T. 

Instead of equations (31) the equations 

u = W(x, t) = 7,(x, t), | 

Up =— U(x, t)=Ayr ist] 


will be used. 
Substituting equations (36) into equations (5) with 


equations (5) reduce to Professors Mindlin and Herrmanns equations! 


o?U 
Oe? 


| 37) 
au aw aw ( 
ME il taa eG) cory cakes ar | 


aw aU 
—8(A+6) U—4015 +6] = 0 


By eliminating U from equations (37) the single equation in W 


CW  Eew a OW a E(l—v) ow (1+ 7)(1—2»r)ea ow 
oa rete Fe lad ora tl 8 dx? or 8a (l+v) dxt 4E ant 


0 (38) 


is obtained. 
The first two terms of equation (38) are the terms given by the Elementary Theory of com- 
pressional waves propagation while the other three terms are the corrections for radial inertia and 
radial shear. The frequency equation corresponding to the waves defined by equation (38) is 
Poe (ae) 07), 1+ 2776? (1 —1?) 
Se Aies>) a2) niet AL) 2) 


Zaet a/E a2 
x= a 5 Co = Oo 5 = i is 
In order to discuss the problem of compressional waves in the case of plane strain, let the compo- 
nents of the displacement be 


4 


x 


=i (39) 


where 


u = »,(x, t) , (40) 


v= vA,(x,t). 


The motion of the bar will be characterized by equations (7). By eliminating /, in equations (7) one 
obtains the single equation 


eu, E eu sl —) Oty E 2du | (3—r) , dtu 
te = 9 a8 @" EB Ot 20 +0" Ox yey ew (41) 
The frequency equation corresponding to equation (41) is 
2 a? 6? 
4 L=p 2 27 67 (3 —) Daas on (42) 
iat ln ita) 


1 See R. D. Mindlin and G. Herrmann, 1. c. 
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7. Conclusions. The results obtained in this paper in the particular case of lateral vibrations of 
bars show an encouraging agreement between the values of c/c) and c,/cy given by the approximate 
theory based on the assumption of internal constraints and the exact theory derived from the use 
of the equations of the Mathematical Theory of Elasticity. 

This approximate theory which will be referred to as the ,,Theory of Internal Constraints“ is in 
the dynamic case completely contained in the constraint equation (1) and in the application of 
Hamiltons Principle. Accordingly the concept of Shear Coefficient is not used. In the general case 
of wave propagation in elastic straight rods this theory unifies a number of separate engineering 
treatments of the problem. Moreover, the same theory can be applied to the study of vibrations of 
curved bars taking into account the effects of shear and of rotatory inertia as has been shown in a 
previous paper. 

The mathematical simplicity of the theory and its degree of accuracy justify its use in dealing 
with engineering problems in vibrations of curved or straight bars for which more exact theories 
cannot be used because of their mathematical complexities. 

The author wants to express his best thanks to Mr. E. C. Zachmanoglou, student in the Depart- 
ment of Aeronautical Engineering, and to Mr. R. V. Milligan, graduate student in the Department 
of Mechanics, at Rensselaer Polytechnic Institute for the valuable help and assistance given to him 
in the preparation of the numerical tables and graphs which are presented in this paper. 
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Anschrift des Verfassers: Prof. Enrico Volterra, Troy, New York (USA) Rensselaer Polytechnic Institute. 1 
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Beitrag zur Theorie der elastischen Ringe mit Beriicksichtigung 
der Woélbbehinderung 


Von A. Bjérklund 


1. Einleitung. Die Theorie der stark gekriimmten offenen oder geschlossenen Kreisringe mit 
beliebigem Querschnitt und allgemein formulierten Belastungsannahmen ist schon friiher von 
F. K. G. Odqvist } behandelt worden. Im folgenden wird diese Theorie soweit erweitert, daB auch 
auf die Kinwirkung der Querschnittsyerwélbung bei Torsion Riicksicht genommen werden kann. 
Beim Aufstellen der Verformungshypothesen und der Grundgleichungen ist vorausgesetzt, daB, da 
eine exakte Theorie fiir den gekriimmten Balken fehlt, die fiir die Verdrehung des geraden Balkens 
geltende Theorie? *, mit gewissen Anderungen auch in diesem Falle angewandt werden kann. 
Die erhaltenen Ergebnisse sind also — besonders bei starker Kriimmung — nur als eine 
erste Naherung anzusehen. Eine Durchrechnung von praktischen Beispielen zeigt jedoch, daB man 
durch Beriicksichtigung der Querschnittsverwélbung oft erhebliche Abweichungen bei den berech- 
neten Verformungen und Spannungen gegeniiber den Werten ohne Beriicksichtigung der Verwél- 
bung! erhalt. Dies gilt hauptsachlich fiir offene Profile. Bei véllig geschlossenen Querschnitten 
dirfte die Wélbwirkung, ebenso wie beim geraden Balken, im allgemeinen wohl keine praktische 
Bedeutung haben. 


Um einen Vergleich mit der Odqvistschen Lésung zu erleichtern, schlieBt sich die folgende Dar- 
stellung dieser soweit wie méglich an. 


—— Ty 


Abb. 1. Allgemeines Koordinatensystem fiir ein Kreisringstiick. Abb. 3. Momentenzeiger am Bogenelement r, dy. 


2. Gleichgewichtsbedingungen. Es werden zwei bewegliche Koordinatensysteme (Abb. 1) ein- 


gefiihrt, das eine mit x = r dp als Bogenlange der Schwerpunktslinie, das andere mit x, = Tr; dy als 
Bogenlange einer Linie durch die Schubmittelpunkte der Querschnitte. Die SchnittgréBen M, , 


M, (Biegemomente), M, (Torsionsmoment), Q, T (Schubkrafte in axialer bzw. radialer Richtung) 
und N (Normalkraft), sowie die auf der Langeneinheit angreifenden, streckenweise stetig verteilten 


1 F. K.G. Odqvist, Ing.-Arch. 22 (1954), S. 98. . 
2 7. B. F. W. Bornscheuer, Stahlbau, 21 (1952), S. 1., S. PIONS Ap) (1953), Sour 
3 N. J. Hoff, Journ. Royal Aeronaut. Soc., 47 (1943), S. 64 (Appendix). 
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auBeren Krafte und Momente p, q, t, m,, m, sind auf das Bogenelement r, dy zuriickgefiihrt gedacht 
(Abb. 2 und 3). 

Wenn die AnfangsgréBen fiir p = 0 der betreffenden SchnittgréBen mit Ny, Ty, Qo, Mi», Meo 
und M,,, bezeichnet werden, ergeben die Gleichgewichtsbedingungen unter der iiblichen Voraussetzung 
der vernachlassigbaren Verformungen: 


N= N,cosp— Ty sing + [p r, sin (p — A) —t r, cos (p—A)] dd, (1) 
0 

T = T, cosy + Ny sing — | [p r, cos (p— A) + tr, sin (p—A)] dd, (2) 
0 

@=Q—fand, (3) 


0 


” ” 
M, = (M,, + Qo 7) cos p— M,, sin p— Qo 4 + far da — [ (m,r, + q 1?) cos(p—A) dA, (A) 
0 0 

D 
My = (Mi +Q 1) sing + Myo cos p— | (m,r, + qr?) sin (p—A) dd, (5) 

f 
M, = M,, + Tyr, sng + Nor (1— cosy) — f mr, dA 

0 


— f {p rsin(p—A) + #73 [1— cos (p—A)} a. (6) 


In diesen Integralen sind p, q, t, m, und m, als Funktionen der Integrationsveranderlichen / 
aufzufassen. 


3. Verformungshypothese, Spannungsverteilung. Die Entwickelung der Lésung geschieht unter 
folgenden Annahmen: 

1) Die Verdrehung jeder Querschnittsflache erfolgt, ebenso wie beim geraden Balken, um den 
Schubmittelpunkt, d. h. um die Koordinatenachse x,, die also die Hauptdrehachse bildet. 

2) Die Ringdehnung ¢ in der x,-Richtung setzt sich aus zwei Anteilen ¢, und ¢, zusammen. 
Dabei ist ¢, die Dehnung, die mit der Annahme des Ebenbleibens der Querschnitte nach Bernoulli, 
€, diejenige, die durch die Wélbbehinderung bei Torsion erhalten wird. 

3) Die Dehnung ¢, und die dadurch hervorgerufenen Spannungen werden im Prinzip der Theorie 
fiir den geraden Balken ? entnommen. Wie aus dem folgenden ersichtlich, werden die in Frage 
kommenden Gleichungen ? nur insoweit erginzt, daB, wenn der Verdrehungswinkel beim geraden 
Balken mit @ bezeichnet wird, sein Differential dQ durch dj}—vwy dp ersetzt wird. Die ¢, ent- 
prechenden Spannungen geben somit schon ein geschlossenes Gleichgewichtssystem. 


Abb. 4. Verschiebungsvektor. Abb. 5. Winkelanderungen. 


Mit Riicksicht auf die Annahme 1) werden simtliche VerschiebungsgréBen, ebenso wie die 
duBeren Belastungen, auf die Hauptdrehachse bezogen. Die positiven Richtungen gehen aus Abb. 4 
und 5 hervor. Entsprechend der Theorie von Odgqvist erhalt man 

~ 1 3 : vy’ 
reer Bah esa a) Sih , (7) 
wobei 


1 =rT+%>, bp) = Gp = Uys %=f—Cy (8) 
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und 
, du Votes dv 
dp’ dp 
ist. 


Nach Bornscheuer? gilt fiir den verdrehten geraden Balken mit dem Verdrehungswinkel 0 
Eo = — WwW dx? P) (9) 


wenn W die fiir den jeweiligen Punkt des Querschnitts geltende ,, Hauptverwélbung* ist. 


Fiir ein Element r, dp eines steifen, unbelasteten Ringes(Abb. 5) gilt jedoch, da die Vektoren#, 


y, x in diesem Fall an beiden Enden einander gleich sein miissen und die Winkel #, y, x und ihre 
Ableitungen klein sind, 


dij—ypdp=0, (10a) 
dy+8dp=0, (10b) 
d= 04 (10c) 


Dies wird deutlich dureh Projektion der Winkel auf die Richtungen der Tangente und der Nor- 
malen des Elementes sowie auf die Normale der Ebene des Ringes. 

Das Element erhalt somit gema8 (10a), ohne verformt zu werden, durch y einen Beitrag d? zum 
Verdrehungswinkel. Die durch das Torsionsmoment erhaltene Vergréferung des Verdrehungs- 
winkel betragt dann d} — y dy. Wird also dO = di — y dp gesetat, so gibt (9) mit dx, = r, dp 

hes = (11) 


ry 


wenn Ableitungen der Winkel nach g mit Strichen bezeichnet werden. 


Die Gleichungen (10b) und (10c) entsprechen den Gleichungen (28) bzw. (29), die fiir das be- 
lastete Element gelten. 


GemaB dem Hookeschen Gesetz bekommt man sodann fiir die den Dehnungen ¢, und ¢, ent- 
sprechenden Spannungsverteilungen 
oO, =-E &\ 5 On Le, (12a, b) 
oder insgesamt 


o=E (e+ e). (12c) 


4, Die Grundgleichungen. Die Aquivalenz auSerer und innerer Krafte im Querschnitt p ergibt 
gemaB Annahme 3) 


Ne ondh (13) 
NM oi, dE, (14) 
M, = [fo dF. (15) 
Werden die Bezeichnungen 
(ich saay (16) 
Cy dF (17) 
1 + n/r j 
gc oe 18 
| 1+ y/r J, (18) 
CdF 1 
| ee ve 


eingefiihrt, so erhalt man 


ae =+(F +2). (20) 


eee (21) 


SB SBE gp (22) 
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Wird auBerdem 


l es 

iD) a 

a (i= Dy (23) 
tenet | 


gesetzt, so wird durch Einsetzen von (7), (8) und (12a) in (13) bis (15) und mit (16) bis (22) 


Ne (ul +0) (FP +S) +O +9) DrtQ(FAt+D]4+—Ur+ mr + J), 2) 
ee ——(w +) [Dr+q(Fe+J]—@+y) Tre+2GqDr+Q(Fr+J)] 
—2 [Dry tooJr+lom(Fr + J). (25) 
ee =—(u' +0) [Jr+m(Fe+J]—(O+ 4’) (Dry +o J r+ lym (Fr + J)] 
<= [Jr + 2m Jr +g (Fr + J). (26) 
Werden diese Gleichangen nach den VerformungegréBen aufgelést, so ergibt sich mit (23) 
MEG Br MOE Gr M, am 
Dy ay es ee ee 
JEM Ns Meta Mi 29 
Fur den geraden Balken gilt 
M,=¢K. ECe. (30) 


wobei C der Wélbwiderstand bezogen auf den Schubmittelpunkt und GK die Verdrehungssteifigkeit 
des Balkens bedeutet. 
Mit d(O=d)—ydy und dx, =r, dp erhalt man 
MeO Re eee 


3 
Ty ry 


(31) 
Tn Gleichung (31) mu8 bei dieser Annaherung C und GK nach den fiir den geraden Balken geltenden 


Formeln berechnet werden. 
Zusammen mit (27) bis (29) und (1) bis (6) erhalt man damit die erforderlichen Gleichungen fiir 
die Berechnung von u, v, ? und yw als Funktionen von g. Die Spannungsverteilung o ist somit 


durch die Gleichungen (7), (11) und (12a) bis (12c) bestimmt. 


5. Integration der Grundgleichungen. Im folgenden wird der Querschnitt als iiberall gleich vor- 
ausgesetzt. 


Folgende Bezeichnungen werden eingefiihrt: 


1 iy r r r GK , 
pepo gy Pe ez, Pe py Bes GR: Bite ht 0, pene ee 


Als Anfangsbedingungen fiir p = 0 wird fiir die Verformungsgréen angenommen 


u(0) = up, (0) =r, v0) =r, w0)= wm, w'(0) = (0) =y, 
00) =O) , 90) =, 070) = 0; « (33) 
Aus (29), (32) und (33) erhalt man nach Integration 


QV —Pp P 


ene. n | M, a—(a% + By n) i M, di —n (0 — iba + Bi ns | Nan. (34) 


0 0 
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p A, 2 A, 
vate thin | | My dt day — (0% + Ban) [aay 
0 0 


0 0 


We A, 
no - Bi So Bano) | [raw 
0 6 


Fiir die in (34) und (35) und im folgenden auftretenden Integrale gilt mit (1) bis (6) 


0 


Pp Ay 
{ Ndd di, = N, (1— cos g) — Ty (py — sin yy) 
06 


p 
+ fpr [p — A— sin (p— A)] da—frr, [1— cos (p—A)] dd, 
0 


0 

gp A 

ily 'Ndidh d/, = MN ~—sing) — T(z g—1 + cosg) 
0 


| 
? 


ae (p —2)?—14 cos (p—A)| da — | try{p—A— sin (p—A)] 42; 


0 
? 

JM, dd = (M,, + Qo 7) (1— cos y) + M,, sin p 

0 


ae / (mr; + q7}) [1— cos (p—A)] da, 


0 


1 


ph 
| [Mp dA dh, = (My + Q 71) (p — sin y) + M,, (1 — 008 @) 
00 


ain (m,7r, + ri) [p—A— sin (p—A)] a, 
0 
gp A, A, l 
| | | M, dA dA, dy = (Muy + Qn) (s 92— 1 + 08 2) + Myo (p— sin 9) 
OO 80 


” 
=| (m, ntar) |e @—ae—] + cos (p—A)| dd, 
6 


jm. dA = M.y 9 + Ty r (1— cos p) + Nor, (p— sin ¢) 


’ ” 
— f m.r,(p—) di— | pri [1 —cos(p—A)] da— | tri[p—A—sin (p—A)] dad, 
0 6 


0 
aa 1 1 
| M, di dd, = 3 Mov? + Tyr, (py — sing) + N (gel + ¢08 9} 
@ 9 
gy — A)? dA— [ prifp—’—sin(o—A)] dA 


0 


— 


pais 
0 


|i rily a —1 + c05 (p—A)| d) , 


” 9 s 
[Nda=N, sin y — Ty (1 — cos 9) + fpr, [1 — cos (p — A)] dA— ane sin (p — A) di, 
0 0 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 
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Ay As 


Einsetzen von (35) in (27) gibt 


Ve QP 
usm [edi +(x% +feti—2 Aion + Bora) | Na 
6 d 


P 
[| | M. dh da, di, 
0 0 0 


1 1 1 ; 
= oet Tyr, (q 92-1 + cos) + Non (G oo + sing) 


” P 
Uh of 1 
ae m, 7, (p— ap aa— | Pils @—A—1 + cos (p—A)| dj 


0 


— [vi L@—as—@—a 4 sin (pA) aa. (44) 


Uta fan) [ Meat (02 — Ait Ba) [ Ma dh , (45) | 


wobei das erste Integral aus 


P 


Q As Pp 
—(»+Bsn) | [ [ Mededn, ddy—n (Pate Pare) | 


erhalten wird. 


P 


Aves 
[et=netynet bin | | | Median a, 
0 O70:246) 


Ay Ay 


Ay 
| Ndidd, di, (46) 
0 


ORO n0 0 


Aus (28), (31) und (32) erhalt man die Differentialgleichung 
OW + (1—«) 0’ —ed =e f,M,—ep,M,—e6dN—ey,M, —f$.M,+4+ 8, M7+6N"”, (47) 


welche die Lésung 


0 = C, Gin Yep + C, Cof fey + C, sing + C, cos y 


ae 


ne 


tlie | Gin Ve (p—A) es sin (p — A) vi 
te ‘lie | M, Gin je (yp —A) ua—|M, sin (p — A) vi 


eae et ace p 
— el fe Neieto—a da— | Nein(p —2) wi 


—ey; i | M, Gin Ye (p — A) dA — | M; sin (p — A) vi (48) 


tale M! Gin Ye (p — A) ua | me: sin (p— A) vi 
~alie | My Gin ye (p—A) ua | Me; sin (p — A) vi 
Bite he |?” Gin ye (p—A) u— fw sin (p — A) al 


0 
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hat. Werden die Anfangsbedingungen (0) = 0), 0” (Oyen OVSs0, (0) = Wp erfillt, so 
werden aus (28), (31) und (48) und mit Riicksicht auf 


j FO Gin fe (p —A) dA =— FP” Gin Pew + Ve { FO) Gof fe (p—A) aa, (49a) 
j FO Gof Ve (p— A) dd = FU-) — FY oj peg + ye i FU-) Gin fe(p—A) da, — (49b) 
[RO sin (py — A) dA =— FG 1) sing + [ae ) cos (p-—A) dd, (49c) 
0 

| FG) cos (p — A) dd = FU) — FG cos y — il "FG-) sin (p— A) di, (49d) 
0 0 


ae Lésungen erhalten: 
Ve (Yo — 99) Gin Pep + (Bo + I —By Meo + Bx Myo — 5 M) Cof Veg 
+ (€Y% + Bo) sin p + (€ By — 99 + B; Mio — fy pmb eng 


= aoe | 
oe tere 1) |p, fv sin (p — A) aia sin (y — A) Ui+0 | Nsin( (p— A) da i 


a a li M, Goj Ve (p — A) aa— | M, cos (p-— A) wal, 


O = ale mH) Cof Veo t+ Ve(e +%—p, Mz) + By Myo— 0 M) Gin Vey 
+ (€ Wp + Bo) cos yp — (€ Fy) — (toes Myo +0) sing 
+ (e+ 1) lf [ M. cos(p — A) dA — pa. [ M, cos (p — A) ) d+ d | Neos (p 


=e 
—ey,|Je [min fe (p—A) df M, sin (p—) aa] : 


a {22 —09 Gin Pep + € (By + HG — Br Mey + By Myy— 0 No) Gof Ye o | 


= (eat 90) ait 9 — (CO — OF FB Mag Ps Mat 6 N;) 008 9 
+ (e+1)|B, M.— P.M, + 0N—p, | Mesn(p A) a 


1 
e+ 1 


U 
{ 
+B. M, sin (p —2) d—9 | Nsin (yp —1) aa| | 
elem Coj Ve (p—A) a+ JM, cos (p— A) in| ; | 


{(vo— 96) ) of Vep— z= Wo + % “8. Mig + fo Mpp — 6 Ny) Gin Ve p 
+ (€Y% + 9) cos p— (€ Fy — Ho prio abies Myo + 6 No) sing 
Sete +14, fm. cos (py — A) ) da —B, | M, cos (p— A) ta+5 ['N cos(p—2) aa| 


+e ae Gin ye (P—A) a — | M, sin (py — 4) a). 


v=o 


? ae eet 1 
— (€ Yo + 9) sin p — (¢ Dy — Do + B1 Mio — Bo Moo + 5 No) cosy 


+ (e+ 1)|fM.— PM, + 6N—p, | M, sin (y — Ah) da (50b) 


ae {V= (ro — 96) Gin Vey — (Oo + 86 — By Meo + Bx Moo — 3.) Gof Vey 


+ Py | My sin (9 — A) dh— 6 i N sin (p — A) al 


+eyy lf a, Coj Ve (p —A) uf M, cos (vy — A) in| 


29 
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Unter Ausnutzung von (1) bis (6) kénnen die in (48a), (48b), (48c), (50a) und (50b) auftretenden 
Integrale nach Vereinfachung wie folgt geschrieben werden: 
yp 
| N sin(y — A) dA => Nyy sing = T, (sin p — 7 cos 7) 


0 yp 


6 0 
us 
| N cos (vy — A) eae (sin p + 9 cos 9) — > Ty p sin ~ 
? ’ 
: 1 5 
| p ry (p—A) sin (p — A) dA—-> [ tr, [sin (p — A) + (p— A) cos (p—A)] dA, (52) 
) 6 


~ 
[ Ms sin (p —)) dh => (Mio + Qo 71) (sin y — @ cos @) +5 Myo p sin yp 
6 : p . 
<x | (mn + ard [sin(o—2)— (=A) 008 (P—1)] dd, (53) 
j 
ve 
| M, cos (vy — A) dA = = (Mig + Qo71) psing +> a (sin p -+ @ cos 7) 
) 


: | 
—7| (sri, 5, ct gr) Glee Ay ei (ee) 


Pp 
| M, sin (p — A) dA = M,, (1 — cos ¢) Le Ty r, (sin p — —7 cos @) 


0 Pp 


-Ny Wr, Ji cosp— a9 sin g}] — | m, r, [1 — cos (p — A)] da 

et ae (55) 
ar [en [sin (gp — A) — (p — A) cos (p — A)] dA 

a 
— | t rj [1 — cos (p — A) : (vy — A) sin (py —A)] dd, 

0 


8 


| M. cos (p g— A) di = M,, sing +> = Tyr psing +N, r, (sin y — @— cos ¢) | 
y ? y | 


— | m, r, sin (~ —) da— + | pr; (es A) sin (p — A) dA 


| (56) 
0 9 
; 
1 ae oe 
—y [| tri bin @—2)—— (PA) c08 (P— A Mh, 
0 
ioe : 1 ‘ 1 ; 

| sin (p — A) dd = > (Muy + Qa) Sin p — ~} Myy (sin y — 7 cos @) 
0 


— Qyr, (1— cose) + | qrj [l—cos(p—A)] dA | (57) 
0 : 
2 


(m, 7, + ¢ 71) (p—A) sin (p—A) di, 


tole 


i=) 
— 


a [px [sin (p= 1) (p— 1) cos(p—A)] u—z | tr, (p—A) sin (p—A) da, (51) 
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? 

4 1 : 
| M, cos (p — A) dA = a (Miy + Qo 75) (sin p + ¢ cos g) — poy sing 
' Co Pp 

: 2 Ub ve ‘ 9 
—Qorisin p+ qrisin (p—A) da — 5. | (mr, + 474) [simp —2) + (p—A) cos (p—A)] da. (58) 
0 0 
Pp 


Je 


| M,Sin Ve (p — A) dk = —— (Mig + Qo 11) (Cof Ve p— cos ?) 


5 ite 
l ; ms ramet ~_/[7- 
Saar ine (Gin Ve vp —Ye sing) — 7 Q, 7, (Gof Ve ~—1) 
ae o | qr? |Coj Ve (p—A)— 1] aa (59) 
0 


Pp 


ar | (m7, + q 73) [Gof |e (p— A) — cos (p—A)] dd, 


| M, Goj Ve (p — A) dA = | (Mi + Qo 7) (Ve Gin Ve w + sing) 


j 
1 a 1 as 
sq Moo (Cof foe are 28) r, GinJe p 
f - (60) 
a = | qr? Sin Ve (p—A) da 
é 5 ; 
taf - = 
ave / (m,r,+ q 72) [Ve Gin Ve (p — A) + sin (p—A)] dd. 
0 
SchlieBlich ist ‘ 
w=wt+ry [v d , (61) 
woraus mit (50a) : 
ry 1 , ; - i , ps 
BOR ge aa fe (Yo— Fo) Sin Yep—— (+08, Mz +B, Myo— 5 No) (Gof Ve py — 1) 
+ (€ Yo + Bo) sin p — (€ Jy — Bg + Bi Meo — Bo Myo + 6 No) (1 — cos g) 
i ” 
+ (e+ 1) a. f M, sin (p — A) dA— f, | M, sin (p—A) dd (62) 
0 


¢ 


+6 [Nn sin (p— A) a| 


—éEy, [2 fa, [1 — Coj ie @—A)] dA + iM, [1 — cos (py — A)] al} 


folgt. 
Die in (62) auftretenden Integrale kénnen mit Hilfe von (51), (53), (55), (58), (60) und der Forme] 


Y 

iE M, dA = (Mig + Qo 7;) sin p — My, (1 — cos ~) —Q) 1, 

: ? ’ 

== | qri(p—A) ds —| (m,r, + q ri) sin (p—A) dd (63) 
0 0 
berechnet werden. 

Aus (34) bis (63) konnen samtliche Verformungen als Funktionen von y und von den 14 Anfangs- 
werten — Ny, Ty, Qo. Mig, Myo» Mz und (33) — gefunden werden. 

Die Lésung schlieBt das Auftreten von konzentrierten Lasten und Momenten innerhalb des 
betrachteten Gebietes aus. Falls solche vorkommen, mu8 der Ring in Intervalle zwischen den Un- 
stetigkeitspunkten der auBeren Belastung aufgeteilt werden. 

29% 
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Bei statisch unbestimmten Problemen, bei welchen die Anfangswerte voneinander abhangig 
sind, miissen ebensoviele Randbedingungen wie statisch Unbestimmte beriicksichtigt werden. 
Von diesen diirften nur diejenigen, die zur Bestimmung von 0, und &y erforderlich sind, einer Er- 
klarung benétigen. 

Die Querschnittsverwélbung ist beim geraden Balken fiir jeden Punkt der Querschnittsflache 


proportional zur GréBe O’ des betrachteten Schnittes. Dieser Verwélbung entsprechen die | 


Spannungen o, nach (12b). Bei jeder Ubergangsstelle zwischen zwei Intervallen mu also mit 


dO = di— y dp 


| Ew vr / 
Tee =) = 


sein, falls es sich um einen kontinuierlichen Ring handelt und die Indizes 1 bzw. 2 die zusammen - 

fallenden Endpunkte der Intervalle bezeichnen. Infolge der Kontinuitatsbedingung y, = y, ver- 

einfacht sich (64) zu =8 (64a) 
i= Wy c 


Ew vr / 
Cia) (65) 


Fehlen auferdem konzentrierte angreifende Biegemomente und Normalkrafte, so wird (0;); = (01). 
und y, = y, , wodurch (65) sich vereinfacht zu 


(¥—y), = (8 —Y¥)e» (64) | 


eae (65a) | 


Fir ein freies, unbelastetes Balkenende gilt 


Im folgenden werden zwei der Falle, die von Odqvist behandelt sind, durchgerechnet. 


6. Der geschlossene Ring mit konzentrierten umstiilpenden Momenten. Der Ring sei von n gleich 
groBen konzentrierten Momenten M mit gleichem Zentriwinkelabstand 2 y, = 2 2/n (Abb. 6) be- 
lastet. Aus Gleichgewichtsgriinden ist M 

M,o= (67) 


~~ 2sin Pi 


sowie JVy— 1, — 0, — i, — 0. | 
Zur Bestimmung der tibrigen 9 Anfangswerte stehen folgende Randbedingungen zur Verfiigung: 
1) 30 0s 2) (Gy) = 0, 3) Tay, 
HG) — 0a 9) (oi) — 0, Oe Yor = 
)pP)=9, 8) Mp)=H8(—H). 9a) WG) =8(—G), 9b) O'(G)) = 8'(— H) - 
Die Bedingungen 1) und 2) in (34) eingesetzt 
ergeben 


a fy 


r sin 
Mz mary ae - a Moo s (68) 
Die Momente M,, M, und M, konnen dann aus 
(4) bis (6) berechnet werden. 

Die Bedingungen 1), 3) und 4) in (45) ein- 
gesetzt ergeben uv). Die Bedingungen 6), 9a) 
und 9b) in (48a) und (48c) eingesetzt ergeben 
&) =0, was auch direkt ersichtlich ist, wenn man 

Inher beachtet, da} p = 0 ein Symmetrieschnitt ist. Der 
Wert 0) = 0 sowie die Bedingungen 6), 7) und 8) 
in (48b) und (50a) eingesetzt ergeben das folgende System zur Lésung von #) und % : 


Je (8 + Fo — Bi Meo + Bo Moo) Gin Ven: — (€ 39— Bo + By Mio — Bo Myo) sing, | 
[ , 1 p 
-(e +1) (Pr M,z sin Gy are Ds Myo (sin gy + GY cos a) 


| 


, a mens 1 : 
+ 07 [15 Myo (Gin Ve — Ve sin g,) + 1 Myo 6in g.—p, 008 9,)| = 0, 
J Z ee (69) 
Ve (8p + 96 — Bi Meo + Bo Myo) Sin Ve p — (¢ By— Oy + fi Mz) — By Moo) sin py 


—~ 


: il ; 
Bert) Ai My sin gy — oF Bs Myo (sin py + py cos ) 


= 


1 le = Sars 1 3 
Ev | TE +) Moo (Gin Ve Py —ye sin () + zy Moo (sin P1 — Y Cos 0) —0. 


HY — = 0 ; (66) | 
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Die Bedingung 5) in (62) eingesetzt ergibt dann wy. Damit ist der Verformungszustand vollstaindig 
bestimmt, und man erhalt ¢ aus (12c). 


7. Der geschlossene Ring mit gleichmafig verteilter Last quer zur Ringebene und mit aquidistanten 
Stiitzen. Die Belastung wird als konstant angenommen (Abb. 7). Uber den Stiitzpunkten tritt ein 
statisch unbestimmtes Biegemoment M, auf, dagegen verschwindet dort aus Symmetriegriinden das 
Torsionsmoment. Eine Gleichgewichtsbetrach- 


tung ergibt My 


No= 1,=9,=M,5=0, Z Tt 


sowie 


Mio= 47; (sing, —¢, cosp,)— My sing, , 


| (70) 


M,o= 471 (Pi8ing:— 1 + eosqy)—M, cosg, ,| 297 Gr) 
Zur Bestimmung von M, sowie der iibrigen 
8 unbekannten Anfangswerte gelten die gleichen 
Randbedingungen wie im vorigen Beispiel. Abb. 7. 
Die Bedingungen 1) und 2) in (34) eingesetzt ergeben ; 
M, = qr; (1— 9, ctg 9) » (71) 
wodurch mit (70) 
opt MB 
Mane yap (ae = 1) (70a) 
erhalten wird. 
Aus (4) bis (6) folgt . 7) 
a sing __ 
to ie ers | 
Be a A Ee , (72) 
My = 97 (Ho 1) | 
Vien J 


Die tibrigen Randbedingungen werden wie im vorigen Beispiel behandelt. Aus (48a) bis (48c), 
(50a) und (50b) erhalt man 


b= {( + + By Myx) Gof 2p + (e Oy — 9% — a Myc) cos 


—(e+ 1) fala Mo sin p + an(re sin p + cos p— 1) 
_ (73a) 
= 1 

ey; Lar (Moa + q 73) (Gof Ye p— cos y) +— qr; (Gof Ye p—1) 


ar = (Moo + 917i) ¢ sin p— qr (1— cos ®)| . 
oe 4 We (8p + 0 + By Myo) Gin Ve p— (€ Do — 06 — Bo Myo) sin p 


— Ee [M,o (sin p + —@ cos p) — qr; (sing — ¢ cos @) | 


13b 
ie (73b) 


Sgtties ror DUR tie = ae 
—eé, ee ered) (GinJep—Je sing) + yz (Gin fep—Je sin 7) 


1 9 : 
<7 (Moo + 471) (in p—o cos) 


eee! 
o eon 


‘ (8 + BY + By Myc) Gof Ve p— (e 9) — 96 — By Myo) 008 | 
1 j . 
+ (e+ 1) Pa|Mno (zp sing—cosg) +z arp sing} 


z Es (73c) 
—ey; Lepr Moot an) (Gof e p— cos g) + 47? (Gof Ye y— cos) | 


1 9 . 
— =z (Moo + 971) p sin || , 
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v= spa { ye ot 80 + Ba Meo) Gin Pe 9 — (€ Bs — 8) — By Moa) sin 9 
£518, [Myo (sin g + 08 9) — 9 7i (sing —¢ 008 @))] 
ent tye + an (Gn fey — fining) +42 (Gin Ven—¥er) | 
sole + 471) (sin p— ¢ cos y) — 97; (P— sin || 
y=: ay (8) + OY + By Myo) Coj Ye p— (e 9p — 9% — Bo My) cos 
(ea) fo| Moo (> 9 sing — cos) +34 ry sin 9 
(74b) 


1 3 = 1 9 = 
| en sa (Moo + 971) (Gof Vey — c08 y) +975 (Gof Yep — J) 
1 9 : 
+5 (Myo + 972) sin p— 4 (1— e089) | 


Um die Abweichungen der Ergebnisse der oben entwickelten Methode von denjenigen nach 
Odqvist zu verdeutlichen, wird der obige Fall mit folgenden Annahmen zahlenmaBig durchgerechnet : 


Der Ring sei auf 6 Stiitzen gelagert, d. h. 2 mp, = 7/3. Das Querschnittsprofil sei doppelsym- 
metrisch und entspreche einem I 16, r=r,=100cm. Nach Bornscheuer! gilt fiir dieses Profil 
K=7,ll em’, C = 3100 cm®, ®yax = 27,9 cm?. Zur Vereinfachung wird angenommen, dab 7 
gegen r bei Berechnung von J und J vernachlassigt werden kann, d. h. daB die Werte fiir den geraden 
Balken Geltung haben. GemaB (17) bis (19) und (23) gilt dann 


DEO CN) Si On ee ee ODO ea 
und gemaB (32) 
1 
Bi epee 05 LO 695 + (kgem)7+, 4, = 36, 568 = (kgem)=4*5) c= 03206" 


Aus (70a) erhalt man M,, = 0,047 198 qr’, und aus (73b), (74a) und den Bedingungen 6) bis 9) 
qr ” 2 
Jy = 0,011 565 T, dg =— 0,40 29747 


Aus (12a) erhalt man mit (27) bis (29) und (32) 


My 
ip a= i) 


was in diesem Fall (r= r,) identisch mit der Biegespannung ist, die nach Odqvist berechnet werden kann. 


07 = 


C p E, (75) 


Die Spannungsverteilung o fiir p = 0 bis 30° ist aus Abb. 8 ersichtlich, wo die o, entsprechende 
Spannung mit gestrichelten Linien dargestellt ist. Die aus der Wélbbehinderung herstammenden 
Spannungen o, sind offensichtlich erheblich gréSer als die reinen Biegespannungen nach (75). 
Gleichzeitig werden auch die Verformungen in hohem MaBe beeinfluBbt. 


Die Gleichungen (53) und (54) der Odqvistschen Theorie ergeben fiir diesen Fall mit den Rand- 
bedingungen wy, = y(,) = 0 


een a as Payne 2 
= qr? Meg SIN P| — Y — Yj ctg py) ; (76) 
B= 0, cos p— qr? (y, + | Pi i _— 
0 cos p— a7? (M1 + Pr) ring, ? Sin p+ cosp—1}, (77) 
und 
y =— DB sin yp— qr? (y, + fy) bee (p cosy + sin y)— sin elm M. (78) 


1 IF’, W. Bornscheuer, Stahlbau, 22 (1953), S. 34. 
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Die Ergebnisse fiir 0 und p berechnet aus (77) und (78), sind in Abb. 9 mit gestrichelten Linien 
zusammen mit den Werten aus (73) und (74) eingetragen. 


Berichtigung bei der Korrektur. VeranlaBt durch eine briefliche Mitteilung von Herrn 
D. Riidiger an Herrn F. Odquist méchten wir hervorheben, da® die Gleichung (7) oben einer 
Korrektur bedarf, die in gewissen Fallen von Bedeutung sein kann. Es ist das Glied wu/r bei 
der Verdrehung der Bogentangente zu beriicksichtigen, was dazu fiihrt, da® im letzten Gliede 
der grofen Klammer in (7) sowie in den folgenden Gleichungen (24) bis (29) v’ durch v” — u’ 
zu ersetzen ist. Man sieht dann, daf die Kombination von (27) und (29) im Specialfall D — 0 
zu der bekannten Resal-Winklerschen Gleichung fiir v fiihrt. Die Gleichungen fiir 0, y und w 
bleiben unberiihrt. Hingegen sind die Gleichungen (34), (35) und (45) durch folgende zu ersetzen: 


4 7) 
v! =v, cos p— vy sing + (B.r + Bolo) { M, cos (p—A) dd — (Px by + Byr) { Mz cos (p — A) dd 
0 iv) 


+ (Bi 11 S0 + B25 — 2 Bil M9 — Bar no) [ Nos (vy — A) dd, (34) 
() 


f +e f p ; 
vy =v, sing + v, cosy + (B,r + Bolo) { M, sin (y — A) dA — (6,6, + Byr) { Mz sin (y — 2) di) 
( ' 


) ) 


+ (Biri So + Boo — 2 Bi lo 0 — Bz r yo) f N sin (p— A) dd, (35) 
0 


y P ‘ r q 
Uu = Uy + v, cos p—tysing + (6, r+ Bylo) [ M, cos (g — A) dA — (By Cy + Byr) | M, cos (wy — A) dd 
U0 
a 4 r) he y 4 
= (Bx a.Cp > Be. o5 — 2 Pr boe— Pst No) | N cos (py —A) dA —B, r, { M, dd 
0 0 
Y y 


é (« : + Bs r) | M, di + (« = TPs ott pa a “) | Nd). (45) 
( 0) 


) 


| 


Die durchgerechneten Beispiele und zpeziell die Gleichungen (67) bis (78) bleiben auch un- 
veradndert. 

In der oben zitierten Arbeit von F’, Odqvist' sind dieselben Korrekturen zu machen. Bei dem 
dort behandelten Beispiel eines Ringes mit radialen Kraften sind die Gleichungen (67) und (68) 
entsprechend zu korrigieren. 


(Eingegangen am 13. Mai 1955.) 


Anschrift des Verfassers: Civilingenjér A. Bjorklund, Stockholm So (Schweden), Tjarhovsgatan 22. 


1 KO: Odquist, Ing.-Arch. 22 (1954) S. 98. 
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Die Berechnung von Tragfliigelprofilen aus der Druckverteilung 


Von R. Eppler 


1. Einleitung und Ubersicht. Die Tragfliigelprofile mit vorgegebener Druckverteilung wurden 
schon seit langer Zeit fiir vorteilhaft erachtet!»2. Erst als jedoch die Ergebnisse der Grenzschicht- 
theorie und ihre experimentelle Bestatigung sowie die Annaherung der Fluggeschwindigkeiten an 
die Schallgrenze genaue und praktisch wichtige Anforderungen an die Druckverteilungen stellten, 
wurden die Berechnungsverfahren so weit durchgebildet, daB sie praktsch zu vielen Beispielen ver- 
wertet werden konnten. 

Die bekannteste und am meisten verwandte Methode wurde von Th. Theodorsen entwickelt. 
Sie geht aus von dem bekannten ,. Theodorsen-Verfahren‘* zur konformen Abbildung beliebiger 
Fliigelprofile auf den Kreis, das in bekannter Weise erst das Profil auf eine kreisahnliche Figur ab- 
bildet und diese durch einen IterationsprozeB vollends auf den Kreis. Mit der mittelbaren Ab- 
bildungsfunktion laBt sich die Geschwindigkeitsverteilung am Profil angeben. Die von Theodorsen 
durchgefiibrte Berechnung von Profilen mit vorgegebener Geschwindigkeitsverteilung? stellt genau 
die Umkehrung des ersten Verfahrens dar. Seine Formeln bilden die Grundlage fiir ein Jterations- 
verfahren zur Erzielung vorgegebener Geschwindigkeitsverteilung und zur Berechnung der zu- 
gehoérigen Profile. 

Der Rechenaufwand dieses Verfahrens ist jedoch schon bei einem einzelnen Iterationsschritt 
recht erheblich. AuBerdem wurde die Iteration bis jezt nur fiir symmetrische Profile bei verschwin- 
dendem Ansteliwinkel durchgefiihrt. Da aber die berechneten Profile nicht nur bei einem einzigen 
Anstellwinkel bzw. Flugzustand giinstige Geschwindigkeitsverteilungen haben sollten®, muf ten 
parallel zur Iteration immer noch die Geschwindigkeitsverteilungen bei einigen nicht verschwindenden 
Anstellwinkeln mit verfolgt werden, was wieder zusatzlichen Aufwand erforderte. Die Berechnung 
gewolbter Profile ware mit diesem Verfahren zwar prinzipiell auch méglich gewesen, hatte aber 
noch ein Vielfaches an Aufwand gegeniiber den symmetrischen Profilen erfordert. Deshalb wurden 
die unsymmetrischen Profile nicht mehr mit der konformen Abbildung behandelt, sondern durch 
lineare Uberlagerung einer Wirbelverteilung gewonnen. Diese Naherung ist fir geringe Wol- 
bungen gut brauchbar, bei starkeren kommen jedoch zu groBe Ungenauigkeiten ins Spiel, wie wir 
hier spater noch sehen werden. 

Insgesamt waren das von Theodorsen entwickelte Verfahren als erstes seiner Art und noch mehr 
die in groBer Zahl damit erzielten Ergebnisse sehr wertvoll; es bestand aber schon lange der dringende 
Bedarf nach einem bequemeren Verfahren, das méglichst auch noch die unsymmetrischen Profile 
exakter erfassen sollte. 

Die erste Forderung fiihrte zu einigen konsequenten Naherungsverfahren. Das neueste davon 
diirfte dasjenige von EF. Truckenbrodt® sein, das an einige Arbeiten von F’. Riegels? anschlieBt. 
Dieses Verfahren ist im Aufwand sehr ertraglich und hat sich im praktischen Gebrauch gut bewahrt. 
Abgesehen von der unkontrollierbaren Giite der Naherung ist jedoch mathematisch unbefriedigend, 
da es von vorn herein nur von einer frei gewahlten Geschwindigkeitsverteilung ausgeht. Da eine 
solche prinzipiell nicht ganz willkiirlich ist, sondern zu einem geschlossenen Profil gehéren muB, 
kénnen dadurch Unvertraglichkeiten in Form von Selbstiiberschneidungen oder unbeabsichtigt 
stumpfen Hinterkanten auftreten. Diese miissen durch Probieren oder durch praktische Erfahrung 
so klein gehalten werden, da} sie ohne weitere Rechnung vollends korrigiert werden kénnen. Des- 
gleichen ist man auch beztiglich der Dicke des entstehenden Profils auf die Erfahrung angewiesen, 
da man auch diese Grobe der vorgegebenen Geschwindigkeitsverteilung nicht ansehen kann. 


F. Weinig, Z. angew. Math. Mech. 9 (1929), S. 507. 
A. Betz, Luftfahrtforschung 11 (1934), S. 158. 
Th. Theodorsen, Theory of Wing Sections of Arbitrary Shape. NACA-Report No. 411 (gsi): 
Th. Theodorsen, Airfoil-Contour Modification Based on e-Curve Method of Calculation Pressure Distri- 
bution. NACA ARR No. L4 GO5, 1944. Verkiirzte Wiedergabe in I. H. Abbott, A. E.von Doenhoff and 
L.S. Stivers, Summary of Airfoil Data. NACA-Report No. 824 (1945). 

® Die physikalischen Griinde hierfiir sind in NACA-Report 824, a. a. O. dargestellt. 

6 E. Truckenbrodt, Ing.-Arch. 19 (1951), S. 365. 


7 F. Riegels, Ing.-Arch. 16 (1948), 5. 373 (I. Mitteilung), 17 (1949), S. 94 (II. Mitteilune 
Saag (Berichtigung). ! asec 
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In der vorliegenden Arbeit wird nun wieder mit konformer Abbildung exakt gearbeitet, es wird 
jedoch der Umweg iiber die kreisahnliche Figur des Theodorsen-Verfahrens eingespart. Dadurch ver- 
einfachen sich alle Formeln so sehr, da® im Endeffekt kein wesentlich groBerer Rechenaufwand 
mehr notig ist, als beim Truckenbrodtschen Naherungsverfahren. Die Méglichkeiten der Geschwin- 
digkeitsvorgabe sind noch allgemeiner als dort, auch bei unsymmetrischen Profilen. Die Schlie- 


Bungs- und Vertraglichkeitsbedingungen werden sorgfaltig erfiillt und die Profildicke auf einfache 
Weise reguliert. 


2. Der Ansatz und die grundlegenden Formeln. Die konforme Abbildung, die zur Berechnung 
der unendlich ausgedehnten Parallelstr6mung um eine Kontur benétigt wird, bildet das 
AuBere dieser Kontur auf das 
AuBere des Kreises ab. Soll die ie 4 
Kontur in der z-Ebene (Abb. 1) 
mit 

z=x+1y (1) 
liegen und die Kreisebene mit 
€ bezeichnet werden, dann mu 
fiir die Abbildungsfunktion z(¢) 


BGO) ==100" 5 


Pa = 00) = reell (2) og V=1 
gelten, damit in z- und €-Ebene 
die gleichen Anstrémbedingun- 
gen vorhanden sind. Ist die 
Abbildungsfunktion 2(€) be- 
kannt, dann kann man jede 
Strémung der €-Ebene auf die 
z-Ebene iibertragen. Soll in 
der z-Ebene eine Strémung 
vom Anstellwinkel ~« um die 
Kontur entstehen, dann mu 
wegen (2) als Grundlage in der 
€-Ebene eine ebensolche Stré- 
mung um den Kreis genommen 
werden. In der z-Ebene ver- 
laBt die Strémung das Profil 


an der Hinterkante, wodurch AY, le 

eine Zirkulation J‘ und ein Auf- 

trieb bedingt ist. Wir kénnen : @ 
= a Sd 


das Profil der z-Ebene so ge- 
dreht voraussetzen, dai die 
Funktion 2z(€) unter den Be- 
dingungen (2) die Hinterkante Abb. 1. Skizze der eingefiihrten Ebenen. 

des Profils nach =1 abbildet. 

Dann muf in der €-Ebene der hintere Staupunkt bei € = 1 liegen. Wir erreichen dies in bekannter 
Weise mit dem komplexen Potential 


FQ) =O +i = Clee +o) ais 


Qi 


—_ 
oe 
aes 


der Parallelstrémung vom Anstellwinkel a um den Kreis mit der tiberlagerten Zirkulation 

F=Aasina-C. (4) 
Mit der Abbildungsfunktion 2(¢) und (3) hat man das komplexe Potential der z-Ebene als mittel 
bare Funktion und kann damit dort alle interessierenden GréBen berechnen. Fiir z(¢) machen wir 
ohne weitere Zwischenabbildung direkt den Ansatz 


= 3 (C4 :) 5 (a, +b) C-* (5) 


wobei die Potenzreihe in |¢| > 1 gleichmaBig konvergent sein soll. Der Ansatz (5) erfillt sicher die 
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Bedingungen (2). Sind alle a,, b, = 0, dann entsteht in der z-Ebene als Bild des Kreises der gerad- 
linige Schlitz bzw. die Strémung um die ebene Platte. Nicht verschwindende a,, b, ergeben gewisse 
Abweichungen von der ebenen Platte. Sind sie nicht groB, dann entstehen flache Profile. 

Nehmen wir zunachst an, die a,, b, seien bekannt, dann folgt aus (5) mit € = e'” fiir das in der 
z-Ebene entstehende Profil die Parameterdarstellung 


x = cosy—> (a, cos vy + 6, sin yg) , (6) 
v=0 
y= > (a, sinvg—b, cosyg). (7) 


o 


i 


Ein zur x-Achse symmetrisches Profil entsteht, wenn alle b, = 0 sind. Beim Ubergang von P nach 
— andert dann y sein Zeichen, wahrend x erhalten bleibt. Der Geschwindigkeitsbetrag v in der 


z-Ebene ist allgemein 


ie 
adF|_ (|d¢ ; (8) 
~ | GB | Wak 

dt 


Wir wahlen in (3) die Konstante C = i , dann wird die Anstrémgeschwindigkeit 


V=volo)=1, (9) 
und die Geschwindigkeit v ist identisch mit dem sonst haufig gebraéuchlichen Geschwindigkeits- 
verhaltnis v/V. Fiir die Geschwindigkeit am Profil setzen wir z = e'? und erhalten aus (8) mit (3) 
und (5) 


Pxcna — eit e—2ip)+isina e—i¢ 


(10) 


vel?) = 


dx + idy 
ie? dp 

Wir kennzeichnen die Ableitung nach gy mit einem Punkt, nach x mit einem Strich, die anderen 

Ableitungen werden als Differentialquotienten ausgeschrieben. Dann folgt aus (10) 


sin (p — a) + sino (11) 
ee 

Die Betragstriche lassen wir bei dieser reellen Gré®Be weg. Sie bekommt damit einfach das Vor- 

zeichen des Zahlers. Die Wurzel bleibt positiv. Formen wir danach (11) noch um zu 


v(e ?) —— 


sin (~—4)+ sina [sin (py —&) + sina] cosO 
eel aa = i : 12 
ae V1 + y®? (— 49) eat) a 


dann muf die Wurzel und deren Reziprokwert cos 0 (9 = Steigungswinkel des Profils nach Abb. 1) 
das Vorzeichen von — x(p) haben. Das Minuszeichen ist dabei unwesentlich. Es wird eingefiihrt, 
damit cos © auf der Oberseite des Profils positiv wird. 

Die Gleichungen (6), (7) und (11) oder (12) geben bei bekannten a,, b, Profil und Geschwindigkeit 
in Abhangigkeit von ¢. 

In (11) ist der Nenner unabhangig von «. Andert man also in der ¢-Ebene a, dann muB sich v 
nach 


~, sin (y—a)+ sina 
(9 x) sin (py — a) + sind 


u(y, 4) (13) 


ebenfalls andern. In diese Gleichung geht das Profil nicht ein. 
Der Auftrieb in der z-Ebene ist nach dem Kuttaschen Satz, wenn o die Dichte des str6menden 
Mediums ist, 
A=oVI, (14) 


wobei J’ auch in der z-Ebene den Wert (4) hat, da sich die Zirkulation bei der konformen Abbildung 
nicht andert. Der Auftriebsbeiwert ist dann 


A : 1 
= 4asina-—. 
2 yy (15) 


Cg = 
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Die Lange | des Profils hat dabei nach (6) bei der ebenen Platte den Wert | = 2, bei dickeren Pro- 
filen wird sie etwas kleiner sein, da dort der Auftriebsbeiwert allgemein etwas gréfer ist!. Der Auf- 
trieb verschwindet bei a = 0. Die fiir die Abbildung der Hinterkante getroffene Festlegung he- 
wirkt also eine zur x-Achse parallele Nullauftriebsrichtung. 


Fur das Moment My der Luftkrafte bei verschwindendem Auftrieb erhalten wir ebenfalls in 
einfacher Weise, z. B. nach H. Glauert?, wenn wir es in mathematisch negativem Sinn positiv zahlen, 


M, =<: 2b, (16) 
und fur den entsprechenden Momentbeiwert 
: Paohles 226, 
moQ-_ [2 2 (17) 


Die abgeleiteten Formeln sind nicht neu. Beispielsweise haben F. Riegels und H. Wittich® schon 
versucht, mit ihnen das Problem der vorgegebenen Kontur ohne die Theodorsensche Zwischenab- 
bildung zu lésen. Dieser Versuch brachte allerdings nur bei symmetrischen Konturen Gewinn. 
Bei den unsymmetrischen ist die Konvergenz des auf (6) und (7) aufgebauten Iterationsverfahrens 
fraglich und schlecht, so daB das Theodorsen-Verfahren hier allgemein vorzuziehen ist. Die Schwie- 
rigkeiten bei der Bearbeitung des vorgegebenen Profils haben F. Riegels dazu veranlabt, Gleichung 
(12) fiir symmetrische Profile zu linearisieren, was einem Quell-Senken-Verfahren gleichkommt?. 
Dieses wurde mittels Wirbelverteilungen und Korrekturen auf unsymmetrische Profile ausgedehnt® 
und diente auch E. Truckenbrodt als Grundlage zu seiner Bearbeitung der vorgegebenen Geschwindig- 
keitsverteilung ®. 

Wir werden im folgenden erkennen, daB® die gesamten Schwierigkeiten, welche die Einsparung 
des Theodorsenschen Umwegs iiber die kreisahnliche Figur beim Problem der vorgegebenen Kontur 
mit sich bringt, bei der Umkehrung vollkommen wegfallen. Bei der Aufgabe der vorgegebenen 
Geschwindigkeit ist eine Linearisierung iiberfliissig. Die Einsparung der kreisahnlichen Figur, die 
beim Problem der vorgegebenen Kontur andere Nachteile mit sich brachte, wirkt sich jetzt voll 
aus. Man spart sehr viel Aufwand und hat bei der Vorgabe der Geschwindigkeit noch besonders 
vorteilhafte Méglichkeiten. 


3. Die Vorgabe der Geschwindigkeitsverteilung. Wenn wir nun die abgeleiteten Formeln dazu 
beniitzen, Profile mit giinstigen Geschwindigkeitsverteilungen zu berechnen, dann wollen wir von 
vorn herein auf das numerische Endergebnis hinarbeiten und uns demnach auch ganz von den 
praktischen Anforderungen leiten lassen, die an dieses Ergebnis gestellt werden. Es geht darum, die 
in der vorangehencen Ziffer eingefiihrten Koeffizienten a,, 6, so zu bestimmen, daf} durch (6) 
und (12) eine vorgegebene Geschwindigkeitsverteilung dargestellt wird. Fiir diese Aufgabe bietet 
Gleichung (12) einige ausgezeichnete Méglichkeiten. 

Wir folgern aus (12) und (6) 
sin (y — 4) + sina 


— x(y) = sing = (va, sin yy— yb, cos vq) = mG) . (18) 


wobei wir 

o(y) =V1+y? V9) (19) 
eingefiihrt haben. Wahlen wir in (18) einen Anstellwinkel « = a* und eine Funktion w(p), dann 
kénnen wir durch harmonische Analyse sofort Koeffizienten a,, b, bestimmen, die mit (6), (7) und 
(12) ein Profil und eine Geschwindigkeitsverteilung liefern. Bei dem in (18) gewahlten a* ist diese 
nach (12) speziell as 

7% \ ee 2) 

olga") = ae (20) 
bei einem anderen a folgt sie hieraus mittels (13). Fiir die GréBe y’ in der Wurzel ist x direkt durch 
(18) gegeben, y aus den a,, b, mittels 


¥(~) = Sy (va, cosvgy + vb, sin v@) (21) 
y=1 


zu berechnen. 


1 W. Weinberger, Luftfahrtforschung 17 (1940), S. 3. 

2 H.Glauert, Die Grundlagen der Tragfliigel- und Luftschraubentheorie, Berlin 1929. 

3 F. Riegels, H. Wittich, Zur konformen Abbildung beliebiger Profile auf den Kreis und Bestimmung der 
Druckverteilung in ebener inkompressibler Stromung, Bericht 42/A/20 der AVA Gottingen, 1942. 

4 F, Riegels, FuBnote 7 von S. 436, I. Mitteilung. 

5 F. Riegels, FuBnote 7 von S. 436, II. Mitteilung. 6 E. Truckenbrodt, FuBnote 6 von S. 436. 
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Es fragt sich nun nur, wie man durch die Wahl von co(g) erreichen kann, da die in Parameter- 
form (6) und (20) dazu entstehende Geschwindigkeitsverteilung eine im Voraus festgelegte Form 
erhalt. Diese Frage ist leicht zu beantworten, solange wir es mit relativ flachen Konturen zu tun 
haben, wie es Fliigelprofile praktisch immer sind. Bei ihnen sind wegen (7) die a,, b, klein gegen 1, 
also ist nach (6) 

x(p)  cosp. (22) 

Zudem ist y’ aufer im Bereich der Profilnase klein, wir haben also fast am ganzen Profil nach (20) 

v(p, a*) ~ o(¢) - (23) 

Durch (22) und (23) ist schon eine Naherung fiir die Verteilung v(x, a*) gegeben. Gleichung (6) be- 

wirkt demgegeniiber noch eine kleine Verschiebung in x-Richtung, die man nach wenigen Beispielen 

gut von vorn herein beurteilen kann. Der EinfluB des Wurzelfaktors in (18) ist qualitativ ebenfalls 

klar. Er ist nur in der Umgebung der Fliigelnase stark. Dort bewirkt das groB werdende y’ 

ein Absinken der Geschwindigkeit, das dort ohnehin wegen der Nahe des vorderen Staupunkts 

vorhanden ist. Die durch die Wahl von w(q) entstehende Geschwindigkeitsverteilung ist also von 
vorn herein weitgehend zu iibersehen. 

Wir gehen also von der Funktion ~(y) aus und berechnen daraus ein Profil und die zu ihm (im 
Rahmen der noch zu klarenden Rechengenauigkeit) exakt gehérende Geschwindigkeitsverteilung. 
Weicht sie zu sehr von dem ab, was wir urspriinglich beabsichtigt haben, dann kénnen wir immer 
durch einen kleinen Zusatz in der Funktion @(p) noch Korrekturen anbringen. 

Man kénnte so ein ganzes Iterationsverfahren aufbauen. Praktisch kommt es jedoch immer nur 
an wenigen Stellen des Profils auf eine Verbesserung der entstehenden Verteilung an, etwa da 
wo konstante Geschwindigkeit entstehen soll, in deren Bereich dann die Stabilitat der laminaren 
Grenzschicht sehr empfindlich ist. In allen anderen Fallen kann man sich mit der entstehenden, 
nach (22) und (23) von vorn herein ungefahr festzulegenden Geschwindigkeit begniigen. Wir be- 
schreiben deshalb kein Iterationsverfabren, sondern beschranken uns auf eine Korrektur, die bis- 
her stets ausgereicht hat. 

DaB die Geschwindigkeit v(x) auf diesem Weg nicht direkt vorgegeben wird, ist kein Nachteil 
des Verfahrens, denn sie kann ohnehin nicht ganz beliebig vorgegeben werden. Sie mu zu einem ge- 
schlossenen Profil gehéren, also noch gewisse Bedingungen erfillent, damit aus ihr iberhaupt ein 
sinnvolles Profil berechnet werden kann. Man muf also freie Parameter in die Rechnung mitnehmen, 
die nachtraglich festgelegt werden, und die somit auch Abweichungen von der urspringlichen 
Vorgabe mit sich bringen. 

Bei unserem Vorgehen miissen sich die genannten Bedingungen auf die Ausgangsfunktion w(~) 
iibertragen. Es wird sich als besonderer Vorteil erweisen, da sie hier sehr einfach werden und 
leicht zu erfiillen sind. 

Selbstverstandlich miissen wir von w(y), beziehungsweise von der in (18) damit gebildeten 
Funktion 


sin (py — a) + sin a 
h(g) =. (24) 

eine gewisse Glattheit verlangen. Wir bleiben geniigend allgemein und vermeiden alle weiteren 
Schwierigkeiten, wenn wir verlangen, dafs h(p) stiickweise analytisch ist und als héchste Singu- 
laritaéten Knicke, d. h. Sprungstellen in der ersten Ableitung haben soll. Damit werden die a,, b 
von der Ordnung 1/y? ?, und alle auftretenden Reihen konvergieren gleichmafig. i 

Die ganze Stetigkeitsforderung wird an sich spater illusorisch, wenn wir von den Fourier-Reihen 
auf endliche Polynome tibergehen, welche h(y) durch eine iiberall analytische Funktion inter- 
polieren. Praktisch ist aber die Glattheit von h(q) auch dann noch wichtig, weil von ihr die GréBen- 
ordnung der vernachlassigten héheren a,, 6, und damit die Zahl der bei geniigender Genauigkeit 
mitzunehmenden Koeffizienten abhangt. Auch hierfiir wird sich unsere Forderung als ausreichend 
erweisen. 

Weitere Bedingungen miissen sich vor allem auf die SchlieBung des entstehenden Profils bezichen. 
Eine solche Bedingung ist bei uns in (18) schon enthalten. Dort beginnt die Reihe mit dem Glied 
y = 1, also mu das Absolutglied 


v 


i bers h(y) dp = 0 (25) 


1 Vergleiche beispielsweise auch A. Betz, FuBnote 2 von S. 436. 
® H. v. Mangoldt — K. Knopp, Einfiihrung in die hohere Mathematik, Bd. 3, S. 508, Leipzig 1944, 
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sein. Ware (25) nicht erfillt, dann wiirde das durch Integration aus x(~) zu gewinnende x(~) nicht 
periodisch, das entstehende Profil wiirde sich nicht schlieBen. Wir nennen deshalb (25) die 
SchlieBungsbedingung. 

Diese Bedingung kann aber noch nicht alles enthalten, denn sie ist bei symmetrischen Profilen, 
bei denen x(p) eine ungerade Funktion ist, automatisch erfiillt, und auch dort ist die Geschwindig- 
keitsverteilung nicht willkiirlich. Es liegt nahe, ebenso wie bei x(p) auch bei y(~) die Periodizitat 
nachzupriifen. Sie erweist sich aber gleichzeitig mit derjenigen von x(q) wegen (6) und (7) als er- 
; fiillt. Fur die SchlieBung des Profils reicht also (25) aus. Dagegen kann in einer anderen Hinsicht 
bei der Wahl von «(q) noch eine unbeabsichtigte Wirkung auftreten. Wir wollen normalerweise 
an der Hinterseite des Profils eine Spitze erzielen, an der die Strémung glatt, also mit einem end- 
lichen, nicht verschwindenden v abflieBen kann. Entsprechend werden wir (0) + 0 vorgeben. 
In (18) erkennen wir aber, da8 trotzdem immer x(0) = 0 ist, waihrend im allgemeinen y(0) + 0 
wird. Fiir das Gegenteil ist jedenfalls keinerlei Ursache vorhanden. Wir erhalten also bei g = 0im 
allgemeinen y= 0, und nach (19) v= 0. Dies gibt genau den Staupunkt, den wir vermeiden 
wollten. Er kann sogar leicht noch so liegen, da®B die senkrechte Tangente im Rahmen einer 
 schleifenférmigen Selbstiiberschneidung des Profils auftritt. Um dies alles zu vermeiden, miissen 
wir bei y = 0 aufer dem endlichen ~(0) auch noch 


5(0) = Siva, =0 (26) 


vy=1 


verlangen. Wir nennen dies die » Spitzenbedingung“. 


_ Das Aussehen der Hinterkante bei Beriicksichtigung der Spitzenbedingung wird durch die GréBe 
¥(0)/x(0) bestimmt. Ist die entstehende Kontur symmetrisch, dann verschwindet (0) als ungerade 
Funktion, wenn ihre Reihe gut genug konvergiert, insbesondere wenn sie nur endlich viele Glieder 
hat. Wir bekommen dann verschwindenden Hinterkantenwinkel, wie es bei dem nicht verschwin- 
denden ~(0) zu erwarten ist. 

Bei unsymmetrischen Konturen kann y'(y = 0) trotz (26) noch unendlich werden. Hat namlich 
h(g) bei g = 0 einen Knick, was durch eine Sprungstelle in @(q) entstehen kann, dann bekommt 
x(q) bei p = 0 eine Sprungstelle und ‘y(~) eine logarithmische Singularitat. Wegen y = 0 entsteht 
aber keine stumpfe Hinterkante, sondern eine senkrecht stehende Spitze am Profil. Man konnte 
dies vermeiden, indem man fiir w(g) bei py = 0 Stetigkeit verlangt. Praktisch ist dies aber nicht 
notig, weil die entstehende Singularitat sehr schwach ist und nur eine kleine Umgebung der Hinter- 
kante einnimmt. Vor allem, wenn man nachher mit endlichen trigonometrischen Polynomen ar- 
beitet, ist sie effektiv nicht mehr zu spiiren. Man nimmt sie in Kauf, ebenso wie dies bei der 
Theorie der Wirbelverteilungen geschieht. 

Die Bedingung (26) ist genau das, was bei Truckenbrodt' nur ungefahr durch die Erfahrung er- 
fillt wird. Man kénnte diese Forderung auch dort noch einfiihren. Allerdings erfordert dies zu- 
satzlichen Aufwand, so daB bei Beriicksichtigung der Spitzenbedingung das dortige Verfahren nicht 
mehr wesentlich einfacher ist als das vorliegende. 

Die bis jetzt abgeleiteten Bedingungen fiir Stetigkeit, SchlieBung und Spitze garantieren diese 
speziellen Eigenschaften. Streng genommen kénnen unter diesen Bedingungen immer noch un- 
mégliche Profilformen entstehen. Mit einigem Geschick kann man ~(p) beispielsweise so vorgeben, 
daB ein Profil mit negativer Dicke entsteht. Diese Méglichkeiten schlieSen wir aber nicht mehr 
gesondert aus. Sie sind im praktischen Gebrauch leicht zu vermeiden. Wir gehen also grundsatzlich 
von a* und @(qy) aus, priifen nur mittels (25) und (26) SchlieBungs- und Spitzenbedingung nach 
und arbeiten dann in der eingangs dieser Ziffer beschriebenen Weise. 

Nach (14) und (15) ist durch den gewahlten Anstellwinkel o* der Auftrieb und der Auftriebs- 
beiwert c, festgelegt. Man wahlt also mit a* von vorn herein den Auftrieb, den die aus w(g) ent- 
stehende Geschwindigkeitsverteilung v(x, «*) liefern soll. Gleichung (25) zeigt dann, ob die Funk- 
tion «(q), die nach (23) ungefahr mit v(p, a*) iibereinstimmt, richtig gewahlt ist. Ist beispielsweise 
a* positiv, dann mu die Geschwindigkeit v im Mittel auf der Unterseite kleiner sein als auf der 
Oberseite. Dies wird durch (25) automatisch geregelt. Wir kénnen also leicht ein Profil entwerfen, 
dessen Geschwindigkeitsverteilung bei einem gewissen c, vorgegebene Eigenschaften hat, was 
praktisch oft verlangt wird. 

Wir haben aber in unseren Formeln eine noch niitzlichere Méglichkeit. In (18), (24) und (25) 
ist grundsatzlich nur von den dort stehenden Funktionen von @ Gebrauch gemacht. Rein formal 


1 K. Truckenbrodt, FuBnote 6 von 5. 436. 
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kénnen wir in diesen Gleichungen auch «* als Funktion von p einfiihren, ohne am Wesen dieser 
Gleichungen und damit am ganzen Verfahren etwas zu andern. Praktisch wichtig ist dabei nur 
eine stiickweise konstante Funktion a*(g). Wir kénnen sie trotz unserer Stetigkeitsforderung 
zulassen, da wir an den Sprungstellen von «*(y) mit Sprungstellen von w(p) stetiges h(y) erzeugen 
kénnen. Dies wird spater keine Schwierigkeiten machen. 

Die nicht iiber das ganze Profil konstante Wahl von a* erscheint zunachst vielleicht absurd, 
mindestens zwecklos. Sie bedeutet aber natiirlich nicht, daB das entstehende Profil gleichzeitig 
unter verschiedenen Anstellwinkeln angestrémt wird. Vielmehr wird auch mit «*(p) nach (18) eine 
konkrete Funktion x(p) und damit ein festes Profil definiert. An diesem Profil kénnen wir dann 
nach (12) mit der Funktion x(q) fiir jedes (selbstverstandlich am ganzen Profil konstante) a eine 
Geschwindigkeitsverteilung berechnen. Diese ist, wenn wir x(p) nach (18) in (12) einsetzen, 


sin (p — 4) + sina o(() 27 
v(p, 0) ~ sin [y —a*(~)] + sin a*(¢) /1- aa y? ; a) 


Stimmt das gewahlte a in einem g-Bereich mit a*(p) iberein, dann ist dort (27) identisch mit (20). 
Dort stimmt also v(p,a) ungefahr mit w(q) iiberein. Da a nach (14) wieder den Auftrieb fest- 
legt, kénnen wir es somit erreichen, da bei einem vorgegebenen c,y langs eines bestimmten Profil- 
stiicks (beispielsweise der Oberseite mit 0 <y <q, ~ 2) die Geschwindigkeit v(p, c,) ungefahr 
mit der dortigen Vorgabe «(q) iibereinstimmt, bei einem anderen c,, die gleiche Eigenschaft langs 
eines anderen Profilstiicks (z. B. der Unterseite) auftritt. Unsere Bedingungen (25) und (26) sorgen 
auch hier wieder automatisch dafiir, daB alle nach (27) berechneten Geschwindigkeitsverteilungen 
mit den jeweils durch (14) dazu festgelegten Auftrieben vertraglich sind. Diese Art der Geschwindig- 
keitsvorgabe, die bisher noch nirgends méglich war, wird sich als sehr wertvoll erweisen. Einige 
typische Anwendungen sind in Ziff. 5 enthalten. 


4, Die praktische Durchfiihrung. Fiir die praktische Durchfiihrung unseres Verfahrens bené- 
tigen wir vor allem Rechenformulare zur harmonischen Analyse und Synthese, da wir diese Rechen- 
gange mehrmals durchfiihren miissen. Wir nehmen fiir beide Operationen das Rungesche Verfahren. 
Bei ihm wird die zu analysierende Funktion im Intervall (0,2 2) an M dquidistanten Punkten be- 
rechnet, wobei M méglichst durch 3 und 4 teilbar sein soll, und die M = 4 N Ordinaten P,, werden 
durch ein trigonometrisches Polynom 


2N 
P(y) = » (d, cos yg 4+ ¢, sin y¢) (28) 
mit 4 N Koeffizienten interpoliert. Die Koeffizienten berechnen sich dabei mittels 
4N 
1 am 
4, = 3H 2s Pa c0s(v 7) (year. ee Ve 
4 
1 
dy ss ZN om Pn : 
; (29) 
1 
daw = Fy Dy Pm (— 1)" 
4.N 
te CLE : mm 
C, IN Py sin (99) (vy 1, 28N) 


also durch einfache Produktsummenbildung. Dabei kann man wegen der Periodizitat der tri- 
gonometrischen Funktionen eine Menge von Faltungen durchfiihren, so da beispielsweise bei 
N = 12 die langste zu berechnende Produktsumme anstatt 48 nur 6 Glieder hat. Fiir die gesamte 
Operation existieren fertige Rechenformulare, so da® wir hier nicht darauf eingehen!. Der Rechen- 
aufwand ist dabei sehr ertriglich. Bei N — 12, wie die Beispiele berechnet wurden, ist eine voll- 
stindige Analyse einschlieBlich kompletter Proben in etwa 4 Stunden leicht zu schaffen. Ganz 
ahnlich verlauft die Synthese einer Funktion P(g) an 4 N aquidistanten Punkten. 

Die Anwendung dieser Methoden bedeutet einerseits, daG man an Stelle der vorgegebenen Funktion 
hi) in (18) mit einem trigonometrischen Polynom arbeitet, das an 4 N Punkten mit h(p) tiberein- 
stimmt, andererseits, dafS man Profil und Geschwindigkeitsverteilung zu diesem Polynom nur an 


1 A. Hufmann, Rechnerische Verfahren zur harmonischen Analyse und Synthese, Berlin 1938. 
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4 N Punkten berechnet. Beides gibt mit N = 12 sicher geniigende Genauigkeit. Es ist aber keiner- 
lei Schwierigkeit, auch mit anderen N-Werten zu arbeiten. 

Die Verwendung der erwahnten Rechenformulare ist kiinftig zu verstehen, wenn wir von Ana- 
lyse oder Synthese sprechen. Die hauptsichliche Rechenarbeit ist mit ihnen, wie wir sehen werden, 
zu bewaltigen. 

Die Funktion w(y) setzen wir, um SchlieBungs- und Spitzenbedingungen erfillen zu kénnen, 


in der Form 
K 
1 
= ky; 
mG) > vil) 


an. Dabei sind zur Erfiillung von (25) und (26) mindestens K = 2 Parameter k; notwendig. Wir 
nehmen aber meist noch mindestens einen weiteren mit, um die Dicke des Profils und andere 
Eigenschaften regeln zu kénnen. Die y,(~) miissen so gewahlt werden, da die durch (30) definierte 
Funktion w(y) die gewiinschten Figenschaften bekommt und da zu einer unstetigen Wahl a*(~) 
eine stetige Funktion h(y) nach (24) entsteht. Wie dies im Einzelfall geschieht, ist an den Bei- 
spielen der nachsten Ziffer zu erkennen. 
Gehen wir mit (30) in (18) ein, so folgt 
2N 


(30) 


K 

sing — )») (va, sinyp— vb, cosyq) = {sin |p—a*(y)] + sina*(p)} Sy ki yi(y) - (31) 
y=1 i 

Wir fiihren die Analysen der den einzelnen k; proportionalen Anteile getrennt durch. Wir setzen 


2N 


2 (dri cos yp + c,; siny gy) = {sin [p — a*(y)] + sin a*(y)} yi(y) = hie) 


und berechnen die so definierten d,;,c,; durch harmonische Analyse. Damit werden die a,, b, 
wegen (31) 


(32) 


K 
l—a,= Ski; , (33) 
ie i=l 
=, DS key; (le Sater arlN. icy (34) 
pal 
K 
cad ts Pe (yea le 2. ZN). (35) 
Fa 
Aus (25) und (26) folgt dann 
K 
> kpdg: = 0 (36) 
=r 
und 
2N K K 
—= Skio¢: t1l— Ska, =0. (37) 
y=2 i=1 j= 
Mit 
2N 
§; = De Cyi (38) 
v=l1 
folgt aus (37) r 
os k; 5 = i e (39) 
i= 


Die Gleichungen (36) und (39) beniitzen wir dazu, zwei von den Parametern k; in den iibrigen aus- 


zudriicken, etwa mittels 


K 
tA, i k; p) 
ay (40) 
k, = B,+ » Bk; , 
i=3 
wo die 
A = dos pee So doi —— si do, 
08; dog — $y.4o1° ‘ 51d) 2— S2do1 : 
—d si dy, — 8, doi (41) 
B 01 B; s= | 


= ° 
$1 dy — Sy doy 


$1 yg — 82 do 


30 
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aus den durch die Analyse berechneten GroBen d,; und s; leicht zu bilden sind. Damit lassen sich 
die a,, b, aus (33) bis (35) in der Form 


K 
sane Shab | 
to 


K (42) 
b, = 1) 4+ Dk; be , | 
3 
mit 
al) Sly By ee) 
— ya = Aje,, + Bo cys (pean Sane) 
— a= A; ¢, + Beeia a ey: (vy == 1, aoe ee) s Pi (43) ° 
yb = Ayd,, + Bod,» (PN eee i 
y b@) =A; d,, +B; d,. 44,5 (vee; 2 IN) 


darstellen. Man kann also auch bei der Synthese der Funktionen x(y) und y(¢g) die den tberzahligen 


Parametern k; proportionalen Anteile getrennt summieren. Wir setzen 


2N 
x(p) = cos p— ¥ (al cos vp + 5 sin vg), 
v=0 


2N 

x (9g) = 2 (a cos vg + b sin vq) (i= 3,...K), (44) 
2N A . 

yp) = (ai) sin vy p — Df) cos v ~) (=.0, Sree ce 4) 


v=0 


und haben dann nach (6), (7) und (42) 


x(p) = #(g) — Sky (g) , 
a (45) 
MD OP) aia Ie yg) . 


Man nimmt also die tiberzahligen Parameter bis ins Ergebnis mit. Man bekommt auf diese Weise 
ein Grundprofil x(q), (qr) und kann dazu mit den k; jeweils noch ganze Scharen von (auto- 
matisch geschlossenen und spitzen) Profilen zusatzlich berechnen, ohne weitere Analysen und Syn- 
thesen durchfiihren zu miissen. Dies ist vor allem giinstig, wenn man bei Profilserien die Dicke 
verdndern will. 


Die Geschwindigkeitsverteilungen solcher Profilserien sind zwar nicht mehr linear in den k;,, 


aber sie lassen sich ebenfalls einfach berechnen. Da an den Profilen immer die Druckverteilung 
interessiert, fiir die nach der Bernoullischen Gleichung 


Ee Nec ee (46) 


gilt, geht man sofort auf v°(p) zu, wofiir nach (11) mit (18) 


1 
U2 Po *(p) == 
EON Dice ore = FEC (47) 
orp) ' {sin (p—a*(p)] + sina*(y)}? 
gilt. Dabei gibt der zweite Term des Nenners die bei kleinem y’ kleine Abweichung zwischen w 


und v. Fir die praktische Berechnung von 1/wkann man entweder nach (40) und (41) die k; aus- 
rechnen und in (30) eingehen, oder man kann auch hier Funktionen 


yO = Ayy, + Boye, vy) = Ary, + Biyg +4; (48) 


berechnen, mit welchen 


1 = 
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wird. Analog setzt man die Funktion ¥(g) aus einzelnen Anteilen zusammen, indem man 


2N 
yO = Sv al) cosy + v b®sin y 7) (50) 


ill 


einfuihrt und am besten gleich Funktionen 
rey ee pyre k , 
oO) — nlp olan a) 


bildet, wenn man mit mehreren Werten k; rechnen will. Der Nenner von (51) ist dabei schon friiher 
berechnet. Mit 


(51) 


K 
g(p) = g(p) + ak g(¢) (52) 
erhalt man dann v*(~,a*) nach (47) in dem einfachen Rechnungsgang 
5 1 
M0) <7 53) 
(ae) + 8°(9) 


Die gesamten Funktionen mit oberem Index i hangen nur von y;, y, und y, ab. Dies 1a8t sich 
an den 4,, B; leicht verfolgen. Will man also einmal noch eine zusatzliche Funktion YK +1 ein- 
fiihren, dann andern sich alle vorangehenden Funktionen nicht, man hat lediglich neue Funktionen 
aK+1),y(+1) und g(k+1) zu berechnen und kann dann die Zusitze in x((), y(g~) und im Nenner von 
v?(p,o*) anbringen. Es macht dabei auch keine wesentliche Miihe, die vorangehenden Parameter 
ks ...kx noch zu andern. 

Die nachtragliche Einfithrung einer zusatzlichen Funktion YK +1 ist vor allem bei einer even- 
tuellen Verbesserung der entstehenden Verteilung u(y, «*) wichtig. Weicht etwa die zu einer Wahl 
der y; und k; mit i < K berechnete Verteilung u(y, a*), fiir welche nach (20) 


1 1 


OP) vp, a*) 14 y? (54) 
gilt, zu sehr von einer beabsichtigten Verteilung v(q) ab, fiir welche 
Bee ae 
Bp) o(y) * (ae) B(p) VI 5? ») 
gelten sollte, dann wahlt man, da y’ und y’ sicher wenig voneinander abweichen, 
A (xen) =yx (9) = (sey El oe (56) 


und geht damit in den Ansatz (30) ein. Entsprechend der Herleitung von (56) nimmt man an, 
daB die Korrektur yx 4, mit kx,,;= 1 am besten arbeitet. Man bemerkt aber im. praktischen 
Gebrauch, da®B damit die Korrektur gern tiber das Ziel hinausschieBt. Man kann dies mit 
kleineren Werten kx,; vermeiden, deren giinstigste GréSe man leicht durch Ausprobieren an 
wenigen Stellen ermitteln kann. 

Berechnet man symmetrische Profile, dann ist x(y) eine ungerade Funktion, und man muf sie 
nur im halben Intervall vorgeben. Dann verschwinden alle b, und d,, und man hat bei allen Ana- 
lysen und Synthesen reine Sinus- oder Cosinusprozesse, die durchweg nur etwa den halben Aufwand 
der gemischten erfordern. Die SchlieBungsbedingung (34) lauft in diesem Fall leer. Wir kénnen 
dann nur mittels (39) einen der Parameter k; in den anderen ausdriicken. Damit wir die abgeleiteten 
Formeln nicht andern miissen, setzen wir in diesem Fall y,(~) = 0 und driicken k; durch (39) in 
den k; mit i > 3 aus. Wir erhalten dann statt (41) 


ee te Ae Baers (57) 
il il 
alles andere bleibt vollstandig erhalten. 

Oft ist es auch bei unsymmetrischen Profilen mdglich, wenigstens einzelne y; als ungerade Funk- 
tionen zu wahlen und auch so Arbeit zu sparen. Alle diese Einzelheiten sind im praktischen Gebrauch 
leicht zu ubersehen. 

Es sei nur noch bemerkt, daB die Geschwindigkeit v(y,a*) nach (53) nicht unbedingt mit be- 
rechnet werden mu, wenn man etwa bei diinnen Profilen die Geringfiigigkeit des Wurzeleinflusses 
in (20) iitbersehen kann, oder wenn man auf hohe Genauigkeit keinen Wert legt. Die alleinige 


30S 
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Berechnung von Profilen oder Profilserien ist dann wirklich bequem zu erledigen. Bei Bedarf 
kann man dann v(p, 4*) immer noch nachtraglich genau berechnen, oder wenigstens den Wurzel- 
einflu8 graphisch mit dem in (12) eingefithrten Steigungswinkel 0 des Profils abschatzen. Man kann 
also das beschriebene Vorgehen bis zum Naherungsverfahren vereinfachen und von ihm aus auf- 
bauend sich erst den genauen Fehler verschaffen und ihn gegebenenfalls durch eine Korrektur be- 
seitigen. eee 

An den Beispielen der folgenden Ziffer ist nun der gesamte Rechnungsgang einschlieBlich der 
sinngerechten Wahl der y; nochmals zu verfolgen. 


5. Beispiele. A. Als erstes Beispiel wollen wir die Berechnung einer symmetrischen Profilserie 
durchfiihren, bei welcher auf der Profiloberseite bei einem Anstellwinkel von «* = — 0,9° von der 
Nase bis etwa 80%, der Fliigeltiefe konstante Geschwindigkeit herrschen soll. Der dahinter folgende 
Druckanstieg soll erst steiler, dann flacher sein. Der physikalische Hintergrund dieser Aufgaben- 
stellung, die mit zu dieser Arbeit angeregt hat, wurde an anderer Stelle beschrieben 1, z 

Um die konstante Geschwindigkeit zu erzielen, wahlen wir im Bereich der Oberseite a«*(y) =— 0,9 
und y,(~) = 1. Nur in der Nahe der Fliigelnase, also bei yp = a kénnen wir diese Festsetzung nicht 
gebrauchen, denn dort ware sonst nach (31) x(p =a + 20*) =0, was wegen der Symmetrie 
erst bei p =z eintreten soll. Bei p =a + 2 .a* liegt ein Staupunkt, wie man in (11) leicht fest- 
stellt. Dort muB vy = 0 und damit w = 0 vorgegeben werden. Dies geschieht am einfachsten, indem 
man dort direkt h(p) vorgibt, da diese GroBe dort stetig sein mu. Wir wahien also 


(t (0<@ < 165°), 
a Asin p (165° <p < 180°) (58) 


sin (py — &*) + sin a* 


Damit bei y = 165° Stetigkeit entsteht, muB 2 = 0,88 057 gewahlt werden. ; 

Gegen die Hinterkante des Profils muB erfahrungsgemaB der Druck p ansteigen, also nach der 
Bernoullischen Gleichung (46) die Geschwindigkeit fallen. Soll dies bei 80%, der Fligeltiefe be- 
ginnen, dann miissen wir 
von x= 0,6 ab zu k, y,(p) 
einen Zusatz k, y,(p) ad- 
dieren, der bei x = 0,6 
verschwindet, damit 1/q@ 
immer stetig ist. Den in 
x erst steilen, dann fla- 
chen Anstieg des Drucks, 
also den entsprechenden 

Abfall von v? beziehungs- 
a f weise jw? erreichen wir 
Abb. 2. Skizzen der aus dem Ansatz des oe entstehenden Verteilungen 1/a(¢), (7) durchieines linearent ne 
stieg von y, nach Abb. 2. 

Ubertragen wir ihn auf pm mittels (22), dann ergibt sich fiir y,(y~) der Ansatz 


v,(~) = cos y — cos q, , | 


~, = are cos 0,6. | 


(59) 


Damit ist der gesamte Ansatz erledigt. Er ist nur in der Nahe der Nase nicht ganz befriedigend, 
da dort der Wurzelfaktor in (20) ein im voraus nicht zu iibersehendes starkes Absinken der Ge- 
schwindigkeit v(m, a*) gegeniiber @(p) mit sich bringt. Dort nimmt man die Geschwindigkeit so, 
wie sie auf Grund des Ansatzes entsteht und ist zufrieden, wenn sie nicht iber den konstanten, 
durch k, y, (p) erzielten Wert ansteigt. Dies ist immer garantiert, da der Wurzelfaktor nur ein 
Absinken von v(p, a*) bewirkt. Der gleiche Mangel haftet allen zuriickliegenden Verfahren an. 
Er storte bisher nicht, da er auf eine kleine Umgebung der Profilnase beschrankt bleibt. 


Mit den Ansatzen (58) und (59) lauft nun das Verfahren ganz automatisch: 


1) Wir berechnen die Funktion sin (p— a*) + sin a* und damit h,(y) und h,(~) nach (32) an 
2 N aquidistanten Punkten zwischen =0 und = 180°. 


2) Wir fihren die harmonischen Analysen (32) durch und berechnen die Cy und: “e7s. = Win 
haben hier reine Sinus-Analysen, die wenig Aufwand erfordern. 


} R. Eppler, Z. Flugwiss. 3 (1955), S. 345. 
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3) Wir berechnen nach (38) s, und s,, nach (57) 4) und A, sowie nach (43) die a) und al), 

4) Wir fithren die Synthesen der Funktionen x((p), x()(p), yO), ¥(~), YO(p) und ¥3)(~) 
durch und berechnen aus den beiden letzten g()(y) und g@)(q). ; 

5) Wir wahlen ein k, und berechnen daraus zuerst nach (53) v(p, a*) und wenn wir damit 
zufrieden sind, auch x(q) und y(@). 

6) Weicht v*(y, a*) im konstanten Teil von 1/@ zu sehr von einer Konstanten ab, dan» be- 
rechnen wir mit einem konstanten Wert nach (56) noch ein y4(p), damit entsprechend 1) bis 4) 
Funktionen «(4)(), yO(p), und g“(q), mit ihnen wieder zuerst v*(~, a*). Dabei versuchen wir es 
zunichst an wenigen Stellen mit verschiedenen k,-Werten, entnehmen den giinstigsten und be- 
rechnen damit v*(p,a*), a(~) und (9). 

Mit k, kann das Dickenverhaltnis 


9 
— “Ymax 
i= i (60) 
des Profils geregelt werden. Dabei ist 
L= 1 — k, 10) (61) 
mit 
1) = x()(0) — x®(180°) , (62) 


wahrend y,,,, nicht yon vorn herein anzugeben ist, da die y‘'),, nicht an den gleichen Stellen p auf- 
treten, also auch die Stellen, die auf y,,, fiihren, von 6 abhangig sind. Wir bericksichtigen diese 
Abhangigkeit nicht, sondern nehmen statt (60) ein Dickenverhaltnis 
, e 2 Ve 
5 ee (63) 
bei welchem wir y, an einer festen, ungefahr beim Maximum y,,,, liegenden Stelle gy. entnehmen. 
Damit wird 6 nach (60) héchstens geringfiigig gréBer als 6’, wenn ymax, etwas groBer ist als y,. 
Der Unterschied ist aber unwesentlich, da y(y) in der Nahe des Maximums sehr flach verlauft. 
In (63) ist 
Yo = ¥(Pe) + hes ¥O(H) = ¥ + by yO), (64) 
und es folgt 
[ é/ Sz oD) Vy 
8 JOY + 2ZyP ° 
Ganz analog geht man vor, wenn man auch in 6) gleichzeitig mit der Wahl von k, auch noch k, 
so festlegen will, daB ein vorgegebenes Dickenverhaltnis entsteht. Hierfiir gilt dann 


(65) 


y__ 2(yD hg ye? + ha ye’) 
O = 1 — k, 13) = k, 1 (66) 
oder ; 
(1 — hk, I) 0’ —2 (y + ky ¥) 
ks = 1D § 2 yO - (67) 


Wir konnen auf diese Weise Profile mit allen gebrauchlichen Dickenverhaltnissen berechnen. 
Wir verfolgen hier aber nur den relativ hohen Wert 6’ = 0,21 niaher, bei dem der EinfluB der Kor- 
rektur gut zu erkennen ist. Weitere Ergebnisse wurden schon an anderer Stelle bekanntgegeben '. 
Die numerischen Rechnungen wurden mit einem yom Verfasser aufgestellten Rechenschema fir 
Analyse und Synthese mit N = 12 durchgefiihrt. Dies andert am Prinzip nichts. Man kann je nach 
der gewiinschten Genauigkeit ebensogut mit den Rechenformularen von Hufmann® mit N = 9 
oder NV = 18 arbeiten. 

In Tabelle 1 sind y,, y, und die zunachst daraus berechnetev Funktionen x), (3), 0), y(3), of0) 
und g@) angegeben. Fiir y, wurde die Stelle p, = 13+ 7,5° genommen, die nach (65) mit 6’ = 0,21 
auf k, = 1,27521 fiihrt. Dazu ergibt sich das v*(y, a*) der neunten Spalte, das in Abb. 3 tiber x(¢) 
nach (22) aufgetragen ist. Zum Vergleich ist auch m*(y) gestrichelt eingezeichnet. Die Berechnung 
von x(q) nach (45) ist noch nicht notwendig, wenn man nur den Einflu8 des Wurzelfaktors tiber- 
sehen will. Dieser macht bei unserem Profil so viel aus, da wenigstens in dem Bereich, wo die Ge- 
schwindigkeit konstant werden sollte, eine Korrektur angebracht ist. 

Die Korrektur ist in den letzten Spalten von Tabelle 1 zu verfolgen. Die Funktion y, ist nach (56) 
mit einem konstanten 0 zwischen gy = 8-7,5° und w = 18: 7,5° berechnet. Dieser Bereich um- 


1 R. Eppler, FuBnote 1 von S. 446. 
2 FuBnote 1 von S. 442. 
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Tabelle 1. Beispiel A, Zwischenergebnisse. 
ED era PRP emer me aan Mayr i) (0) y@) 
a a  ——————————_LL—— 

0,40 000 0,99 491 0,07 870 0,00 000 0,00 000 

‘ Se 0,39 144 0,98 629 0,08 069 0,00 005 0,00 015 

2 0,36 593 0,96 053 0,03 673 0,00 018 0,00 258 

3 0,32 388 0,91 819 0,09 478 0,00 047 0,00 905 

4 0,26 603 0,86 000 0,10 358 0,00 085 0,02 085 

5 0,19 335 0,78 699 0,10 946 0,00 139 0,03 868 

6 0,10 711 0,70 043 0,10 973 0,00 205 0,06 123 

7 0,00 876 0,60 182 0,09 999 0,00 292 0,08 470 

8 0,00 000 0,49 296 0,08 426 0,00 390 0,10 157 

9 0,37 567 0,06 759 0,00 500 0,11 424 

10 0,25 202 0,04 901 0,00 620 0,12 361 

11 0,12 415 0,03 006 0,00 749 0,12 921 

12 | — 0,00 568 0,01 129 0,00 881 0,13 197 

13 — 0,13 524 — 9,00 732 0,01 010 0,13 254 

14 — 0,26 228 — 0,02 631 0,01 133 0,13 040 

15 — 0,38 459 — 0,04 448 0,01 247 0,12 507 

16 — 0,50 003 — 0,06 116 0,01 340 0,11 736 

UG — 0,60 662 — 0,07 639 0,01 411 0,10 795 

18 — 0,70 243 — 0,09 079 0,01 446 0,09 650 

19 — 0,78 591 — 0,10 321 0,01 440 0,08 284 

20 — 0,85 544 — 0,11 324 0,01 381 0,06 772 

21 — 0,09 995 — 0,12 106 0,01 245 0,05 174 

22 1,00 000 — 0,94 827 — 0,12 681 0,00 995 0,03 492 

23 1,15 832 — 0,97 082 — 0,13 014 0,00 537 0,01 750 

24 0,00 000 0,00 000 — 0,97 853 — 0,13 125 0,00 000 0,00 000 

Fortsetzung. 
@|7,5° g() g3) v?(p, X*) Vs x64) yA gh) 

0 0,00 000 0,00 000 0,56 437 0,00 000 | —0,00 447 0,00 000 0,00 000 
1 0,00 612 0,04 159 0,57 259 — 0,00 454 | — 0,00 004 | — 0,00 364 
2 0,00 622 0,13 043 0,59 254 — 0,00 472 | — 0,00 010 | — 0,00 215 
3 0,00 643 0,17 479 0,63 473 — 0,00 505 | — 0,00 024 | — 0,00 447 
4 0,00 700 0,23 189 0,69 547 — 0,00 548 | — 0,00 050 | — 0,00 477 
5 0,00 774 0,25 502 0,79 961 — 0,00 605 | —0,00 095 | — 0,00 751 
6 0,00 831 0,26 631 0,96 332 — 0,00 668 | — 0,00 162 | — 0,00 864 
7 0,00 890 0,19 993 1,31 201 0,00 000 — 0,00 748 | — 0,00 276 | — 0,01 529 
8 0,00 932 0,12 467 1,42 362 0,00 741 — 0,00 788 | — 0,00 470 | —0,01 944 
9 0,00 972 0,09 518 1,44 657 0,01 405 — 0,00 734 | — 0,00 688 | — 0,01 738 
10 0,00 995 0,05 831 1,46 797 0,02 021 — 0,00 602 | — 0,00 881 | —0,01 384 
ll 0,01 018 0,03 217 1,47 769 0,02 300 — 0,00 413 | — 0,01 019 | —0,00 817 
12 0,01 016 0,01 178 1,48 203 0,02 424 | —0,00 202 | —0,01 100 | —0,00 446 
118} 0,01 004 | —0,00 411 1,48 340 0,02 461 0,00 023 | — 0,01 128 0,00 004 
14 0,00 955 | — 0,03 095 1,48 147 0,02 407 0,00 244 | — 0,01 100 0,00 399 
15 0,00 891 — 0,05 632 1,47 478 0,02 216 0,00 454 | — 0,01 023 0,00 919 
16 0,00 743 — 0,07 840 1,46 484 0,01 931 0,00 638 | — 0,00 896 0,01 339 
17 0,00 559 ; —0,10 201 1,45 010 0,01 507 0,00 788 | —0,00 716 0,02 127 
18 0,00 171 — 0,14 227 1,41 560 0,00 507 0,00 849 | —0,00 489 0,02 469 
19 — 0,00 364 | —0,19 102 1,36 013 0,00 000 0,00 810 | —0,00 311 0,01 781 
20 — 0,01 609 | —0,25 371 1,26 672 0,00 756 | — 0,00 204 0,01 367 
21 — 0,03 725 | —0,35 483 1,09 418 0,00 714 | — 0,00 134 0,01 309 
22 — 0,11 807 — 0,57 720 0,71 246 0,00 686 | — 0,00 081 0,01 564 
23 — 0,40 969 | —1,34 273. 0,18 495 0,00 669 | —0,00 039 0,03 022 
24 — 2,99 769 | — 8,57 134 0,00 546 0,00 000 0,00 665 | — 0,00 000 0,18 123 


faft, wie in Abb. 4 zu erkennen ist, den wesentlichen Teil des Profilbereichs, in dem die Verbesserung 
beabsichtigt ist. Nur eine kleine Umgebung der Nase ist nicht eingeschlossen. Dort ist das Ab- 
sinken der Geschwindigkeit ohnehin physikalisch sinnvoll. Die GréBe + wurde so gewahlt, da® Ve 
an den Grenzen des y-Bereichs ungefahr verschwindet, also die Stetigkeit von 1/a nicht gestért 
wird. Die Funktionen «(*), y() und g( ergeben sich aus y, wieder automatisch. An wenigen Stellen 
wurden dann mit verschiedenen Werten k, nach (53) v2(p, «*) berechnet, wobei jeweils k, nach (67) 
verwendet wurde. Der Wert k, = 0,855 erwies sich als der giinstigste. Er ergab die in Tabelle 2 
niedergelegten Ergebnisse. Die Profilkoordinaten wurden dabei noch auf den iiblichen Mafstab 
transformiert, bei dem die Nase bei X = 0 , die Hinterkante bei X — 100 liegt. In Abb. 4 sind die 
Ergebnisse dargestellt. Die nach (27) berechneten Linien v(y, a) sind fiir « = 0 und « = + 0,9° 
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Tabelle 2. Beispiel A, Endergebnisse. 


Q/7,5° x yy v2(p, 0*) @7,5° | ON Y v*(p, X*) 
0 100,00 0,00 0,5553, ils} 40,11 10,50 1,4679 
1 99,34 0,01 0,5639 14 34,03 10,42 1,4678 
2 97,36 0,21 0,5841 yas 28,17 10,10 1,4679 
3 94,25 0,73 0,6277 16 22,61 9,61 1,4679 
4 90,15 1,67 0,6902 17 17,47 9,00 1,4692 
5 85,42 3,09 0,7995 18 12,94 8,23 1,4657 
6 80,85 4,88 0,9728 19 9,04 7,24 1,4144 
cl 75,36 6,73 1,3588 20 Dadi 6,06 1,3051 
8 710,22 8,02 1,4644 E Pal 3,20 4,75 ES 
9 64,61 8,98 1,4074 |) ee? 1,42 3,30 0,7032 

10 Done 9,69 1,4687 eae 0,36 1,68 0,1768 
AGH 52,63 10,14 1,4682 24 0,00 0,00 0,0051 
i 46,37 10,39 1,4680 dedi =— Mel 1O 2 


mit eingezeichnet. Die letztere entspricht wegen der Symmetrie auch der Profilunterseite 
bei a = a* =—0,9°. Die Konstanz von v*(y,a*) ist auBer in der Nahe der Nase wirklich 
befriedigend. 

Mit der Lange / nach (61) ist mittels (15) die GréBe dc,/dax berechnet, die vielfach interessiert. 
Fiir die praktische Verwendung des Profils, insbesondere fiir die Wélbung mittels Wirbelverteilung 
werden oft noch einige andere Daten verlangt. Die Geschwindigkeitsinderung Av(y) bezogen auf 
eine Anderung des Auftriebsbeiwerts 4c, — 1 berechnet man iiber den Anstellwinkel az, bei dem 
C, = list. Er folgt leicht aus (15). Mit ihm ist Jv(y) = v(¢, az) — v(¢, 0). Darauf gehen wir hier 
jedoch nicht naher ein, da wir auch gewolbte Profile mit unserem Verfahren berechnen. Ebenso ist 
die Angabe des Nasenradius iiberfliissig, da bei 
N= 12 die Punkte in der Nahe der Nase sehr 


dicht liegen. 20 
| 7 =H 
=-09° Unterseite 
16 ~ | if =f 2. 
es 


a=—09° Oberseite 


Las 
0 50 700 
Abb. 3. Der Einflu8 des Wurzelfaktors beim Abb. 4. Das Profil des Beispiels A mit der korrigierten 
Beispiel A. Geschwindigkeitsverteilung. 


Es sei aber noch bemerkt, da® man bei den diinneren Profilen unserer Serie die Funktion y4 
nicht jedesmal neu berechnen muG. Man kann sie vielmehr samt den Funktionen «@), y und g) 
dem Beispiel 6’= 0,21 entnehmen und erhilt damit einen geniigend genauen Ausgleich. Man mu 
dabei auch die GréBe von k, nicht jedesmal neu ausprobiercn, sondern erhalt mit k, = 0,855 (6’/0,21)? 
gute Verbesserungen. 

B. Als zweites Beispiel berechnen wir eine Serie von unsymmetrischen Profilen, die bei einem 
Anstellwinkel «, = 5,6° auf der Oberseite, bei 0, = 7,4° auf der Unterseite von der Nase bis etwa 
60%, der Tiefe konstante Geschwindigkeit haben sollen. Der Druckanstieg auf der Hinterseite soll 
wieder eine abnliche Form wie beim Beispiel A haben. Diese Aufgabenstellung ist, abgesehen von 
der etwas kleineren Lange der Strecken konstanter Geschwindigkeit, das Analogon zum Beispiel A, 
nur sind die Anstellwinkel der konstanten Geschwindigkeit auf Ober- und Unterseite um + 6,5° 
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verschoben. Wir bekommen gewolbte Profile. Wir wahlen analog zum Beispiel A 


ye ene ake (0<pSm+2on), | 68 

alg) =f ee (1+ 20.59 S22), f ie 

at sin (p—2 4,,) 25 ee Oa 69 

”i() i sin (po) oe ain and ea l (69) 

VA) =41— Pos | (70) 
[7 5in (@@ —4,) + sin 2 (v2 <9 <2m). 


Die zweiten Terme in (69) und (70) vermeiden wieder, daB x(p) an den Stellen a + 2 oa, und 7 + 2 a5 
verschwindet und erzeugen bei gy =2-+ 2a, senkrechte Tangente. Der vordere Staupunkt liegt 
allgemein beig=2-+2a. Die Gré8e a, ist also derjenige Anstellwinkel, bei dem der Staupunkt 
an der Stelle senkrechter Profiltangente, der Nase, liegt. Wir wihlen deshalb a, =6,5° in der Mitte 
zwischen a, und a, weil bei diesem Anstellwinkel ohnehin an der Nase ungefahr symmetrischer 
Strémungscharakter (bzw. der,,stoBfreie Kintritt) zu erwarten ist. AuBerdem wablen wir j= 0,87389, 
dann wird die Ubergangsstelle q,, an welcher die beiden Ausdriicke fiiry, die gleichen Werte bekommen, 
Qo = 210°, wihrend die entsprechende Stelle p, von y, etwas groBer als 180° wird. Auf den genauen 
Wert kommt es nicht an, da wir y,(p) nur an den Stellen 180° und 187,5° beniitzen. Zwischen diesen 
beiden Werten liegt dann der Ubergang von der ersten zur zweiten Darstellung von y,(@). 

Den Druckanstieg wahlen wir ebenfalls analog zum Beispiel A. Wir haben dort bemerkt, daf 
x(~) im hinteren Bereich des Profils sich gegeniiber (22) verkleinert. Wir wahlen deshalb 


ys(~) = cosp—cosps, cosy, = 0,3, (0 Sy Ss)» (71) 
und nicht cos gm; = 0,2, wie es an sich entsprechend der Aufgabenstellung nach (22) sein miuBte. 
Auf der Unterseite wahlen wir 
sin (py + 0) — sin a, 
sin (p —a,) + sina, 


Va(P) = (cos p — cos @) (27— 93 Sp S2n), (72) 
an allen nicht erfaften Stellen soll y,(~) verschwinden. Damit wird h,(g) nach (32) mit «*(p) nach (68) 


eine ungerade Funktion, wodurch an der numerischen Rechnung einiges gespart wird. 


Die Stetigkeit von h(p) nach (24) ist gesichert, weil die Sprungstellen von a*(y) jeweils dort 
hegen, wo h(~) verschwindet. Die Beschrinkung von y,(q) auf die Oberseite und von y,(p) auf die 
Unterseite ist vorteilhaft, weil nach (25) wegen der positiven Werte a*(y) ein Geschwindigkeits- 
unterschied zwischen Ober- und Unterseite entstehen mu. Dieser ist besonders gut zu tibersehen, 
wenn fiir Ober- und Unterseite zwei getrennte Parameter k, und k, zur Verfiigung stehen. Das Vor- 


zeichen der Funktionen y,(p) ist, wie man nach (12) leicht nachpriift, so gewahlt, daB alle k; positiv 
werden. 


Die Berechnung der Funktionen x), y und g( lauft genau wie beim vorangehenden Beispiel. 
Wir verzichten hier auf die tabellarische Wiedergabe der Rechnungen. In Abb. 5 ist das entstehende 
Profil x(~), y(@) fiir kj = 0 und k, = 1 dargestellt. Einige Punkte gleicher y-Werte sind mitein- 
ander durch Geraden verbunden. Man erkennt so, auf welche Weise die Profile fiir Zwischenwerte 
von k, entstehen: Man braucht wegen (45) nur auf jeder Verbindungsgeraden k, im Verhaltnis zur 
ganzen Lange abzutragen. Fir k; = 0,5 ist diese Konstruktion durchgefiihrt. Ebenso kénnen wir 
fir alle anderen Werte k, Profile berechnen. Wir greifen wieder einen festen heraus, namlich 
k, = 0,45 037, der auf eine Profildicke 6 < 0,15 fiihrt. Die hiermit berechneten Koordinaten X(qg), 
Y(), die Geschwindigkeiten »?(~, 0*) sowie die nach (15) und (17) berechneten GréBen de,/dax und Cn 
sind in Tabelle 3 angegeben. Die Koordinaten sind dabei so transformiert, daB das Profil von 
X = 0 bis X = 100 geht und ungefahr waagrecht liegt, was durch eine Drehung um 6,5° erreicht 
wurde. In Abb. 6 sind die Ergebnisse dargestellt. Zu v*(p, a*) wurden erginzend nach (27) die 
gesamten Verteilungen v*(~p, a,) und v°(p, 45) berechnet und mit eingezeichnet. Unter dem Diagramm 
der Abb. 6 ist das Profil noch in der urspriinglichen Lage dargestellt, auf welche sich das y’ des 
Wurzelfaktors bezieht. Die gestrichelten Linien geben w*(g) in den Teilen, in denen diese Funktion 
konstant ist. Die entstehende Verteilung v?(~, a*), also auf der Oberseite v>(~, a,), auf der Unter- 


seite v*(p, ay) ist an den Stellen waagrechter Tangente (des nicht gedrehten Profils) identisch mit 
«*(q), an allen anderen Stellen kleiner. 
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Der Einflu8 des Wurzelfaktors ist erheblich. Vor allem riickt das (in der Grenzschichttheorie 
entscheidende) Geschwindigkeitsmaximum auf der Oberseite nach hinten, auf der Unterseite nach 
vorn. Dieser unterschiedliche Einflu8 wurde bisher noch nirgends erfaBt, da sowohl bei Theodorsen 1 
als auch bei Truckenbrodt ? der Wurzelfaktor nur am ungewélbten Profil beriicksichtigt wird. Vor 


Abb. 5. Die beim Beispiel B fir k, = 0, k, = 0,5 und k, =1 entstehenden Profile x(q), y(g). 


(i) 


Abb. 6. Ein spezielles Profil der Serie B samt seinen Geschwindigkeitsverteilungen, 


allem bei starker gewélbten Profilen ist es also wichtig, die zum Profil gehdrenden exakten Ge- 
schwindigkeitsverteilungen zu kennen. Eine Anderung derselben wurde bei diesem Beispiel nicht 
durchgefiihrt, da die Kinfliisse des Wurzelfaktors weitgehend vorauszuschen waren und die ent- 
stehenden Verteilungen auch ohne Korrektur physikalisch giinstig sind. 

1 FuBnote 4 von S. 436. 

2 FuBnote 6 von S. 436. 
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Es sei abschlieBend noch bemerkt, da® keines der beiden behandelten Beispiele mit den bisher 
vorliegenden Methoden hatte direkt angegangen werden kénnen. Dies wurde erst durch die hier 


eingefiihrte Wahl a*(¢) ermoglicht. 


Tabelle 3. Beispiel B, Endergebnisse. 


——_—_——— eee. OOD 
g/7,5° | x x¢ v?(p, a*) | 9[7,5° x yy v?(p, x*) 
0 100,00 0,00 0,0000 26 0,02 — 1,98 0,0006 
1 99,45 0,19 0,9034 PAI 0,56 — 3,28 0,5911 
2 97,94 0.67 0.9452 28 1,93 Sai 0.9275 
3 95,32 1,46 0,9795 29 4,13 — 5,08 0,9507 
4 91,86 Ae yi 1,0331 30 7,24 — 5,67 0,9472 
5 87,72 4,01 1.1019 31 11.06 — 6,07 0,9422 
6 82,97 5,76 1,1975 32 l>you — 6,33 0,9379 
a 77,91 dele Sao = 333) 20,77 — 6,44 0,9345 
8 72,70 9,80 1,4961 || 34 26,40 — 6,39 0,9307 
9 67,41 11,83 1,7304 35 32,44 — 6,20 0,9286 
10 62,29 13,42 1,9492 36 38,78 SD 0,9243 
11 97,02 14,43 1,9593 37 45,18 — 5,31 0,9202 
12 51,48 15,03 1,9529 38 51,62 — 4,53 0,9074 
13 45,95 15,26 1,9373 39 98,18 — 3,33 0,8251 
14 40,36 15315 1,9165 40 64,91 730) 0,7434 
15) 34,83 14,73 1,8865 41 71,76 — 0,90 0,6856 
16 29,50 14,01 1,8507 42 78,36 - — 0,04 0,6432 
17 24,34 13,04 1,8051 | 43 84,39 0,51 0,6119 
18 19,52 11,82 1,7488 | 44 89,74 0,73 0,5888 
19 15,16 10,38 1,6766 | 45 94,10 0,68 0,5724 
20 1h Au 8,78 1,5833 || 46 97,30 | 0,46 0,5603 
Pa 71,84 7,04 1,4531 | 47 99,34 0,16 0,5477 
22 5,07 5,20 1,2710 | 48 100,00 0,00 | 0,0000 
23 2,83 3,34 0,9990 
24 1.29 1,46 0.5785 deq/dx = 1,120 - 2x 
25 0,38 — 0,38 0,1025 Cmo = — 0,1848 


6. Zusammenfassung und Ausblick. Es wurde ein einfaches Verfahren angegeben, das die Be- 
rechnung von Profilen und Profilserien mit vorausbestimmten Eigenschaften der Geschwindigkeits- 
verteilung sehr allgemein erméglicht. Es erlaubt, die Geschwindigkeitsverteilung stiickweise bei 
verschiedenem Anstellwinkel einzufiihren, die Profildicke zu regeln, und es kann ohne Schwierig- 
keiten auf unsymmetrische Profile angewandt werden. Dabei zeigt es sich, da die bisher be- 
kannten Naherungsverfahren einen wesentlichen EinfluB nicht bericksichtigen. 


Das Rechenverfahren liefert zu den berechneten Profilen ihre exakte Geschwindigkeitsverteilung, 
diese stimmt nur nicht exakt mit der Vorgabe iiberein, sondern ist ihr gegeniiber um einen Faktor 


1//1 + y'? = cos O verringert, wo @ der von vorn herein nicht bekannte Steigungswinkel des 
Profils ist. Der EinfluG des cos-Faktors ist fast iberall so klein, daB man die Abweichung in Kauf 
nehmen kann oder sie durch einen einzigen Korrekturschritt beseitigen kann. Nur in der Umgebung 
der Profilnase ist der Einflu8 zu stark. Dort belaBt man, wie bei den friiheren Arbeiten, das Ab- 
sinken der Geschwindigkeit zum Staupunkt hin so, wie es der cos-Faktor ergibt. Dieser letzte 
Mangel wird in einer weiteren Arbeit des Verfassers noch beseitigt werden, die auf anderem Weg 
auch sonst einige tiefergehende Ergebnisse zeitigt. Das dort entwickelte Verfahren ist zwar im 
Aufwand des einzelnen Profils nicht ungiinstiger als das vorliegende, es erlaubt aber keine so ein- 
fache serienweise Regelung der Profildicke mehr. Solange man es also mit serienweisen, nicht sehr 
dicken oder sehr stark gewélbten Profilen zu tun hat, deren Geschwindigkeitsverteilungen nicht 
sehr genau mit den Vorgaben tibereinstimmen miissen, wird das vorliegende Verfahren vorteilhafter 


bleiben. 


(Eingegangen am 17. Mai 1955) 


Anschrift des Verfassers: Dr. R. Eppler, Stuttgart-W, Leibnizstr. 84. 
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